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Semiconductor  Measurement  Technology  : 
EVALUATING  A  CHIP,  WAFER,  OR  LOT 
USING  SUXES,  SPICE,  AND  STAT2 

Janet  C.  Marshall  and  Richard  L.  Mattis 
Semiconductor  Electronics  Division 
National  Institute  of  Standards  and  Technology 
Gaithersburg,  MD  20899 

Abstract 

The  computer  procedure  KEYS  (linKing  softwarE  to  analYze  waferS)  links  SUXES  (Stan- 
ford University  eXtractor  of  modEl  parameters),  SPICE  (a  Simulation  Program  with  In- 
tegrated Circuit  Emphasis),  and  STAT2.  Given  data  points  for  individual  devices,  SUXES 
obtains  the  model  parameters  for  SPICE.  SPICE  predicts  the  behavior  of  an  individual 
device  or  an  entire  circuit.  After  analyzing  each  test  chip  on  a  wafer,  STAT2  determines 
correlation  coefficients  and  generates  wafer  maps  of  selected  parameters.  These  wafer  maps 
are  valuable  to  the  designer,  modeler,  and  process  engineer. 

The  entire  package  accomplishes  the  following:  (1)  standardizes  the  technique  of  running 
SUXES  and  SPICE  in  an  integrated  mode;  (2)  simulates  and  plots  the  characteristic 
curves;  (3)  simulates  and  plots  the  results  of  an  optional  dynamic  circuit  (for  example, 
a  ring  oscillator);  (4)  performs  steps  (2)  and  (3)  for  every  test  chip  on  each  wafer;  (5) 
summarizes  the  results  from  each  chip,  each  wafer,  and  the  lot;  (6)  rank-orders  the  model 
parameters  for  each  wafer  according  to  their  correlation  coefficients  with  respect  to  chosen 
circuit  parameters;  and  (7)  generates  wafer  maps  of  several  quantities.  A  CMOS  19-stage 
ring  oscillator  is  used  to  illustrate  the  capabilities  of  KEYS. 

Key  words:  circuit  simulator;  CMOS;  ring  oscillator;  SPICE;  STAT2;  SUXES;  wafer  maps 
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NOMENCLATURE 


Symbol 

Units 

Identification 

F/m 

Permittivity  of  free  space 
8.85  x  10-12  F/m 

Cox 

F/m 

Static  dielectric  constant  of  SiC>2 
3.9  x  ec  =  34.5  x  10~12  F/m 

F/m 

Static  dielectric  constant  of  silicon 
11.7  x  eQ  =  103.5  x  10  "  F/m 

ERR 

% 

The  average  error  is  the  sum  of  the  absolute 
value  of  the  normal  errors  divided  by  the 
number  of  points.  The  normal  error  is  the 
percent  difference  between  the  measured  data 
points  and  SUXES's  predicted  data  points. 

Ids 

A 

Drain-to-source  current 

k 

J/(K) 

Boltzmann's  constant  1.38  x  10-23  J/(K) 

m 

Designed  channel  length  for  an  n-channel 
MOSFET 

m 

Designed  channel  length  for  a  p-channel 

A  if /TIT"  rp 

A  /TIT™ 
MxlZ 

The  measured  circuit  response,  e.g.,  frequency 

PDIFF 

% 

The  percent  difference  between  the  measured  and 
simulated  circuit  response  is  the  absolute 
value  of  [(SCKT  -  MCKT)/(MCKT)]  x  100% 

WIIClC  11   IS   do  a  ULlllv^UL    L 11 CL  L    b  11                      Ul           Lll  V/  Ul  b 

response  is  the  reference  response. 

q 

c 

Electronic  charge  =  1.602  x  10~19  C 

RMS 

V 

The  root  mean  square  deviation  between  the 

lil             111           1 i                                              1    "  j        1  ■ 

measured  threshold  voltage  curve  and  its  linear 
least-squares  fit  when  plotted  against 
\/|Vbs|  +  tphif  -  ^tphif  where  tphif=0.6  V. 

SCKT 

MHz 

The  simulated  circuit  response,  e.g.,  frequency 

T 

IV 

j.emperaiure 

Vb8 

V 

Body-to-source  voltage 

T  T 

Vds 

■f  T 

V 

Drain-to-source  voltage 

v 

V  gs 

V 

Gate-to-source  voltage 

wn 

m 

Designed  channel  width  for  an  n-channel 
MOSFET 

Wp 

m 

Designed  channel  width  for  a  p-channel 
MOSFET 
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1.  INTRODUCTION 


This  report  describes  a  computer  procedure  called  KEYS  (linKing  softwarE  to  analYze 
wafer S)  that  was  developed  at  the  National  Institute  of  Standards  and  Technology.  KEYS 
links  SUXES  (Stanford  University  eXtractor  of  modEl  parameters)  [1],*  SPICE  (a  Simula- 
tion Program  with  Integrated  Circuit  Emphasis)  [2],  and  STAT2  [3]  to  facilitate  integrated 
circuit  evaluation.  Other  examples  of  packages  for  integrated  circuit  evaluation  have  been 
described  [4,5].  A  simplified  block  diagram  for  KEYS  is  given  in  figure  1. 

SPICE  predicts  currents  and  voltages  at  chosen  nodes  in  a  circuit  and  requires  a  models 
file  which  contains  input  parameters  for  the  devices  in  the  circuit.  The  accuracy  of  these 
input  parameters  and  the  device  models  determines  how  well  the  simulated  currents  and 
voltages  predict  the  measured  currents  and  voltages.  Since  developing  a  models  file  for 
SPICE  is  a  difficult  task,  a  computer  program  called  SUXES  was  developed.  SUXES 
obtains  most  of  these  model  parameters  and  adapts  to  any  model.  Given  initial  estimates 
for  the  parameter  values,  upper  and  lower  bounds  for  these  values,  and  measured  data 
points  for  each  device,  SUXES  optimizes  the  parameter  values  by  varying  them  slightly, 
according  to  an  initial  estimates  file,  an  options  file,  and  a  strategy  file,  so  as  to  decrease  the 
root  mean  square  (rms)  deviation  between  the  predicted  and  measured  data  points.  To  aid 
process  engineers,  modelers,  and  designers,  STAT2  determines  correlation  coefficients  with 
respect  to  selected  parameters  and  generates  wafer  maps  of  selected  parameters.  These 
maps  provide  a  visual  inspection  of  the  positional  variations  of  parameter  values  over  a 
wafer. 

The  flowchart  for  the  computer  procedure  KEYS  is  given  in  figure  2.  To  perform  these 
tasks,  additional  FORTRAN  programs  as  given  in  Appendix  A  were  written.  Command 
files  connect  these  programs  so  that  they  may  run  as  a  single  entity.  The  software  used  in 
KEYS  is  listed  in  table  1. 

In  the  computer  procedure  KEYS,  there  are  three  sets  of  parameters:  (1)  the  estimated 
input  parameters  for  SUXES,  (2)  SUXES's  optimized  parameters,  and  (3)  the  model  input 
parameters  for  SPICE.  Given  the  data  points  for  individual  devices,  SUXES  obtains  the 
model  parameters  for  SPICE.  The  combination  of  SUXES  and  SPICE  gives  a  procedure 
for  parameterizing  experimental  data  from  circuits.  SUXES's  optimized  parameters  may 
not  be  physically  meaningful.  Nonetheless,  in  this  report,  the  assumption  is  made  that 
SUXES's  optimized  parameters  have  the  same  definition  as  the  model  input  parameters 
for  SPICE.  Hence,  the  software  to  be  presented  directly  inserts  these  parameters  into  the 
models  file  of  SPICE. 

The  computer  procedure  KEYS  characterizes  a  chip,  wafer,  or  lot.  In  this  report,  a  19-stage 


*  Certain  commercial  equipment,  instruments,  and  computer  programs  are  identified  in  this 
paper  in  order  to  specify  the  experimental  procedure  adequately.  Such  identification  does 
not  imply  recommendation  or  endorsement  by  the  National  Institute  of  Standards  and 
Technology,  nor  does  it  imply  that  the  equipment  or  program  identified  are  necessarily  the 
best  available  for  the  purpose. 
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4 


(Next  wafer) 


Figure  2.    Flowchart  for  the  computer  procedure  KEYS. 


5 


Table  1.       A  Description  of  the  Programs  Used  in  the  Computer  Procedure  KEYS 


# OF  SOURCE 

PROGRAM     PROGRAM  PURPOSE  AUTHOR 

NAME  STATEMENTS 


SUXES 
(version  I) 


3792 


Obtains  the  model  parameters 
for  SPICE 


SPICE  43414 
(version  2G.5) 


NCAR  Graphics  1764 
MCVAX 
UTILITY 


STAT2  16212 
(version  1.31) 
CRDB 
CMLIB 


Misc.  Files 
(version  I) 


1620 


The  circuit  simulator  which 
obtains  the  I-V  characteristics 
and  simulates  a  dynamic  circuit 


Makes  plots  on  the  specified 
output  devices 


Determines  the  correlation 
coefficients  and  generates 
wafer  maps 


These  include  the  command 
procedures,  the  editor  files, 
the  data  files,  and  the 
FORTRAN  files  needed  to 
link  the  software 


K.  Doganis 
Electrical  Engineering 

Software 
Santa  Clara,  Calif. 
(408)  296-8151 

Univ.  of  Calif. 
Berkeley,  Calif, 
with  revisions  by 

R.  Mattis 
(301)  975-2235 

NCAR  Graphics 

Information 
NCAR,  Boulder,  Colo. 
(303)  497-1201 

R.  Mattis 

Electronics  and  Electrical 
Engineering  Laboratory 
NIST,  Gaithersburg,  Md. 
(301)  975-2235 

J.  Marshall,  R.  Mattis, 

and  M.  Gaitan 
Electronics  and  Electrical 

Engineering  Laboratory 
NIST,  Gaithersburg,  Md. 
(301)  975-2052 
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ring  oscillator  is  used  to  illustrate  the  procedure.  This  example  takes  approximately  11 
CPU  minutes  per  chip  on  a  VAX-11/780.  The  entire  package  accomplishes  the  following: 
(1)  standardizes  the  technique  of  running  SUXES  and  SPICE  in  an  integrated  mode;  (2) 
simulates  and  plots  the  characteristic  curves;  (3)  simulates  and  plots  the  waveforms  of  an 
optional  dynamic  circuit  (for  example,  a  ring  oscillator);  (4)  performs  steps  (2)  and  (3) 
for  every  test  chip  on  each  wafer  if  desired;  (5)  summarizes  the  results  from  each  chip, 
each  wafer,  and  the  lot;  (6)  rank-orders  the  model  parameters  for  each  wafer  according  to 
their  correlation  coefficients  with  respect  to  chosen  circuit  parameters;  and  (7)  generates 
wafer  maps  of  selected  parameters.  The  wafer  maps  are  typically  made  for  the  most  highly 
correlated  parameters.  Therefore,  a  careful  selection  of  the  chosen  parameters  with  regard 
to  the  correlation  coefficients  will  improve  the  usefulness  of  the  wafer  maps  to  process 
engineers,  modelers,  and  designers. 

This  report  contains  three  parts:  Motivation  and  General  Principles,  An  Example  Using 
CMOS  Technology,  and  the  Conclusions.  In  the  first  part,  Section  2.1  describes  the  linking 
of  SUXES,  SPICE,  and  STAT2.  It  is  assumed  that  the  reader  is  familiar  with  SUXES, 
SPICE,  and  STAT2,  so  these  programs  are  only  briefly  described  in  sections  2.2,  2.3,  and 
2.4.  In  the  second  part,  section  3.1  describes  the  test  structure  measurement  methods,  and 
section  3.2  discusses  the  output  which  consists  of  plots,  tabular  listings,  and  wafer  maps. 
Appendix  A  contains  the  listing  of  the  computer  procedure  KEYS.  Appendix  B  presents 
a  method  of  modifying  SUXES 's  strategy  file  to  decrease  the  rms  deviation  between  the 
measured  data  points  and  SUXES 's  predicted  data  points  while  still  providing  realistic 
curves.  Finally,  Appendix  C  describes  the  procedure  to  run  KEYS. 

2.   MOTIVATION  AND  GENERAL  PRINCIPLES 

2.1  Linking  SUXES,  SPICE,  and  STAT 2 

SUXES,  SPICE,  and  STAT 2  are  the  programs  described  in  the  following  sections.  SUXES 
supplies  model  parameters  for  SPICE,  SPICE  simulates  circuits,  and  STAT2  obtains  cor- 
relation coefficients  and  generates  wafer  maps.  KEYS  can  be  used  to  provide  a  smooth 
transition  between  these  programs.  The  following  three  paragraphs  demonstrate  this  and 
the  last  paragraph  describes  how  it  is  done. 

Typically,  the  estimated  input  parameters  for  SUXES  are  manually  inserted  into  SUXES 's 
initial  estimates  file.  The  estimated  input  parameters  vary,  for  example,  with  oxide  thick- 
ness. Hence,  these  values  need  to  be  calculated  for  each  value  of  oxide  thickness.  KEYS 
can  be  used  to  perform  these  calculations  and  eliminate  manual  intervention. 

Typically,  the  model  input  parameters  for  SPICE  are  manually  inserted  into  the  models 
file  of  SPICE.  This  becomes  a  task  which  is  burdensome  since  each  chip  has  its  own  set  of 
parameters.  KEYS  can  be  used  to  insert  directly  SUXES's  optimized  parameters  into  the 
models  file  of  SPICE.  The  values  of  the  model  input  parameters  for  SPICE  are  also  stored 
for  future  correlation  coefficient  analyses  and  possible  wafer  maps. 

Typically,  the  values  for  the  model  input  parameters  for  SPICE  are  manually  inserted  into 
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files  suitable  for  STAT2.  KEYS  can  be  used  to  place  these  values  into  the  appropriate  files 
to  eliminate  manual  intervention.  Thus,  analyzing  chips,  wafers,  and  lots  has  become  an 
easy  task. 

In  smoothing  the  transitions  between  SUXES,  SPICE,  and  STAT2,  KEYS  uses  the  EDT 
editor  (copyright,  Digital  Equipment  Corporation).  Some  editors  can  be  imbedded  in  au- 
tomated programs  and  can  actually  edit  a  file  while  the  program  is  executing.  EDT  is  an 
example  of  such  an  editor.  This  eliminates  human  intervention  and  error.  A  command 
procedure  calls  the  file  to  be  edited,  and  uses  the  editor  commands  in  a  specified  file. 
Thus,  it  places  the  parameter  values  from  SUXES  into  SPICE,  and  it  modifies  command 
procedures  to  reference  the  data  files  for  the  next  chip.  The  linking  of  the  three  soft- 
ware packages  is  accomplished  in  this  manner.  The  procedure  to  run  KEYS  is  given  in 
Appendix  C. 

2.2  SUXES 

SUXES  [1]  obtains  the  model  parameters  for  SPICE  [2,6,7].  A  model-independent  Leven- 
berg-Marquardt  algorithm  is  used  to  minimize  the  rms  deviation  between  the  measured 
and  predicted  data  points.  This  minimization  makes  it  possible  to  perform  a  SPICE 
analysis  more  accurately.  Given  initial  estimates  for  the  parameter  values,  upper  and 
lower  bounds  for  these  values,  and  measured  data  points  for  each  device,  the  parameter 
values  are  optimized  by  varying  them  slightly,  according  to  an  initial  estimates  file,  an 
options  file,  and  a  strategy  file,  so  as  to  decrease  the  rms  deviation  between  the  predicted 
and  measured  data  points.  In  the  computer  procedure  KEYS,  the  parameters  optimized 
by  SUXES  are  indicated  in  table  2. 

In  KEYS  the  following  input  files  are  needed  to  run  SUXES:  the  data  file(s),  the  initial 
estimates  file,  the  options  file,  and  the  strategy  file.  After  the  optimization,  an  extracted 
parameters  file  and  a  plot  file  are  created.  The  extracted  parameters  file  contains  the 
optimized  values  of  the  parameters,  and  the  plot  file  contains  the  predicted  data  points. 

The  data  file(s)  contain(s)  the  measured  data  points.  The  proper  polarities  are  used  for 
the  devices.  There  is  an  exception:  for  j?-channel  MOSFETs,  a  positive  current  must 
be  recorded,  or  the  device  must  be  treated  as  an  n-channel  MOSFET.  The  automated 
procedure  ensures  that  the  proper  polarities  are  used  in  SUXES. 

The  initial  estimates  file  includes  the  initial  guesses  for  the  parameter  values,  the  lower 
bounds,  and  the  upper  bounds  for  these  parameters.  The  parameter  values  are  free  to 
vary  within  these  specified  boundaries.  The  upper  and  lower  bounds  used  in  this  analysis 
were  obtained  using  the  technique  described  in  the  SUXES  report  [1]  or  from  the  method 
described  in  Appendix  B  to  find  an  appropriate  strategy  file.  No  attempt  was  made  to 
quantify  the  values  of  these  bounds  (say  ±10%  or  ±<r  about  the  initial  estimate),  since 
each  parameter  has  its  own  optimization  sensitivity.  The  user  may  alter  these  values 
to  increase  the  number  of  acceptable  chips.  In  KEYS,  an  acceptable  chip  is  one  whose 
optimized  SUXES  values  are  not  pegged  at  the  lower  or  upper  bounds.  In  the  example  with 
CMOS  technology,  there  are  two  initial  estimates  files:  NINPUT.DAT  is  for  the  n-channel 
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Table  2.       A  Brief  Description  of  the  Parameters  Considered  in  this  Analysis 


PARAMETER 

UNITS 

PARAMETER  DESCRIPTION 

LEVEL 

— 

model  index 

TYPE 

— 

type  of  device: 

+1    n-channel  device 
—  1    p- channel  device 

*LD 

m 

lateral  diffusion 

♦  TOX 

m 

oxide  thickness 

*NSUB 

cm-3 

substrate  doping 

*VTO 

V 

zero-bias  threshold  voltage 

♦  KP 

A/V2 

y0.5 

intrinsic  transconductance  parameter 

♦GAMMA 

bulk  threshold  parameter 

♦  PHI 

V 

surface  potential  at  strong  inversion 

♦  UO 

cm2/V-s 

surface  mobility  at  low  gate  voltages 

♦  DELTA 



width  effect  on  threshold  voltage 

♦  VMAX 

m/s 

maximum  drift  velocity  of  carriers 

♦  XJ 

m 

metallurgical  junction  depth 

♦KAPPA 

saturation  field  factor 

♦  NFS 

cm-2 

effective  fast  surface  state  density 

♦  ETA 

— 

static  feedback  effect  parameter 

♦THETA 

V-i 

empirical  mobility  modulation  parameter 

CJSW 

F/m 

zero-bias  bulk  junction  sidewall  capacitance  per  meter 

of  junction  perimeter 

CGDO 

F/m 

gate-drain  overlap  capacitance  per  meter  channel  width 

CGSO 

F/m 

gate-source  overlap  capacitance  per  meter  channel  width 

CGBO 

F/m 

gate-bulk  overlap  capacitance  per  meter  channel  length 

RSH 

n/n 

sheet  resistance  of  the  drain  and  source  diffusions 

JS 

A/m2 

bulk  junction  saturation  current  per  square  meter  of 
junction  area 

TPG 

type  of  gate  material: 

MJ 

MJSW 

PB 

XQC 

LAMBDA 

UEXP 

UCRIT 

NEFF 

NSS 


V/m 
cm-2 


+1    opposite  to  substrate 
—1    same  as  substrate 
0    Al  gate 
bulk  junction  bottom  grading  coefficient 
bulk  junction  sidewall  grading  coefficient 
bulk  junction  potential 

thin-oxide  capacitance  model  flag  and  coefficient  of 

channel  charge  share  attributed  to  drain  (0.0  —  0.5) 

channel  length  modulation  parameter 

critical  field  exponent 

critical  field  for  mobility  degradation 

total  channel  charge  coefficient 

effective  surface  charge  density 


♦  The  parameters  found  via  SUXES,  in  this  analysis. 
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MOSFETs,  and  PINPUT.DAT  shown  in  table  3  is  for  the  p-channel  MOSFETs. 


This  list  contains  all  the  parameters  that  SUXES  can  optimize.  In  this  analysis  a  subset 
of  these  parameters  is  used.  The  methods  used  to  obtain  the  initial  estimates  for  this 
subset  of  parameters  are  given  in  table  4.  The  other  parameter  values  are  initialized  to 
the  default  values  specified  in  this  table. 

The  options  file  selects  the  desired  data  points  from  the  data  file(s),  determines  any  weight- 
ing of  the  data,  and  determines  the  method  of  convergence  or  termination.  In  the  example 
with  CMOS  technology,  all  the  data  points  are  used  in  SUXES.  The  options  file  for  the 
n-channel  MOSFETs  is  called  NOPT.FIL,  and  the  options  file  for  the  p-channel  MOS- 
FETs, shown  in  table  5,  is  called  POPT.FIL.  This  table  also  provides  the  definitions  of  the 
convergence/termination  criteria  used  in  this  file.  Additional  information  can  be  found  in 
the  original  report  on  SUXES  [1]. 

The  strategy  file  selects  the  order  and  groups  of  parameter  values  to  be  optimized.  In  the 
example  with  CMOS  technology,  the  parameters  whose  values  are  optimized  are:  LD,  TOX, 
NSUB,  VTO,  KP,  GAMMA,  PHI,  UO,  DELTA,  VMAX,  XJ,  KAPPA,  NFS,  ETA,  and 
THETA.  Refer  to  table  2  for  the  parameter  descriptions.  The  strategy  file  for  the  n-channel 
MOSFETs  is  called  NSTRAT.FIL,  and  the  strategy  file  for  the  p-channel  MOSFETs,  shown 
in  table  6,  is  called  PSTRAT.FIL.  A  parameter's  value  is  optimized  if  a  "1"  appears  in  the 
column  corresponding  to  the  iteration  of  interest.  A  procedure  is  described  in  Appendix  B 
which  modifies  the  strategy  file  in  order  to  decrease  the  rms  deviation  between  the  measured 
data  points  and  SUXES's  predicted  data  points.  This  procedure  provides  curves  that  are 
consistent  with  the  measured  data. 

After  the  optimization  is  complete,  SUXES  generates  a  plot  file.  In  the  example  with 
CMOS  technology,  the  Ids  values  are  predicted  for  all  combinations  of  Vgs  and  Vas  values. 
Therefore,  if  the  voltages  are  measured,  a  tremendous  number  of  I-V  curves  are  generated. 
To  save  memory,  to  decrease  the  computer  computation  time,  and  to  have  a  plot  file  of 
a  manageable  size,  the  measured  voltage  values  in  the  data  file(s)  should  be  rounded  or 
truncated  before  running  SUXES.  The  automated  procedure  does  this. 

The  capacitive  parameters  (CGDO,  CGSO,  and  CGBO)  are  calculated  after  the  opti- 
mization. The  equations  given  in  table  7  are  used  with  the  appropriate  values  from  the 
optimization. 

2.3  SPICE 

SPICE  performs  dc,  transient,  or  ac  small-signal  analyses  of  circuits  [2,6,7],  Dynamic 
circuits  are  typically  simulated  using  a  transient  analysis  which  predicts  branch  currents 
and  node  voltages  as  a  function  of  time.  A  dc  analysis  determines  the  circuit's  operating 
point,  making  it  possible  to  obtain  the  characteristic  curves.  The  ac  small-signal  analysis 
determines  the  ac  variables  as  a  function  of  frequency.  In  the  example  given,  the  I-V 
characteristics  for  the  p-channel  MOSFETs  and  the  response  of  the  optional  dynamic 
circuit  (the  ring  oscillator)  were  simulated  using  the  SPICE  files  listed  in  tables  8  and  9, 
respectively.  The  files  used  to  get  the  I-V  characteristics  for  the  n-channel  MOSFETs  are 
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Table  3.      A  P-Channel  MOSFET  Initial  Estimates  File  of  SUXES 


*  FILENAME: 


PINPUT.DAT 


*    parameters  initial  lowerbound  upperbound 


LEVEL 

3 

1 

3 

TYPE 

-1.0 

-1.0 

1.0 

LD 

0.45 

0.1 

2.0 

TOX 

635.0 

250.0 

1000.0 

NSUB 

1.2  x  1016 

1.0  x  1013 

1.0  x  1017 

VTO 

-0.85 

-3.0 

-0.1 

KP 

14.02  x  10"6 

1.0  x  10~6 

100.  x  10~( 

GAMMA 

1.16 

0.1 

2.0 

PHI 

0.3618 

0.3 

0.9 

UO 

500.0 

0.1 

1250.0 

UEXP 

0.0 

0.001 

2.3 

UCRIT 

0.0 

1.0  x  104 

1.0  x  106 

DELTA 

0.6749 

-1.0  x  102 

1.0  x  102 

VMAX 

3.24  x  105 

1.0  x  104 

9.9  x  107 

XJ 

0.66 

0.1 

1.5 

LAMBDA 

0.0 

1.0  x  10"6 

1.0  x  10"1 

KAPPA 

2.0 

0.0 

999.0 

NFS 

1.17  x  1012 

1.0  x  1010 

4.0  x  1014 

NEFF 

0.0 

1.0  x  10"2 

1.0  x  101 

NSS 

1.0  x  1011 

0.0 

1.0  x  1014 

TPG 

-1.0 

-1.0 

1.0 

ETA 

0.1210 

-10.0 

10.0 

THETA 

0.2396 

-10.0 

10.0 
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Table  4.       Methods  Used  to  Obtain  Initial  Estimates  for  Input  Parameters  of  SUXES 


PARAMETER        DEFAULT  VALUE  USED  OR  HOW  OBTAINED 


LEVEL 
TYPE 

LD 

TOX 

NSUB 

VTO 


KP 


GAMMA 


PHI 


Default  value  =  3 

iV-channel  MOSFET  default  value 
P-channel  MOSFET  default  value 


1.0 
-1.0 


N- channel  MOSFET  default  value  =  0.38  /mi 
P-channel  MOSFET  default  value  =  0.45  /xm 

Default  value  =  63.5  nm 


NSUB  =  ( 


GAMMA(eox)12 


2(e„)q)[ 


(1)  This  is  measured  at  Vds  =  0.2  V  and  Vbs  =  0.0  V  for  the  n-channel 
devices. 

(2)  This  is  measured  at  Vds  =  —0.2  V  and  Vbs  =  0.0  V  for  the  p-channel 
devices. 

3]  A  five-point  maximum  slope  technique  is  used. 

4)  When  the  maximum  slope  is  found,  a  straight  line  is  extrapolated  to 
the  x-axis  and  0.5  Vd8  is  subtracted  from  the  Vg8  value  to  obtain  the 
threshold  voltage. 
(5)  VTO  =  the  threshold  voltage 

(1)  This  is  measured  at  Va8  =  0.2  V  and  Vb8  =  0.0  V  for  the  n-channel 
devices. 

(2)  This  is  measured  at  Va8  =  —0.2  V  and  Vb8  =  0.0  V  for  the  p-channel 
devices. 

/o\  T/'p    maxslope(L) 

W  ^r  -  Vds(W) 

(4)  An  unreasonable  value  of  KP  causes  the  computer  program  to  abort. 

(1)  GAMMA  =  |slope|  of  the  threshold  voltage  vs. 

[\/|Vb8|+tphif- v/tphif] 
2}  Assume  tphif  =  0.6  V. 

3)  For  the  n-channel  MOSFETs,  measurements  are  taken  at  Vbs  =  —3.0, 
-6.0,  and  -9.0  V  for  Vds  -  0.2  V. 

(4)  For  the  j?-channel  MOSFETs,  measurements  are  taken  at  Vb8  =  3.0, 
6.0,  and  9.0  V  for  Vds  =  -0.2  V. 

(5)  A  five-point  maximum  slope  technique  is  used.  When  the  maximum 
slope  is  found,  a  straight  line  is  extrapolated  to  the  x-axis  and  0.5  Vd8 
is  subtracted  from  the  Vg8  value  to  obtain  the  threshold  voltage. 

(6)  A  recorded  default  value  for  a  measured  threshold  voltage  causes  the 
computer  program  to  abort. 


NSUB 


)  where  nj  =  1.45  x  1010  cm"3  at  300  K 


PHI 

The  simulated  ring  oscillator  frequencies  vary  by  less 
than  0.2  %  for  values  of  n-,  between  1.00  X  1010  cm-3 
and  1.45  x  1010  cm"3. 
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uo 

UEXP 

UCRIT 

DELTA 


UO 


_  KP(TOX) 


VMAX 


XJ 

LAMBDA 
KAPPA 


NFS 


NEFF 

NSS 
TPG 

ETA 


THETA 


This  is  not  a  MOS3  SPICE  input  parameter. 
Default  value  =  0.0 

This  is  not  a  MOS3  SPICE  input  parameter. 
Default  value  =  0.0 

The  value  of  this  parameter  is  preset  to  be  the  value  it  typically  is  from 

previous  SUXES  runs. 

JV-channel  MOSFET  default  value  =  2.685 

P-channel  MOSFET  default  value  =  0.6749 

The  value  of  this  parameter  is  preset  to  be  the  value  it  typically  is  from 
previous  SUXES  runs. 

TV-channel  MOSFET  default  value  =  7.9  x  105  m/s 
P-channel  MOSFET  default  value  =  3.24  x  105  m/s 


XJ  = 


LP 
0.8 


This  is  not  a  MOS3  SPICE  input  parameter. 
Default  value  =0.0 

The  value  of  this  parameter  is  preset  to  be  the  value  that  helps  the 
theory  match  the  experiment. 
JV-channel  MOSFET  default  value  =  0.80 
P-channel  MOSFET  default  value  =  0.20 

The  value  of  this  parameter  is  preset  to  be  the  value  it  typically  is  from 
previous  SUXES  runs. 

iV-channel  MOSFET  default  value  =  2.28  x  1012  cm-2 
P-channel  MOSFET  default  value  =  1.17  x  1012  cm-2 

This  is  not  a  MOS3  SPICE  input  parameter. 
Default  value  =  0.0 

Default  value  =  1011  cm-2 

JV-channel  MOSFET  default  value  =  1.0 
P-channel  MOSFET  default  value  =  -1.0 

The  value  of  this  parameter  is  preset  to  be  the  value  it  typically  is  from 
previous  SUXES  runs. 

JV-channel  MOSFET  default  value  =  0.1306 
P-channel  MOSFET  default  value  =  0.1210 

The  value  of  this  parameter  is  preset  to  be  the  value  it  typically  is  from 
previous  SUXES  runs. 

JV-channel  MOSFET  default  value  =  0.3612  V"1 
P-channel  MOSFET  default  value  =  0.2396  V-1 
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Table  5.       A  P-Channel  MOSFET  Options  File  of  SUXES 


* 


FILENAME: 


POPT.FIL 


* 


*  Nsign  =  Number  of  significant  digits 

*  Maxfn  =  Maximum  number  of  function  evaluations 

*  Delta  =  Norm  of  gradient  bound 

*  Epsil  =  Sum  of  squares  difference 

*  Linit  =  Marquardt  parameter  initial  value 

*  Lscal  —  Marquardt  parameter  scale 

*  Luppr  =  Marquardt  parameter  upper  value 

*  Fcdsw  =  Forward/central  difference  switch 


Criteria 

Nsign=4  Maxfn=200  Delta=1.0  x  10-9  Epsil=1.0  X  10~ 
Linit=0.1  Lscal=2.0  Luppr=1000  Fcdsw=0.1 


RangVd      min=— 5.5      max=0.00  incr=0.0 


RangVg      min=— 5.5      max=0.00  incr=0.0 


* 


* 


RangVb 


min=0.0 


max =0.0 


incr=1.0 


* 


* 


Weight 
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Table  6. 


A  P-Channel  MOSFET  Strategy  File  of  SUXES 


FILENAME: 


PSTRAT.FIL 


parameters 

first 

second 

third 

fourth 

fifth 

iteration 

iteration 

iteration 

iteration 

iteration 

LEVEL 

0 

0 

0 

0 

0 

TYPE 

0 

0 

0 

0 

0 

LD 

0 

0 

1 

0 

0 

TOX 

0 

1 

0 

0 

1 

NSUB 

0 

0 

1 

0 

0 

VTO 

0 

1 

0 

0 

1 

KP 

0 

0 

1 

0 

0 

GAMMA 

0 

1 

0 

0 

1 

PHI 

0 

0 

0 

1 

0 

UO 

0 

0 

0 

1 

0 

UEXP 

0 

0 

0 

0 

0 

UCRIT 

0 

0 

0 

0 

0 

DELTA 

1 

1 

0 

0 

1 

VMAX 

0 

0 

0 

1 

0 

XJ 

0 

0 

1 

0 

0 

LAMBDA 

0 

0 

0 

0 

0 

KAPPA 

0 

0 

0 

1 

0 

NFS 

0 

1 

0 

0 

1 

NEFF 

0 

0 

0 

0 

0 

NSS 

0 

0 

0 

0 

0 

TPG 

0 

0 

0 

0 

0 

ETA 

1 

0 

0 

0 

0 

THETA 

1 

0 

0 

0 

0 
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Table  7.       The  Methods  Used  to  Obtain  Several  Input  Parameters  for  SPICE 


PARAMETER  DEFAULT  VALUE  USED  OR  HOW  OBTAINED 


CJSW 

Default  value  =  0.6  x  10~9  F/m 

CGDO 

CGDO  =  e-$g- 

where  LD  is  the  lateral  diffusion 

CGSO 

CGSO  =  CGDO 

CGBO 

CGBO  =  ftgg 

where  GO  is  the  design  rule  of  the  gate  overlap 
=  4  fim 

RSH 

j\T-channel  MOSFET  default  value  =  11.63  il/a 
P-channel  MOSFET  default  value  =  61.57  il/a 

JS 

Default  value  =  5.0  X  10-5  A/m2 

MJ 

Default  value  =  0.5 

MJSW 

Default  value  =0.5 

PB 

Default  value  =  0.9  V 

XQC 

Default  value  =  0.45 
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Table  8.       A  SPICE  File  to  Obtain  the  I-V  Characteristics  for  a  P-Channel  MOSFET 


*  FILENAME  =  PMODEL.FIL 


.MODEL  MODP  PMOS  LEVEL=3  LD=0.450UM  TOX=0.0631UM  NSUB=0.789E+15 
+  KP=10.30U  GAMMA=0.186  PHI=0.766  UO=205.7  DELTA=3.025  VMAX=315790. 
+  XJ=0.874UM  KAPPA=25.8  NFS=0.300E-f  12  ETA=0.0798  THETA=0.1628 
+  CJSW=0.6E-9  CGDO=0.246E-09  CGSO=0.246E-09  CGBO=0.109E-08 
+  RSH=61.57  JS=5.0E-5  TPG=-1.0  MJ=0.5  MJSW=0.5  PB=0.9  XQC=0.45 
+  VTO=-0.715V 


*  FILENAME  =  IVCHARP.FIL 

*  NOTE:      Use  with  PMODEL.FIL  to  make  IVCHARP.CEL 
# 

*  P-CHANNEL  MOS  OUTPUT  CHARACTERISTICS 
* 

M2  1  2  0  0  MODP  L=2.33UM  W=6.33UM  AD=42.25P  AS=42.25P  PD=26U 

+  PS=26U  NRD=1.0  NRS=1.0 

* 

VD  4  0 
VG  2  0 
VIDS  4  1 
♦ 

.OPTIONS  NODE  NOPAGE 
.DC  VD  0  -5  -.5  VG  0  -5  -1 
.PLOT  DC  I(VIDS) 
.END 
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Table  9.       A  SPICE  File  for  a  19-Stage  CMOS  Ring  Oscillator 

*  

*  FILENAME  =  NMODEL.FIL 

* 

.MODEL  MODN  NMOS  LEVEL=3  LD=0.408UM  TOX=0.0558UM  NSUB=0.778E+15 

+  KP=29.87U  GAMMA=0.302  PHI=0.591  UO=574.5  DELTA=1.659  VMAX=794207. 

+  XJ=0.475UM  KAPPA=24.2  NFS=0.228E+13  ETA=0.0939  THETA=0.1767 

+  CJSW=0.6E-9  CGDO=0.252E-09  CGSO=0.252E-09  CGBO=0.124E-08 

+  RSH=11.63  JS=5.0E-5  TPG=1.0  MJ=0.5  MJSW=0.5  PB=0.9  XQC=0.45 

+  VTO=0.784V 

* 

*  FILENAME  =  PMODEL.FIL 
♦ 

.MODEL  MODP  PMOS  LEVEL=3  LD=0.450UM  TOX=0.0631UM  NSUB=0.789E+15 

+  KP=10.30U  GAMMA=0.186  PHI=0.766  UO=205.7  DELTA=3.025  VMAX=315790. 

+  XJ=0.874UM  KAPPA=25.8  NFS=0.300E+12  ETA=0.0798  THETA=0.1628 

+  CJSW=0.6E-9  CGDO=0.246E-09  CGSO=0.246E-09  CGBO=0.109E-08 

+  RSH=61.57  JS=5.0E-5  TPG=-1.0  MJ=0.5  MJSW=0.5  PB=0.9  XQC=0.45 

+  VTO=-0.715V 

* 

♦ 

*  FILENAME  =  INVERT. FIL 
* 

.SUBCKT  INVERT  12  3 

Ml  100  1  0  0  MODN  L=2.33UM  W=2.33UM  AD=12.5P  AS=12.5P  PD=15U 
+  PS=15U  NRD=2.0  NRS=2.0 

M2  100  1  3  3  MODP  L=2.33UM  W=6.33UM  AD=42.25P  AS=42.25P  PD=26U 
+  PS=26U  NRD=1.0  NRS=1.0 
Rl  100  2  1332 
CI  2  0  0.043P 
.ENDS  INVERT 


*  FILENAME  =  X19.FIL 

*  NOTE:      Use  with  NMODEL.FIL,  PMODEL.FIL  and  IN  VERT. FIL  to  make  X19.CEL 

♦ 

XI  1  2  3  INVERT 
X2  2  4  3  INVERT 
X3  4  5  3  INVERT 
X4  5  6  3  INVERT 
X5  6  7  3  INVERT 
X6  7  8  3  INVERT 
X7  8  9  3  INVERT 
X8  9  10  3  INVERT 
X9  10  11  3  INVERT 
REND  11  200  1110 
CEND  200  0  .044P 
X10  200  12  3  INVERT 

XII  12  13  3  INVERT 
X12  13  14  3  INVERT 
X13  14  15  3  INVERT 
X14  15  16  3  INVERT 
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X15  16  17  3  INVERT 
X16  17  18  3  INVERT 
X17  18  19  3  INVERT 
X18  19  20  3  INVERT 
X19  20  400  3  INVERT 
RCIR  400  1  334 
CCIR  1  0  .073P 
RPP  1  101  1.175 
CPP  101  0  .322P 
CPROBE  101  0  .IP 
* 

.IC  V(l)=5  V(2)=0  V(4)=5  V(5)=0  V(6)=5  V(7)=0  V(8)=5 

.IC  V(9)=0  V(l0)=5  V(11)=0  V(12)=5  V(13)=0  V(14)=5  V(15)=0 

.IC  V(16)=5  V(17)=0  V(18)=5  V(19)=0  V(20)=5 

.IC  V(3)=5 

♦ 

VCC  3  0  5V 
* 

.OPTIONS  VNTOL=10UV  LIMPTS=500  ITL4=50  ITL5=0 

.TRAN  1.0NS  200NS  UIC 

.PLOT  TRAN  V(l) 

.PLOT  TRAN  V(2) 

.PLOT  TRAN  V(4) 

.END 
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similar  to  those  of  the  p-channel  MOSFETs.  MOSS  was  the  model  used  in  SPICE  2G.5 
for  all  the  simulations. 

To  perform  the  simulations,  a  models  file  is  needed  to  characterize  the  devices.  It  is  as- 
sumed that  SUXES's  optimized  parameters  have  the  same  definition  as  the  input  model 
parameters  for  SPICE.  Hence,  these  parameters  are  directly  inserted  into  the  models  file 
of  SPICE.  SUXES  supplies  most  of  the  model  parameter  values.  Several  additional  pa- 
rameters (CGDO,  CGSO,  and  CGBO)  are  calculated  solely  from  the  optimized  SUXES 
values  according  to  the  equations  given  in  table  7.  In  the  computer  procedure  KEYS,  the 
parameters  LEVEL  and  TYPE  are  predetermined  (as  specified  in  table  4)  and  not  allowed 
to  be  changed  by  SUXES.  Several  parameters  are  used  in  SPICE  but  cannot  be  optimized 
with  the  current  version  of  SUXES.  They  are:  CJS W,  RSH,  JS,  MJ,  MJSW,  PB,  and 
XQC.  These  parameters  are  considered  noncritical  [8].  The  parameters  mentioned  above 
are  briefly  described  in  table  2,  and  the  default  values  used  in  this  analysis  are  listed  in 
table  7. 

In  addition  to  the  models  file,  SPICE  requires  that  the  user  supply  the  circuit  elements 
and  connections,  any  voltage  or  current  sources,  the  initial  conditions,  the  type  of  analysis, 
the  output  desired,  and  other  pertinent  SPICE  information. 

To  obtain  realistic  simulations,  especially  as  device  dimensions  are  reduced,  it  becomes 
necessary  to  include  parasitic  elements  in  the  simulations.  The  parasitic  elements  included 
in  the  example  given  in  this  report  are  the  source  and  drain  diffusion  resistances,  the 
polysilicon  resistances  and  capacitances,  the  metal  resistances  and  capacitances,  and  the 
probe  capacitance. 

The  two  most  common  errors  that  occur  when  simulating  a  circuit  are  the  inability  of  the 
values  of  the  node  potentials  in  the  circuit  to  converge  and  the  internal  timestep  being  too 
small  in  the  transient  analysis.  To  avoid  these  two  problems,  the  following  can  be  done: 

(1)  Initialize  each  node  in  the  circuit  to  an  appropriate  starting  potential. 

(2)  Use  the  options  on  the  .OPTIONS  card  as  shown  in  table  9. 

(3)  If  all  else  fails,  modify  the  models  file.  Perhaps  a  different  set  of  parameters  is 
needed  or  the  values  for  the  existing  parameters  are  unacceptable.  For  comparison, 
refer  to  the  models  file  included  in  table  9. 

(4)  For  additional  help,  consult  reference  [9]. 

2.4  STAT 2 

STAT2  is  used  to  determine  correlation  coefficients  [3,10,11]  and  to  generate  wafer  maps 
for  KEYS.  The  sample  correlation  coefficient,  r,  is  a  measure  of  the  similarity  of  the 
spatial  variation  of  two  sets  of  data.  Consider  the  paired  observations  (a?i,yi),  (a?2 ,2/2)) 
. . .,  (xn,yn)  on  two  quantities.  If  a  large  value  of  Xi  implies  a  large  value  of  then  the 
quantities  are  said  to  be  positively  correlated,  and  r  is  positive.  If  a  large  value  of  Xi 
implies  a  small  value  of  y»,  then  the  quantities  are  said  to  be  negatively  correlated,  and  r 
is  negative.  If  a  large  value  of  X{  implies  nothing  about  yi,  then  x  and  y  are  said  to  be 
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uncorrelated,  and  r  is  near  zero.  The  statistic  r  is  given  as  follows: 

r  =        YA=Axi-x){yi-y)  (1) 
^/E,■=l(^-*)2]E.■=l(3/i-y)2], 

where  x  and  y  are  the  arithmetic  means  of  the  X{  and  yi  values  of  the  sample  summed 
over  the  n  points.  Note  that  r  must  take  on  values  in  the  range  [—1,1].  Stronger  positive 
or  negative  correlation  is  indicated  by  r  approaching  +1  or  —1,  respectively. 

The  output  of  SUXES  produces  a  set  of  files  readable  by  STAT2.  There  is  one  file  for  each 
parameter  of  interest.  An  example  of  such  a  STAT2-readable  file  is  shown  in  table  10.  Site 
location  information  can  be  understood  with  the  help  of  figure  3.  The  first  16  records  of 
the  file  contain  the  location  information  for  the  data  values  which  follow.  The  first  of  the 

16  records  gives  the  location  of  the  first  sites  in  rows  1  through  4.  The  second  record  gives 
the  location  of  the  last  sites  in  rows  1  through  4.  In  the  example  of  table  10,  the  first  row 
of  the  wafer  has  test  sites  in  columns  4  through  6,  the  second  row  of  the  wafer  has  test 
sites  in  columns  2  through  8,  etc.  The  third  and  fourth  records  give  the  location  of  the 
first  and  last  site  in  row  5.  When  all  rows  have  been  characterized,  remaining  positions 
within  the  first  16  records  of  the  file  contain  zeroes.  The  data  values  are  given  in  records 

17  through  the  end  of  the  file.  They  are  ordered  by  rows,  such  that  values  from  the  first 
row  of  sites  are  first,  followed  by  values  from  the  second  row  of  sites,  and  so  on  for  as  many 
rows  of  sites  as  are  present. 

In  order  for  STAT2  to  calculate  correlation  coefficients,  it  is  necessary  to  create  a  database 
and  enter  the  processed  output  of  SUXES  into  the  database.  The  database  is  created  by 
the  CRDB  program  (a  program  which  is  adjunct  to  and  supplied  with  STAT2)  [3]  such 
that  data  from  all  test  sites  are  included  in  the  database.  A  sequence  of  commands  such 
as  that  shown  in  table  11  directs  STAT2  to  enter  the  processed  output  of  SUXES  into 
the  database.  Each  of  the  STAT2-readable  files  is  assigned  and  read  (STAT2  commands 
ASG  and  REA),  values  less  than  10-12  are  excluded  (XLT),  statistics  are  calculated  for 
the  remaining  sites  (PRS),  replacement  values  are  calculated  for  the  excluded  sites  (AXP), 
and  the  data  set  is  written  to  the  database  (WDB).  A  total  of  43  parameters  are  thus 
processed.  Because  the  first  record  in  the  database  serves  as  a  directory,  writing  of  actual 
data  begins  at  record  2  or  entry  number  2  and  ends  at  entry  number  44. 

Having  constructed  the  database,  correlation  coefficients  are  calculated  using  the  STAT2 
GET  and  SDB  commands.  As  a  default,  the  measured  frequency,  simulated  frequency,  and 
the  percent  difference  are  used  as  reference  data  sets  (database  entries  4,  5,  and  6)  with 
which  all  other  data  sets  are  correlated.  The  resulting  correlation  coefficients  are  sorted  in 
KEYS  by  magnitude.  Table  12  contains  an  example  of  the  output.  This  table  shows  at  a 
glance  the  parameters  which  are  most  influential  in  determining  the  measured  frequency, 
simulated  frequency,  and  percent  difference. 

It  is  often  useful  to  see  the  spatial  variation  of  one  or  more  parameters  over  the  wafer.  Using 
a  sequence  of  commands  such  as  that  shown  in  table  13,  STAT2  can  produce  gray-tone 
maps  similar  to  the  one  shown  in  figure  4. 
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Table  10. 


Example  of  a  File  Suitable  for  Input  to  STAT2* 


z 

1 
1 

£ 

0 

Q 

o 

Q 

y 

o 

o 

4 

u 

u 

u 

/» 
0 

u 

0 

0 

U 

u 

0 

0 

u 

U 

u 

u 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.000E+00 
0.000E+00 
0.000E+00 
0.000E+00 
0.990E+01 
0.940E+01 
0.930E+01 
0.960E+01 
0.100E+02 
0.103E+02 
0.000E+00 
0.105E+02 
0.930E+01 
0.970E+01 
0.950E+01 
0.930E+01 
0.940E+01 
0.980E+01 
0.112E+02 
0.000E+00 
0.900E+01 
0.930E+01 
0.990E+01 
0.104E+02 
0.109E+02 
0.113E+02 
0.000E+00 
0.990E+01 
0.112E+02 


Starting  column  locations  for  rows  1  through  4 
Ending  column  locations  for  rows  1  through  4 
Starting  column  location  for  row  5 
Ending  column  location  for  row  5 
Zero  values  indicate  the  array  of  test  sites 
does  not  go  beyond  row  5.  The  64 
numbers  in  the  first  16  records  of 
the  file  provide  for  site  locations 
in  an  array  having  up  to  32  rows 
and  32  columns. 


Data  value  (for  example,  MCKT)  for  site  at  row  1,  column  4 
Data  value  for  site  at  row  1,  column  5 
Data  value  for  site  at  row  1,  column  6 
Data  value  for  site  at  row  2,  column  2 
Data  value  for  site  at  row  2,  column  3 


!  Data  value  for  site  at  row  4,  column  8 
!  Data  value  for  site  at  row  5,  column  4 
!  Data  value  for  site  at  row  5,  column  5 
!  Data  value  for  site  at  row  5,  column  6 


*Text  beginning  with  an  exclamation  point  is  explanatory  and  not  part  of  the  actual  file. 
The  original  file  had  nine  rows  of  sites.  In  the  interest  of  brevity  the  table  shows  an 
example  having  five  rows  and  nine  columns  as  shown  in  figure  3. 
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Table  11.  Example  of  a  STAT2  Command  Sequence  Used  to  Build  a  Database  and  Cal- 
culate Correlation  Coefficients* 


ASG  W.DAT 
REA  1  0 
XLT  1.0E-12 
PRS 
AXP 

WDB,1,1,1,1,2 
ASG  L.DAT 
REA  1  0 
XLT  1.0E-12 
PRS 
AXP 

WDB,1,1,1,1,3 

ASG  MCKT.DAT 

REA  1  0 

XLT  1.0E-12 

PRS 

AXP 

WDB,1,1,1,1,4 
ASG  SCKT.DAT 
REA  1  0 
XLT  1.0E-12 
PRS 
AXP 

WDB,1,1,1,1,5 

ASG  PDIFF.DAT 

REA  1  0 

XLT  1.0E-12 

PRS 

AXP 

WDB,1,1,1,1,6 


Process  the  file  W.DAT  as  explained  in  the  text.  W.DAT  is 
of  the  type  described  in  table  10.  This  becomes 
entry  number  2  in  the  database. 


!  Process  the  file  L.DAT  in  a  similar  manner  to  produce 
!  database  entry  number  3. 


Process  the  file  MCKT.DAT.  This  contains  the  measured 
frequency  information  and  becomes  database  entry 
number  4. 


!  Process  the  file  SCKT.DAT.  This  contains  the  simulated 
!  frequency  information  and  becomes  database  entry  5. 


Process  the  file  PDIFF.DAT.  This  contains  the  percent 
differences  between  the  measured  and  simulated 
frequencies  and  becomes  database  entry  6. 


ASG  PRMS2.DAT 

REA  1  0 

XLT  1.0E-12 

PRS 

AXP 

WDB,1,1,1,1,44 
GET  4 

SDB  ALL  COEF  0.0 
GET  5 

SDB  ALL  COEF  0.0 
GET  6 

SDB  ALL  COEF  0.0 
END 


Continue  processing  data  files  in  a  similar  manner,  concluding 
with  the  file  PRMS2.DAT  which  becomes  database 
entry  number  44. 


Calculate  the  sample  correlation  coefficients  for  database 
entries  4  (the  measured  frequency),  5  (the 
simulated  frequency),  and  6  (the  percent 
difference),  against  all  other  entries  in 
the  database. 


'Text  beginning  with  an  exclamation  point  is  explanatory  and  not  part  of  the  actual  file. 
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Table  12.  A  Rank- Ordering  of  the  Magnitudes  of  the  Ten  Highest  Correlated  Parame- 
ters with  Respect  to  (a)  the  Percent  Difference  Between  the  Measured  and  Simulated  Ring 
Oscillator  Frequencies,  (b)  the  Measured  Frequency,  and  (c)  the  Simulated  Frequency.  An 
n-  or  p-channel  device  designation  appears  before  most  of  the  parameters  to  distinguish 
between  the  device  types. 


(a)  The 

/ 1   \  mi 

(b)  The 

/    \  mi 

(c)  The 

Percent 

Measured 

Simulated 

Difference 

Frequency 

Frequency 

Corr. 

Corr. 

Corr. 

Parameter* 

Coef. 

Parameter** 

Coef. 

Parameter 

Coef. 

1 

PDIFF 

1.000 

MCKT 

1.000 

SCKT 

1.000 

2 

NKP 

0.915 

SCKT 

0.745 

NKP 

0.979 

3 

NUO 

0.909 

NKP 

0.660 

NUO 

0.972 

4 

SCKT 

0.886 

NUO 

0.659 

PDIFF 

0.886 

5 

PKP 

0.809 

NKAPPA 

0.618 

PKP 

0.877 

6 

NLD 

0.730 

PKP 

0.616 

MCKT 

0.745 

7 

NCGBO 

0.706 

PTOX 

0.511 

NSUB 

0.713 

8 

NSUB 

0.701 

NSUB 

0.451 

NLD 

0.698 

9 

NCGDO 

0.698 

NGAMMA 

0.404 

NCGBO 

0.678 

10 

NGAMMA 

0.672 

NLD 

0.385 

NCGDO 

0.670 

*Remodeling  these  parameters  can  change  the  percent  difference  between  the  measured 
and  simulated  ring  oscillator  frequencies. 

** Wafer  maps  of  these  parameters  may  suggest  strategies  to  improve  process  control. 


Table  13.       Example  of  a  STAT2  Command  Sequence  Used  to  Make  a  Gray- Tone  Map 


ASG  MCKT. DAT 

REA  1  0 

PRS 

XLT  1.0E-12 

AXP 

PRS 

MP4,8,1 ,1,0,0,0 

MEASURED  FREQUENCY  (MHz)  -  LOT  YY027  WAFER  Y 
END 
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3.   AN  EXAMPLE  USING  CMOS  TECHNOLOGY 


In  this  part  of  this  report,  an  example  of  the  application  of  KEYS  is  presented.  In  par- 
ticular, a  twin-tub,  n-substrate,  CMOS  bulk  process  using  four-micrometer  design  rules 
was  used  in  this  example  to  manufacture  65  test  chips  on  a  4-in.  (10.16-cm)  wafer.  A 
19-stage  ring  oscillator  was  selected  as  the  demonstration  circuit.  MOSFETs,  both  n- 
and  j)-channel,  with  the  same  dimensions  as  the  ring  oscillator  MOSFETs,  were  used  to 
measure  the  I-V  data  points  for  SUXES. 

This  part  is  divided  into  two  sections.  Section  3.1  describes  the  test  structure  measurement 
methods  and  section  3.2  discusses  the  output  which  consists  of  plots,  tabular  listings,  and 
wafer  maps. 

3.1   Test  Structure  Measurement  Methods 

Three  types  of  test  structures  were  probed  on  each  test  chip:  a  ring  oscillator,  cross-bridge 
resistors  [12,13],  and  n-  and  p-channel  MOSFETs  with  the  same  dimensions  as  the  ring 
oscillator  MOSFETs.  A  manual  probe  station  and  oscilloscope  were  used  to  probe  the 
ring  oscillators,  and  the  measurements  for  the  other  test  structures  were  taken  with  an 
Accutest  Model  3000  parametric  test  system. 

The  ring  oscillator  was  selected  as  the  demonstration  circuit.  It  is  a  simple  circuit  and  can 
be  used  to  compare  frequencies  among  chips.  Since  the  frequency  is  sensitive  to  the  model 
parameters,  it  is  also  a  measure  of  how  "reasonable"  the  critical  parameters  are  to  the  user. 
Critical  parameters  are  denned  here  to  be  those  parameters  which  have  the  greatest  effect 
on  the  ring  oscillator  frequency  [8] .  If  the  user  is  not  satisfied  with  the  percent  difference 
between  the  simulated  and  measured  frequencies,  then  the  value  of  a  critical  parameter  is 
considered  "unreasonable"  to  the  user. 

The  19-stage  ring  oscillator  is  shown  in  figure  5.  The  designed  channel  lengths  of  the  n-  and 
p-channel  MOSFETs  and  the  designed  channel  widths  of  the  n-channel  MOSFETs  are  4 
fj.m.  The  designed  channel  widths  of  the  p-channel  MOSFETs  are  8  ^m.  A  cross- sectional 
view  of  an  inverter  is  shown  in  figure  6.  The  n- wells  are  not  contacted  on  the  top  side  of 
the  wafer,  creating  a  small  voltage  drop  (approximately  0.1  V).  This  is  due  to  the  doping 
differences  between  the  substrate,  set  at  5.0  V,  and  the  n-well.  A  crude  estimate  of  this 
voltage  drop  assumes  an  abrupt  change  in  the  doping  concentrations  between  the  n-type 
well  and  n-type  substrate.  This  estimate  is: 


voltage  drop  = 


/NWELL\ 
SUB  J 


(2) 


where  NWELL  is  the  n-well  doping  concentration,  and  SUB  is  the  n-substrate  doping 
concentration.  This  estimate  for  the  voltage  drop  decreases  the  simulated  ring  oscillator 
frequency  by  0.3%.  To  minimize  interference  with  the  actual  signal,  low  capacitance  probes 
(0.01  pF)  are  used  to  probe  the  output  of  the  ring  oscillator.  Care  was  taken  to  minimize 
effects  that  could  alter  the  measured  results.  A  5-V  voltage  swing  was  obtained.  Start-up 
circuitry  is  not  incorporated  into  the  design  of  this  ring  oscillator;  however,  a  comparison 
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Figure  5.    Block  diagram  of  the  19-stage  CMOS  ring  oscillator. 
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I 

of  the  measured  waveforms  with  the  simulated  waveforms  reveals  that  the  fundamental 
oscillation  frequency  is  obtained. 

Four-micrometer  cross-bridge  resistors  were  probed  automatically  to  determine  sheet  re- 
sistance and  linewidth  values  for  n+  diffusion,  poly,  and  metal.  A  cross-bridge  resistor  is 
shown  in  figure  7.  For  details  on  how  to  obtain  the  sheet  resistance  and  linewidth  values 
from  these  test  structures,  see  references  [12,13].  The  polysilicon  linewidth  values  obtained 
from  this  measurement  are  representative  of  the  magnitude  and  changes  in  the  MOSFET 
channel-lengths  used  in  SPICE. 

The  third  group  of  test  structures  are  the  n-  and  p-channel  MOSFETs  with  the  same  di- 
mensions as  the  ring  oscillator  MOSFETs.  These  devices,  of  which  the  p-channel  MOSFET 
is  shown  in  figure  8,  were  used  to  measure  the  I-V  data  points  for  SUXES  and  to  obtain 
the  initial  estimates  for  the  parameters  NSUB,  VTO,  KP,  GAMMA,  PHI,  and  UO.  These 
parameters  are  defined  in  table  2.  The  voltage  values  used  to  obtain  these  measurements 
are  given  in  table  14. 

3.2  The  Output 

KEYS  can  produce  output  in  the  form  of  plots,  tabular  listings,  and  wafer  maps.  Each 
form  is  discussed  in  the  subsections  that  follow. 

3.2.1   The  Output  Plots 

For  each  chip,  plots  are  made  for  the  n-  and  p-channel  MOSFET  I-V  characteristics  and  a 
simulated  ring  oscillator  node  potential  as  a  function  of  time.  These  plots  verify  that  the 
parameter  values  are  satisfactory. 

The  I-V  characteristics  include  three  separate  sets  of  curves:  the  measured  I-V  data, 
SUXES's  predicted  I-V  data,  and  SPICE's  simulated  I-V  data.  Any  deviations  will  be 
visually  apparent  to  the  circuit  designer.  A  plot  for  a  p-channel  MOSFET  is  shown  in 
figure  9.  The  open  circles  refer  to  the  measured  I-V  data  points.  SUXES  predicts  the  I-V 
data  points  using  its  optimized  parameters  in  SPICE,  and  the  symbol  for  these  data  is 
a  filled  circle.  Finally,  a  solid  line  refers  to  the  simulated  I-V  data  points  obtained  from 
SPICE  using  the  models  file  discussed  in  section  2.2. 

For  the  whole  wafer,  good  agreement  is  seen  among  the  three  curves  for  the  p-channel 
MOSFETs.  This  is  also  the  case  for  the  n-channel  MOSFETs.  For  the  wafer,  the  average 
percent  deviation,  for  all  the  acceptable  chips,  between  the  measured  I-V  data  points  and 
SUXES's  predicted  I-V  data  points  for  the  j>-channel  MOSFETs  is  3.2%.  This  value  is 
3.4%  for  the  n-channel  MOSFETs.  An  acceptable  chip  is  one  whose  optimized  SUXES 
values  are  not  pegged  at  the  lower  or  upper  bounds. 

In  KEYS,  SPICE  is  also  used  to  simulate  the  response  of  an  optional  dynamic  circuit. 
The  measured  waveforms  are  compared  with  the  SPICE  simulated  waveforms.  Figure  10 
displays  a  simulated  ring  oscillator  waveform.  The  frequency  and  period  are  also  recorded 
on  these  plots.  The  ring  oscillator  is  oscillating  in  its  fundamental  mode  [14]  and  the  node 
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Figure  7.     Four-/xm  cross-bridge  resistor  for  determining  the  polysilicon  linewidth. 
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Figure  8.  P-channel  MOSFET  of  the  same  dimensions  as  the  p-channel  ring  oscillator 
MOSFETs. 
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Table  14.       The  Voltage  Values  Used  to  Measure  the  N-  and  P-Channel  MOSFETS 


Measurements 

I-V  Data  Point  Needed  for  the 

Measurements  Initial  Estimates 


n- channel  p-  channel  n- channel  p-  channel 


Device 

dimensions 

Ln=4  /im 
Wn=4  /im 

Lp=4  /im 
Wp=8  /xm 

Ln=4  /xm 
Wn=4  /xm 

Lp=4  /xm 
Wp=8  /xm 

vb8 

0.0  V 

0.0  V 

0,-3,-6,-9  V 

0,3,6,9  V 

0.1,0.2,0.3, 
0.5,1,1.5, 
2  ?  2 . 5  ?  3  ^ 
3.5,4,4.5, 
and  5  V 

-0.1,-0.2, 
-0.3,-0.5, 
-1,-1.5, 
-2,-2.5, 
-3,-3.5, 
-4,-4.5, 
and  -5  V 

not 
applicable 

not 
applicable 

vd8 

0.1,0.2,0.3, 
0.5,1,1.5, 
2,2.5,3, 
3.5,4,4.5, 
and  5  V 

-0.1,-0.2, 
-0.3,-0.5, 
-1,-1.5, 
-2,-2.5, 
-3,-3.5, 
-4,-4.5, 
and  -5  V 

0.2  V 

-0.2  V 

Vchuck 

0.0  V 

0.0  V 
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voltages  undergo  a  5-V  voltage  swing.  For  the  wafer,  the  average  percent  difference,  for  all 
the  acceptable  chips,  between  the  simulated  and  measured  ring  oscillator  frequencies  was 
8.6%. 

The  titles  for  figures  9  and  10  indicate  the  wafer  and  the  chip  location  on  the  wafer.  For 
the  I-V  characteristics,  an  n-  or  p-channel  device  designation  is  included.  For  example, 
PWYR4C7  is  a  plot  (fig.  9)  of  the  I-V  characteristics  for  a  p-channel  MOSFET  on  the  chip 
located  at  row  4  column  7  on  wafer  Y,  and  WYR5C8  is  a  plot  (fig.  10)  of  a  ring  oscillator 
node  potential  on  the  chip  located  at  row  5  column  8  on  wafer  Y. 

3.2.2   The  Tabular  Listings 

KEYS  generates  two  types  of  tabular  listings:  a  comprehensive  summary  and  rank-order- 
ings  of  the  correlation  coefficients.  The  comprehensive  summary  is  shown  in  table  15.  This 
table  provides  the  measurement  results,  the  optimized  values  of  SUXES's  parameters,  the 
simulation  results,  the  calculated  parameters  for  each  chip,  the  wafer  averages,  and  the  lot 
averages  when  appropriate.  The  comprehensive  summary  in  combination  with  the  plots 
verifies  that  the  I-V  data,  SUXES,  and  SPICE  are  working  properly. 

Preceding  the  results  from  each  chip,  a  "+"  or  "— "  sign  appears.  A  chip's  results  are  either 
accepted  or  rejected.  If  the  values  for  the  n-  and  p-channel  MOSFETs  lie  between  specified 
lower  and  upper  bounds  in  SUXES's  initial  estimates  file,  it  is  considered  an  "acceptable" 
chip  as  defined  by  the  user.  Preceding  the  results  of  an  "acceptable"  chip,  a  "+"  sign 
appears.  If  a  parameter's  value  is  pegged  at  one  of  these  boundaries,  then  preceeding  the 
results  for  that  chip,  a  "— "  sign  appears.  Assuming  the  I-V  data  look  reasonable  to  the 
user,  a  rejected  chip  can  become  an  accepted  chip  by  modifying  the  strategy  file.  This  is 
described  in  Appendix  B.  Chips  with  unreasonable  linewidth  and  frequency  measurements 
as  determined  by  the  user  are  removed  from  the  database. 

The  wafer  averages  follow  the  results  from  all  the  chips.  The  number  of  acceptable  chips 
out  of  the  total  number  of  chips  is  recorded.  The  format  of  the  wafer  averages  is  similar 
to  the  format  for  a  single  chip.  The  values  are  the  averages  of  the  corresponding  values 
from  the  acceptable  chips.  Throughout  this  table,  the  n-channel  MOSFET  information  is 
listed  before  the  p-channel  MOSFET  information. 

The  lot  averages,  when  appropriate,  follow  the  results  from  the  last  wafer.  The  total 
number  of  wafers  is  recorded.  The  format  for  the  lot  averages  is  similar  to  that  of  the 
wafer  averages.  The  values  are  the  corresponding  average  values  from  the  wafer  averages. 

For  brevity,  table  15  provides  the  results  from  five  chips  as  opposed  to  all  the  chips.  Hence, 
in  the  complete  listing,  the  average  values  would  correspond  to  the  wafer  averages.  For 
the  entire  wafer,  the  average  percent  difference  between  the  simulated  and  measured  ring 
oscillator  frequencies  was  8.6%.  Different  strategy  files  can  change  this  average  percent 
difference. 

The  second  type  of  tabular  listing  is  a  rank-ordering  of  the  spatial  correlation  coefficients 
with  respect  to  chosen  parameters.    The  rank-orderings  for  this  example  are  given  in 
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Table  15.  A  Tabular  Summary  of  Each  Chip's  Results  Similar  to  the  Wafer  Summary. 
For  brevity,  this  listing  was  made  for  five  chips.  The  n-channel  MOSFET  information  is 
listed  before  the  p-channel  MOSFET  information. 
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table  12.  These  listings  are  with  respect  to  the  following  three  chosen  parameters:  the 
percent  difference  between  the  measured  and  simulated  ring  oscillator  frequencies  (PDIFF), 
the  measured  ring  oscillator  frequency  (MCKT),  and  the  simulated  ring  oscillator  frequency 
(SCKT).  In  these  tables,  an  n-  or  j?-channel  device  designation  appears  before  most  of  the 
parameters  to  distinguish  between  the  device  types.  The  exceptions  are  the  parameters 
which  do  not  require  a  designation.  They  are:  PDIFF,  MCKT,  and  SCKT. 

The  parameters  at  the  top  of  the  rank-ordering  in  table  12(a),  if  reoptimized,  can  change 
PDIFF.  To  reoptimize  the  parameters,  the  strategy  file  may  be  altered  according  to  the  pro- 
cedure given  in  Appendix  B.  Remodeling  these  same  parameters  can  also  change  PDIFF. 
This  could  have  the  effect  of  decreasing  the  rms  deviation  between  the  measured  and 
predicted  I-V  data  points,  and  it  could  increase  the  number  of  acceptable  chips. 

Table  12(b)  lists  the  correlation  coefficient  of  each  parameter  with  respect  to  MCKT.  Given 
a  minimum  number  of  acceptable  chips  (as  predetermined  by  the  user),  a  wafer  map  of 
MCKT  and  wafer  maps  of  the  highly  correlated  parameters  with  respect  to  MCKT  may 
suggest  strategies  to  improve  process  control.  An  example  of  a  wafer  map  is  given  in 
figure  4  for  MCKT.  This  is  the  third  form  of  output. 

3.2.3  The  Wafer  Maps 

Wafer  maps  may  be  generated  from  the  results  of  the  correlation  coefficient  rank-orderings. 
If  a  minimum  number  of  acceptable  chips  (as  predetermined  by  the  user)  are  obtained, 
the  wafer  maps  are  generated  by  STAT2.  In  this  example,  wafer  maps  are  made  of  the 
following  parameters:  (1)  MCKT;  (2)  SCKT;  (3)  PDIFF;  (4)  the  two  highest  correlated 
parameters  with  respect  to  PDIFF;  excluding  the  maps  already  generated;  and  (5)  the 
two  highest  correlated  parameters  with  respect  to  MCKT,  excluding  the  maps  already 
generated.  An  example  of  a  wafer  map  is  given  in  figure  4  for  the  parameter  MCKT.  This 
map  along  with  wafer  maps  of  the  highly  correlated  parameters  with  respect  to  MCKT 
may  suggest  strategies  to  improve  process  control. 

The  user  determines  the  minimum  number  of  acceptable  chips.  If  this  value  is  too  low, 
meaningless  wafer  maps  can  be  generated.  Altering  the  strategy  file  (using  the  procedure 
in  Appendix  B)  for  the  parameters  at  the  top  of  table  12(a)  may  increase  the  number  of 
acceptable  chips.  Other  techniques  may  be  just  as  viable.  Meaningful  wafer  maps  are  an 
invaluable  tool  in  evaluating  and  improving  the  process.  Distinctive  patterns  may  exist 
which  may  be  caused  by  limitations  in  the  processing  equipment  or  the  manner  in  which 
such  equipment  is  used  to  handle  wafers.  For  example,  the  doping  is  at  a  minimum  at  the 
center  of  the  wafer  when  it  is  placed  in  the  middle  of  a  fully  loaded,  slotted  quartz  diffusion 
boat.  In  this  case,  the  flow  of  gas  is  maintained  parallel  to  the  face  of  the  wafer.  When 
this  flow  is  not  parallel  to  the  face  of  the  wafer,  other  patterns  may  result.  One  example 
is  two  doping  minima  [15]. 

4.  CONCLUSIONS 

This  report  describes  a  computer  procedure  called  KEYS  that  can  be  used  to  characterize 
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a  chip,  wafer,  or  lot.  KEYS  links  SUXES,  SPICE,  and  STAT2  according  to  the  flowchart 
in  figure  2.  The  software  used  in  KEYS  is  given  in  table  1. 

The  benefits  of  this  procedure  are  as  follows:  (1)  the  technique  of  running  SUXES  and 
SPICE  becomes  standardized;  (2)  a  considerable  amount  of  time  is  saved;  (3)  manual  in- 
tervention and  error  are  reduced;  (4)  the  procedure  can  be  modified  for  different  dynamic 
circuits;  (5)  the  positional  information  and  the  comprehensive  output  are  valuable  to  pro- 
cess engineers,  modelers,  and  designers;  and  (6)  modifying  KEYS  for  updated  versions  of 
SUXES,  SPICE,  and  STAT2  is  relatively  easy  as  long  as  the  authors  of  these  programs  do 
not  modify  their  input-output  formats. 

An  example  using  CMOS  technology  is  also  given.  Some  of  the  limitations  of  this  procedure 
for  the  example  are  as  follows:  (1)  the  program  aborts  if  the  I-V  data  are  not  characteristic 
of  a  typical  MOSFET;  (2)  the  boundary  values  in  SUXES  are  arbitrary;  (3)  the  user 
determines  the  minimum  number  of  acceptable  chips,  and  meaningless  wafer  maps  are  an 
indication  this  value  is  too  low;  (4)  considerable  effort  is  required  to  modify  the  strategy 
file,  according  to  the  procedure  given  in  Appendix  B,  to  increase  the  number  of  acceptable 
chips;  and  (5)  the  user  helps  to  determine,  via  the  strategy  file,  how  closely  the  simulated 
results  match  the  measured  results. 

The  output  of  this  software  addresses  problems  with  simulating,  modeling,  and  processing. 
With  regard  to  the  simulations,  the  wafer  summary,  a  portion  of  which  is  shown  in  table  15, 
determines  which  chips  are  "acceptable."  The  I-V  characteristics  and  the  ring  oscillator 
plots  of  these  chips  verify  that  the  parameter  values  are  realistic.  For  the  unacceptable 
chips,  this  information  along  with  the  parameter  values  help  determine  if  the  I-V  data, 
SUXES,  or  SPICE  need  reexamination.  Also,  this  information  helps  determine  if  it  is 
worthwhile  modifying  the  strategy  file  to  make  these  chips  acceptable.  With  regard  to 
modeling  efforts,  the  parameters  at  the  top  of  table  12(a),  if  remodeled  properly,  would 
be  responsible  for  decreasing  the  percent  difference  between  the  measured  and  simulated 
ring  oscillator  frequencies.  These  same  parameters,  if  reoptimized,  may  be  responsible  for 
increasing  the  number  of  acceptable  chips.  With  regard  to  the  processing  efforts,  wafer 
maps  reveal  distinctive  patterns  of  nonuniformity  in  semiconductor  processes.  A  wafer  map 
of  the  measured  frequency  and  wafer  maps  of  the  parameters  at  the  top  of  table  12(b)  may 
assist  to  develop  processes  which  minimize  nonuniformity.  For  example,  ring  oscillator 
frequency  variations  may  be  minimized  by  providing  tighter  controls  on  the  polysilicon 
linewidth  variations,  which  may  be  due  to  variations  in  film  thickness  [15]. 


5.  ACKNOWLEDGMENTS 

Many  thanks  go  to  J.  A.  Horst  who  measured  the  threshold  voltages  and  I-V  data  points, 
developed  the  Accutest  computer  software  to  obtain  these  measurements,  and  transferred 
these  data  to  the  VAX.  We  appreciate  the  software  made  available  to  us  by  M.  Gaitan. 
For  comments  on  this  report,  we  thank  H.  S.  Bennett,  A.  R.  Hefner,  and  L.  W.  Linholm. 
The  editorial  assistance  of  J.  A.  Gonzalez  was  invaluable.  The  help  or  support  of  all  others 
involved,  including  J.  Albers  and  C.  L.  Wilson,  is  appreciated. 


40 


6.  REFERENCES 


[I]  Doganis,  K.,  and  Dutton,  R.  W.,  SUXES:  Stanford  University  Extractor  of  Model 
Parameters  (Users  Manual),  U.S.  Army  Research  Office  Contract  No.  DAAG-29-80-K- 
C-0013  and  DARPA  Contract  No.  MDA903-79-C-0257,  Integrated  Circuits  Laboratory, 
Stanford  University,  Stanford,  CA  (November  1982). 

[2]  Nagel,  L.  W.,  SPICE2:  A  Computer  Program  to  Simulate  Semiconductor  Circuits, 
Memorandum  No.  ERL-M520,  Electronics  Research  Laboratory,  University  of  Califor- 
nia, Berkeley,  CA  (May  1975). 

[3]  Mattis,  R.  L.,  S emiconductor  Measurement  Technology:  A  FORTRAN  Program  for 
Analysis  of  Data  from  Microelectronic  Test  Structures,  NBS  Special  Publication  400- 
75  (July  1983);  Mattis,  R.  L.,  and  Zucker,  R.,  Release  Notes  for  STAT2  Version  1.31: 
An  Addendum  to  NBS  Special  Publication  400-75,  NBSIR  83-2779  (November  1983); 
Mattis,  R.  L.,  Release  Notes  for  STAT2  Version  1.7:  An  Addendum  to  NBS  Special 
Publication  400-75,  NBSIR  86-3333  (March  1986). 

[4]  Nye,  B.,  Sangiovanni-Vincentelli,  A.,  Spoto,  J.,  and  Tits,  A.,  DELIGHT. SPICE:  An 
Optimization-Based  System  for  the  Design  of  Integrated  Circuits,  Proc.  1983  IEEE 
Custom  Integrated  Circuits  Conf.,  pp.  233-238,  May  1983. 

[5]  Scharfetter,  D.  L.,  Tremain  Jr.,  R.  E.,  Oki,  T.  J.,  Doganis,  K.,  and  Chen,  S.,  An  In- 
tegrated IC  Process  Characterization  Facility,  I  EDM  Tech.  Dig.,  246-249  (December 
1983). 

[6]  Vladimirescu,  A.,  and  Liu,  S.,  The  Simulation  of  MOS  Integrated  Circuits  Using 
SPICE2,  Memorandum  No.  ERL-M80/7,  Electronics  Research  Laboratory,  University 
of  California,  Berkeley,  CA  (February  1980). 

[7]  Alexander,  D.  R.,  Antinone,  R.  J.,  and  Brown,  G.  W.,  SPICE2  MOS  Modeling  Hand- 
book, BDM/A-77-071-TR,  The  BDM  Corporation,  Albuquerque,  NM  (May  1977). 

[8]  Cassard,  J.  M.,  A  Sensitivity  Analysis  of  SPICE  Parameters  Using  an  Eleven-Stage 
Ring  Oscillator,  IEEE  Trans.  Electron  Devices  ED  -  31  (2),  264-269  (1984)  and 
IEEE  J.  Solid  State  Circuits  SC  -  19  (1),  130-135  (February  1984). 

[9]  Unit  Software  &  Consulting  Inc.,  SPICE2  Error  Messages  and  Explanations,  (November 
1982). 

[10]  Snedecor,  G.  W.,  and  Cochran,  W.  G.,  Statistical  Methods,  (The  Iowa  State  University 
Press,  Ames,  IA,  1967),  pp.  172-198. 

[II]  Mendenhall,  W.,  Introduction  to  Probability  and  Statistics,  (Duxbury  Press,  Belmont, 
CA,  1971),  pp.  276-279. 

[12]  Buehler,  M.  G.,  Grant,  S.  D.,  and  Thurber,  W.  R.,  Bridge  and  van  der  Pauw  Sheet 
Resistors  for  Characterizing  the  Linewidth  of  Conducting  Layers,  J.  Electrochem. 


41 


II 
I 

Soc.  125  (4),  650-654  (April  1978). 

[13]  Buehler,  M.  G.,  and  Linholm,  L.  W.,  Role  of  Test  Chips  in  Coordinating  Logic  and  Cir- 
cuit Design  and  Layout  Aids  for  VLSI,  Solid  State  Technology  24  (9),  68-74  (Septem- 
ber 1981). 

[14]  Sasaki,  N.,  Higher  Harmonic  Generation  in  CMOS/SOS  Ring  Oscillators,  IEEE 
Trans.  Electron  Devices  ED  -  29  (2),  280-283  (February  1982). 

[15]  Perloff,  D.  S.,  Wahl,  F.  E.,  and  Reimer,  J.  D.,  Contour  Maps  Reveal  Non-Uniformity 
In  Semiconductor  Processing,  Solid  State  Technology  20  (2),  31-36,  42  (February 
1977). 

[16]  Doganis,  K.,  and  Scharfetter,  D.  L.,  General  Optimization  and  Extraction  of  IC  Device 
Model  Parameters,  IEEE  Trans.  Electron  Devices  ED  -  30  (9),  1219-1228  (Septem- 
ber 1983). 


42 


Appendix  A  -  The  Program  Listing  for  the  Computer  Procedure  KEYS 


This  appendix  contains  the  computer  procedure  KEYS.  It  is  ordered  by  file  type.  The  four 
different  file-types  specified  are  command  procedures  (with  a  .COM  or  .XXX  extension), 
FORTRAN  files  (with  a  .FOR  extension),  data  files  (with  a  .FIL,  .DAT,  or  .XXX  exten- 
sion), and  editor  files  (with  a  .WAY  extension).  They  are  listed  in  this  order  and  the  files 
within  each  file-type  are  listed  in  the  order  of  their  occurrence.  The  command  procedures 
are  given  first  and  they  are: 

(1)  KEYS.COM 

(2)  LAST3.COM 

(3)  INKEY.XXX 

(4)  SUXES.XXX 

(5)  SPICEIV.XXX 

(6)  RMSKEY.COM 

(7)  STKEY.COM 

(8)  CORKEY.XXX 

(9)  MAPKEY.XXX 

The  FORTRAN  files  included  in  this  Appendix  are  listed  second  and  they  are: 

(1)  RCKEY.FOR 

(2)  SUKEY.FOR 

(3)  CAPKEY.FOR 

(4)  IVKEY.FOR 

(5)  SUMKEY.FOR 

(6)  RANKEY.FOR 

The  data  files  are  listed  third  and  they  are: 

(1)  NLW.FIL 

(2)  RC.FIL 

(3)  LINE.FIL 

(4)  LOTNUMBER.FIL 
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(5)  NPARIN.FIL 

(6)  PPARIN.FIL 

(7)  NINPUT.XXX 

(8)  PINPUT.XXX 

(9)  NOPT.FIL 

(10)  POPT.FIL 

(11)  NSTRAT.FIL 

(12)  PSTRAT.FIL 

(13)  IVCHARN.XXX 

(14)  IVCHARP.XXX 

(15)  NMODEL.XXX 

(16)  PMODEL.XXX 

(17)  SUBCKT.XXX 

(18)  CIRCUIT.FIL 

(19)  PCHIPS.FIL 

(20)  STEND.FIL 

(21)  COREF.FIL 

The  editor  files  that  are  listed  fourth  are: 

(1)  EDNO.WAY 

(2)  EDSUXN.WAY 

(3)  EDSUXP.WAY 

(4)  EDPAR.WAY 

(5)  EDSPICEN.WAY 

(6)  EDSPICEP.WAY 

(7)  EDLAST.WAY 

(8)  EDNRMS.WAY 

(9)  EDPRMS.WAY 
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(10)  EDRES.WAY 

(11)  EDSKEY.WAY 

(12)  EDMAP.WAY 

(13)  EDSIMF.WAY 

(14)  EDPDIFF.WAY 

(15)  EDMEAF.WAY 

(16)  EDMAPS.WAY 

Refer  to  table  C-4  for  a  more  complete  listing  of  the  file  names. 
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KEYS . COM 


$  ! 

$  !     To  run  KEYS  type:     SUBMIT/noprinter/notify  KEYS 
$  ! 

$   !     FILENAME  =  KEYS.COM 
$  ! 
$  ! 

$   !  CURDIR  =  F$ENVIRONMENT( "DEFAULT") 
$  SET  DEFAULT  'KEYSDIR' 
$  ! 

$  copy  nlw.fil  in.dat 
$  ! 

$  edit/edt/command=edno . way/output=blank. f il  in.dat 

$  copy  blank. fil  switchl.dat 

$  copy  blank. fil  plw.dat 

$  copy  blank. fil  nguest.dat 

$  copy  blank. fil  pguest.dat 

$  copy  blank. fil  nparspice.dat 

$  copy  blank. fil  pparspice.dat 

$  copy  blank. fil  nlwl.dat 

$  copy  blank. fil  plwl.dat 

$  copy  blank. fil  mapl.dat 

$  ! 

$  !     Determine  which  chip  to  consider 
$  ! 

$  run  RCKEY 
$  ! 
$  ! 

$  !     If  the  last  chip  is  being  considered,  set  the  switch  in  RMSKEY.COM 
$  ! 

$  copy/concat  blank . fil , rmskey . com, switchl . dat  in.yyy 
$  ! 

$  edit/edt/command=edlast . way/output=rmskey . com  in.yyy 
$  ! 

$  !  Modify  and  run  the  command  procedure  INKEY  for  the  appropriate  chip 
$  ! 

$  copy/concat  blank. f il, inkey. xxx.switchl.dat, nlw.fil, plw.dat  in.yyy 
$  ! 

$  @last3  2  in.yyy  LLL  WWW 

$  edit/edt/command=edt3.com/output=out.yyy  in.yyy 
$  @last3  3  out.yyy  DDD  BBB  AAA 

$  edit/edt/command=edt3 . com/output=in.yyy  out.yyy 
$  @last3  2  in.yyy  EEE  III 

$  edit/edt/command=edt3 . com/output=out . yyy  in.yyy 
$  @last3  3  out.yyy  QQQ  QQQ  QQQ 

$  edit/edt/command=edt3 . com/output=in.yyy  out.yyy 
$  @last3  3  in.yyy  QQQ  PPP  NNN 

$  edit/edt/command=edt3 . com/output=inkey . com  in.yyy 
$  ! 

$  @ INKEY 
$  ! 
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LAST3 . COM 


Procedure  LAST3.COM  --  Reads  the  last  1,  2  or  3  records  of  a  file 
and  generates  an  EDT  command  file  which  makes  substitutions  based 
on  the  contents  of  those  last  three  records.     PI  contains  the 
number  1,  2,  or  3  and  indicates  how  many  substitute  commands  are 
to  be  produced.     P2  contains  the  name  of  the  file  to  be  edited. 
The  output  file  is  always  EDT3.COM.     P3  contains  the  string  which 
is  to  be  replaced  by  the  string  in  the  last  line.     P4  contains  the 
string  which  is  to  be  replaced  by  the  next  to  the  last  line.  P5 
contains  the  string  which  is  to  be  replaced  by  the  second  to  last 
line . 

On  any  of  the  last  three  lines  which  are  processed,   if  there  are  any 
leading  spaces,  these  spaces  are  removed  and  the  text  on  that  line 
is  left  justified  in  the  EDT3.COM  file. 


$ 
$ 
$ 
$ 
$ 
$ 
$ 
$ 
$ 
$ 
$ 
$ 
$ 
$ 
$ 
$ 

$  WS   :=  "WRITE  SYS$OUTPUT  " 

$  IF  (Pl.EQ.l)    .OR.    (P1.EQ.2)    .OR.    (P1.EQ.3)  THEN  GOTO  LO 
$  WS  "Argument  missing  or  not  in  range  1  to  3" 
$  EXIT 

$  ! 

$  LO: 

$  OPEN/READ  INFILE  ' P2 ' 

$  ! 

$  !  Branch  depending  on  how  many  records  are  to  be  examined. 

$  ! 

$  IF  Pl.EQ.l  THEN  GOTO  LI 

$  IF  P1.EQ.2  THEN  GOTO  L2 

$  IF  P1.EQ.3  THEN  GOTO  L3 

$  ! 

$  !  Case  1:  Generate  EDT3.COM  based  on  just  the  last  record  in  the  file 

$  ! 
$  LI: 

$  READ / END_0 F_FI LE=T00_FEW  INFILE  REC0RD1 
$  Lll: 

$  READ/END_0F_FILE=L12  INFILE  REC0RD1 
$  GOTO  Lll 
$  ! 

$  L12: 

$  CLOSE  INFILE 

$  OPEN/WRITE  OUT FILE  EDT3.COM 
$  ! 

$   !  Left  justify  the  last  record. 
$  ! 

$  LENGTH  :=  ' F$LENGTH(REC0RD1) ' 

$  IF  LENGTH  .LE.   1  THEN  GOTO  L14 

$  FIRST_SPACE  :=  ' F$L0CATE( "  " , REC0RD1) ' 

$  IF  FIRST_SPACE  .NE.   0  THEN  GOTO  L14 

$  REC0RD1   :=  ' F$ EXTRACT (1, LENGTH, REC0RD1 ) ' 

$  ! 
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$  L14: 

$  WRITE  OUT  FILE  "SUBS  /"  P3 '  /  "  RECORD1 '  /  1 :  END" 

$  WRITE  OUT FILE  "T  END - 1 " 

$  WRITE  OUT FILE  "D  . : END" 

$  WRITE  OUT FILE  "EXIT" 

$  CLOSE  OUT FILE 

$  GOTO  AFTER_WRITE 

$  ! 

$  !  Case  2:  Construct  EDT3.COM  using  the  last  two  records  in  the  file. 

$  ! 
$  L2: 
$  NUM=0 

$  READ / END_0 F_F I LE=TOO_FEW  INFILE  REC0RD2 
$  NUM=1 

$  READ / END_0 F_F I LE=TOO_FEW  INFILE  RECORD 1 
$  NUM=2 
$  L21: 

$  READ / END_0 F_F I LE=L2 2  INFILE  REC0RD2 
$  NUM=1 

$  READ / END_0 F_F I LE=L2 2  INFILE  REC0RD1 
$  NUM=2 
$  GOTO  L21 
$  ! 

$  L22: 

$  CLOSE  INFILE 

$  OPEN/WRITE  OUT FILE  EDT3.COM 
$  ! 

$   !  Left  justify  the  RECORD 1 . 
$  ! 

$  LENGTH  :=  ' F$LENGTH(REC0RD1) ' 

$  IF  LENGTH  .LE.   1  THEN  GOTO  L24 

$  FIRST_SPACE  :=  ' F$LOCATE( "  " .RECORDl) ' 

$  IF  FIRST_SPACE  .NE.   0  THEN  GOTO  L24 

$  REC0RD1   :=  ' F$EXTRACT(1, LENGTH, RECORDl ) ' 

$  ! 

$  L24: 
$  ! 

$  !  Left  justify  REC0RD2 . 
$  ! 

$  LENGTH  :=  ' F$LENGTH(REC0RD2) ' 

$  IF  LENGTH  .LE.   1  THEN  GOTO  L25 

$  FIRST_SPACE  :=  ' F$LOCATE( "  " , REC0RD2 ) ' 

$  IF  FIRST_SPACE  .NE.  0  THEN  GOTO  L25 

$  RE CORD 2   :=  ' F$ EXTRACT (1, LENGTH , REC0RD2 ) ' 

$  ! 

$  L25: 

$  IF  NUM.EQ.l  THEN  WRITE  OUTFILE  "SUBS  / "  P3 ' / " REC0RD2 ' /  1 : END" 
$  IF  NUM.EQ.l  THEN  WRITE  OUTFILE  "SUBS  /"P4'/"R£C0RD1'/  1 ;  END" 
$  IF  NUM.EQ.2  THEN  WRITE  OUTFILE  "SUBS  /"  P3  '/" RECORDl ' /  1 :  END" 
$  IF  NUM.EQ.2  THEN  WRITE  OUTFILE  "SUBS  /"  P4 '  /"  REC0RD2  '  /  1 :  END" 
$  WRITE  OUTFILE  "T  END- 2" 
$  WRITE  OUTFILE  "D  . : END" 
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$  WRITE  OUTFILE  "EXIT" 
$  CLOSE  OUTFILE 
$  GOTO  AFTER_WRITE 

$  ! 
$  L3: 

$  ! 

$  !  Case  3:  Construct  EDT3.COM  using  last  3  records  of  file, 
$  ! 

$  NUM=0 

$  READ/END_OF_FILE=TOO_FEW  INFILE  REC0RD3 
$  NUM=1 

$  READ/END_OF_FILE=TOO_FEW  INFILE  REC0RD2 
$  NUM=2 

$  READ /END_OF_FILE=TOO_FEW  INFILE  RECORDl 
$  NUM=3 
$  L31: 

$  READ/END_0F_FILE=L3 2  INFILE  REC0RD3 
$  NUM=1 

$  READ /END_0F_FILE=L3  2  INFILE  REC0RD2 
$  NUM=2 

$  READ/END_OF_FILE=L32  INFILE  RECORDl 
$  NUM=3 
$  GOTO  L31 
$  ! 

$  L32: 

$  CLOSE  INFILE 

$  OPEN/WRITE  OUTFILE  EDT3.COM 
$  ! 

$   !  Left  justify  RECORDl. 
$  ! 

$  LENGTH  :=  ' F$LENGTH (RECORDl ) ' 

$  IF  LENGTH  .LE.   1  THEN  GOTO  L33 

$  FIRST_SPACE  :=  ' F$LOCATE( "  ".RECORDl)' 

$  IF  FIRST_SPACE  .NE.   0  THEN  GOTO  L33 

$  RECORDl   :=  ' F$ EXTRACT (1, LENGTH, RECORDl ) ' 

$  ! 

$  L33: 
$  ! 

$   !  Left  justify  REC0RD2 . 
$  ! 

$  LENGTH  :=  ' F$LENGTH(REC0RD2) ' 

$  IF  LENGTH  .LE.   1  THEN  GOTO  L34 

$  FIRST_SPACE  :=  ' F$LOCATE( "  ",REC0RD2)' 

$  IF  FIRST_SPACE   .NE.   0  THEN  GOTO  L34 

$  REC0RD2   :=  ' F$EXTRACT(1 .LENGTH, REC0RD2) ' 

$  ! 

$  L34: 
$  ! 

$   !   Left  justify  RECORD 3 . 
$  ! 

$  LENGTH  :=  ' F$ LENGTH (RECORD 3) ' 
$  IF  LENGTH  .LE.  1  THEN  GOTO  L35 
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$  FIRST_SPACE  :=  ' F$LOCATE( "  ",REC0RD3)' 

$  IF  FIRST_SPACE  .NE.  0  THEN  GOTO  L35 

$  RECORD 3   :=  ' F$EXTRACT ( 1 , LENGTH , RECORD3 ) ' 

$  ! 

$  L35: 


$ 

IF 

NUM. 

EQ. 

1 

THEN 

WRITE 

OUTFILE 

"SUBS 

/' 'P3'/' ' RECORD3 ' / 

1: 

:  END 

$ 

IF 

NUM. 

EQ. 

1 

THEN 

WRITE 

OUT FILE 

"SUBS 

/' 'P4'/' ' RECORD1 ' / 

1: 

:  END 

$ 

IF 

NUM. 

EQ. 

1 

THEN 

WRITE 

OUTFILE 

"SUBS 

/' 'P5'/' ' RECORD2 ' / 

1: 

:  END 

$ 

IF 

NUM. 

EQ. 

2 

THEN 

WRITE 

OUTFILE 

"SUBS 

/' 'P3'/' ' RECORD 2 ' / 

1: 

:  END 

$ 

IF 

NUM. 

EQ. 

2 

THEN 

WRITE 

OUTFILE 

"SUBS 

/' 'P4'/' ' RECORD 3 ' / 

1: 

:  END 

$ 

IF 

NUM. 

EQ. 

2 

THEN 

WRITE 

OUTFILE 

"SUBS 

/' 'P5'/' ' RECORD 1 ' / 

1: 

;  END 

$ 

IF  NUM.EQ. 

3 

THEN 

WRITE 

OUTFILE 

"SUBS 

/"P3'/'  'RECORD1'/ 

1: 

:  END 

$ 

IF 

NUM. 

EQ. 

3 

THEN 

WRITE 

OUTFILE 

"SUBS 

/"P4'/'  ' RECORD2 ' / 

1: 

:  END 

$ 

IF 

NUM. 

EQ. 

3 

THEN 

WRITE 

OUTFILE 

"SUBS 

/' 'P5'/' 'RECORD3'/ 

1: 

:  END 

$  WRITE  OUTFILE  "T  END- 3" 

$  WRITE  OUTFILE  "D  . : END " 

$  WRITE  OUTFILE  "EXIT" 

$  CLOSE  OUTFILE 

$  GOTO  AFTER_WRITE 

$  ! 

$  AFTER_WRITE: 

$  EXIT  " 

$  ! 

$  TOO_FEW: 

$  WS  "Too  few  records  in  the  file  'P2'" 

$  EXIT 
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INKEY.XXX 


$  ! 

$   !     FILENAME:  INKEY.COM 
$  ! 

$  !  In  the  program  SUKEY,  the  width  is  asked  for  first  then 

$  !  the  length. 

$  ! 
$  ! 

$  SET  DEFAULT  'KEYSDIR' 

$  copy  WAAARBBBCDDDVB.DAT  vbin.dat 

$  ! 

$  ! 

$  !  Consider  the  n-channel  MOSFETs. 

$  ! 
NNN 

$  IF  L  .EQS.   0  THEN  GOTO  PCHAN 
$  copy  nWAAARBBBCDDD.dat  in.dat 
$  copy  nparin.fil  in3.dat 
$  ! 

$  !  Gather  the  initial  estimates  for  SUXES . 
$  ! 

$  run  SUKEY 

III 

EEE 

$  copy  guest.dat  nguest.dat 

$  copy  out.dat  nWAAARBBBCDDDF.dat 

$  copy  out2.dat  nWAAARBBBCDDDS.dat 

$  ! 

$  !     Insert  the  parameters  in  the  SUXES  input  file. 
$  ! 

$  copy/concat  blank. fil , ninput .xxx , nguest . dat  in.yyy 
$  ! 

$  @last3  3  in.yyy  VW  VW  VW 

$  edit/edt/command=edt3 . com/output=out .yyy  in.yyy 
$  @last3  3  out. yyy  VW  VW  VW 

$  edit/edt/command=edt3 . com/output=in. yyy  out. yyy 
$  @last3  2  in.yyy  VW  VW 

$  edit/edt/command=edt3 . com/output=out .yyy  in.yyy 
$  @last3  2  out. yyy  HHH  WV 

$  edit/edt/command=edt3 . com/output=in.yyy  out. yyy 
$  @last3  2  in.yyy  KKK  VW 

$  edit/edt/command=edt3 . com/output=out .yyy  in.yyy 
$  @last3  2  out. yyy  JJJ  VW 

$  edit/edt/command=edt3 . com/output=in.yyy  out. yyy 
$  @last3  2  in.yyy  MMM  WV 

$  edit/edt/command=edt3 . com/output=out . yyy  in.yyy 
$  @last3  2  out. yyy  000  VW 

$  edit/edt/command=edt3 . com/output=in.yyy  out. yyy 
$  @last3  2  in.yyy  RRR  WV 

$  edit/edt/command=edt3.com/output=out.yyy  in.yyy 
$  @last3  2  out. yyy  SSS  VW 
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$  edit/edt/command=edt3 . com/output=in.yyy  out.yyy 

$  @last3  2  in.yyy  TTT  WV 

$  edit/edt/command=edt3.com/output=out.yyy  in.yyy 

$  @last3  2  out.yyy  UUU  WV 

$  edit/edt/command=edt3 . com/output=ninput . dat  out.yyy 

$  ! 

$  !  Modify  and  run  the  command  procedure  for  SUXES  for  the  appropriate 

$  !  chip  and  device  type . 

$  ! 

$  copy/concat  blank. fil , suxes . xxx.plw. dat  in.yyy 

$  ! 

$  @last3  2  in.yyy  WV  WV 

$  edit/edt/command=edt3 . com/output=out . yyy  in.yyy 

$  @last3  3  out.yyy  ggg  aza  ccc 

$  edit/edt/command=edt3 .com/output=in.yyy  out.yyy 

$  edit/edt/command=edsuxn. way/output=suxes . com  in.yyy 

$  ! 

$  copy  ninput.dat  input.dat 

$  copy  nopt.fil  opt. dat 

$  copy  nstrat.fil  strat.dat 

$  ! 

$  @SUXES 

$  ! 

$  !  Gather  the  extracted  parameters. 

$  ! 

$  copy  expar.dat  in. dat 

$  ! 

$  edit/edt/command=edpar.way/output=exparn.dat  in. dat 

$  ! 

$  copy  exparn.dat  expar.dat 

$  run  CAPKEY 

$  ! 

$  copy  parspice.dat  nparspice.dat 

$  ! 

$  !   Insert  the  width  and  length  into  the  SPICE  file  which  obtains  the 

$  !  IV  characteristics. 

$  ! 

$  copy/concat  blank. fil , ivcharn . xxx , nlw . fil  in.yyy 

$  ! 

$  @last3  2  in.yyy  BZB  FFF 

$  edit/edt/command=edt3 . com/output=out .yyy  in.yyy 

$  @last3  3  out.yyy  WV  WV  WV 

$  edit/edt/command=edt3 . com/output=ivchar . f il  out.yyy 

$  ! 

$  copy  ivchar.fil  ivcharn. fil 

$  ! 

$  !  Put  together  the  n- channel  MOSFET  models  file  for  SPICE. 

$  ! 

$  copy/concat  blank. fil .nrnodel . xxx .nparspice . dat  in.yyy 

$  ! 

$  @last3  3  in.yyy  vw  aqa  bqb 

$  edit/edt/command=edt3 . com/output=out . yyy  in.yyy 
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$  @last3  3  out.yyy  cqc  dqd  eqe 

$  edit/edt/command=edt3 . com/output=in.yyy  out.yyy 
$  @last3  3  in.yyy  fqf  gqg  hqh 

$  edit/edt/command=edt3.com/output=out.yyy  in.yyy 
$  @last3  3  out.yyy  iqi  jqj  kqk 

$  edit/edt/command==edt3 .  com/output=in.  yyy  out.yyy 
$  @last3  3  in.yyy  lql  mqm  nqn 

$  edit/edt/comniand=edt3  .  com/output=out  .yyy  in.yyy 
$  @last3  3  out.yyy  oqo  pqp  rqr 

$  edit/edt/command=edt3 . com/output=in . yyy  out.yyy 
$  @last3  1  in.yyy  sqs 

$  edit/edt/command=edt3 . com/output=nmodel . f il  in.yyy 
$  ! 

$  COPY/CONCAT  BLANK . FIL , NMODEL . FIL , IVCHARN . FIL  IVCHARN. CEL. 
$  copy  ivcharn.cel  ivchar.cel 
$  ! 

$  !  Modify  the  command  procedure  which  obtains  the  IV  characteristics 
$  ! 

$  copy/concat  blank. fil , spiceiv.xxx.nlwl . dat  in.yyy 
$  edit/edt/command=edspicen.way/output=out .yyy  in.yyy 
$  @last3  2  out.yyy  gzg  fzf 

$  edit/edt/command=edt3 . com/output=in.yyy  out.yyy 
$  @last3  3  in.yyy  eze  dzd  czc 

$  edit/edt/command=edt3 . com/output=spiceiv. com  in.yyy 
$  ! 

$   !     Simulate  the  n-channel  MOSFET  IV  characteristics. 
$  ! 

$  (§spiceiv 

$  ! 

$  !  

$  ! 

$   !  Repeat  the  above  for  the  P-CHANNEL  MOSFETs . 
$  ! 

$  PCHAN: 
PPP 

$  IF  M  . EQS .   0  THEN  GOTO  CKT 
$  copy  pWAAARBBBCDDD.dat  in. dat 
$  copy  pparin.fil  in3.dat 
$  ! 

$   !  Gather  the  initial  estimates  for  SUXES . 
$  ! 

$  run  SUKEY 

WWW 

LLL 

$  copy  guest.dat  pguest.dat 
$  copy  out. dat  pWAAARBBBCDDDF.dat 
$  copy  out2.dat  pWAAARBBBCDDDS.dat 
$  ! 

$   !   Insert  the  parameters  in  the  SUXES  input  file. 
$  ! 

$  copy/concat  blank. fil , pinput . xxx , pguest . dat  in.yyy 
$  ! 
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$  @last3  3  in.yyy  VW  VW  VW 

$  edit/edt/command=edt3 . com/output=out . yyy  in.yyy 
$  @last3  3  out. yyy  VW  VW  VW 

$  edit/edt/command=edt3 . com/output=in.yyy  out. yyy 
$  @last3  2  in.yyy  VW  VW 

$  edit/edt/command=edt3 . com/output=out .yyy  in.yyy 
$  @last3  2  out. yyy  HHH  VW 

$  edit/edt/command=edt3 .com/output=in.yyy  out. yyy 
$  @last3  2  in.yyy  KKK  VW 

$  edit/edt/command=edt3 . com/output=out .yyy  in.yyy 
$  @last3  2  out. yyy  JJJ  VW 

$  edit/edt/command=edt3 . com/output=in.yyy  out. yyy 
$  @last3  2  in.yyy  MMM  VW 

$  edit/edt/command=edt3.com/output=out.yyy  in.yyy 
$  @last3  2  out. yyy  000  VW 

$  edit/edt/command=edt3 . com/output=in.yyy  out. yyy 
$  @last3  2  in.yyy  RRR  WV 

$  edit/edt/command=edt3.com/output=out.yyy  in.yyy 
$  @last3  2  out. yyy  SSS  VW 

$  edit/edt/command=edt3 . com/output=in.yyy  out. yyy 
$  @last3  2  in.yyy  TTT  VW 

$  edit/edt/command=edt3 . com/output=out .yyy  in.yyy 
$  @last3  2  out. yyy  UUU  WV 

$  edit/edt/command=edt3 . com/output=p input . dat  out. yyy 
$  ! 

$  !  Modify  and  run  the  command  procedure  for  SUXES  for  the  appropriate 
$  !  chip  and  device  type . 

$  ! 

$  copy/concat  blank. fil , suxes .xxx.plw. dat  in.yyy 
$  ! 

$  @last3  2  in.yyy  VW  VW 

$  edit/edt/command=edt3.com/output=out.yyy  in.yyy 
$  @last3  3  out. yyy  ggg  aza  ccc 

$  edit/edt/command=edt3 .com/output=in.yyy  out. yyy 

$  edit/edt/command=edsuxp . way/output=suxes . com  in.yyy 

$  ! 

$  copy  pinput.dat  input.dat 
$  copy  popt.fil  opt. dat 
$  copy  pstrat.fil  strat.dat 
$  ! 

$  @SUXES 
$  ! 

$  !  Gather  the  extracted  parameters. 
$  ! 

$  copy  expar.dat  in. dat 
$  ! 

$  edit/edt/command=edpar.way/output=exparp.dat  in. dat 
$  ! 

$  copy  exparp.dat  expar.dat 
$  run  CAPKEY 

$  copy  parspice.dat  pparspice.dat 
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$  !  Insert  the  width  and  length  into  the  SPICE  file  which  obtains  the 
$  !  IV  characteristics. 

$  ! 

$  copy/concat  blank. fil , ivcharp . xxx.plw. dat  in.yyy 
$  ! 

$  @last3  2  in.yyy  BZB  FFF 

$  edit/edt/command=edt3 . com/output=out .yyy  in.yyy 
$  @last3  3  out. yyy  VW  VW  VW 

$  edit/edt/command=edt3 . com/output=ivchar . f il  out. yyy 
$  ! 

$  copy  ivchar.fil  ivcharp. fil 
$  ! 

$  !  Put  together  the  p- channel  MOSFET  models  file  for  SPICE. 
$  ! 

$  copy/concat  blank. fil .prnodel .xxx.pparspice . dat  in.yyy 
$  ! 

$  @last3  3  in.yyy  wv  aqa  bqb 

$  edit/edt/command=edt3 . com/output=out .yyy  in.yyy 
$  @last3  3  out. yyy  cqc  dqd  eqe 

$  edit/edt/command=edt3 . com/output=in.yyy  out. yyy 
$  @last3  3  in.yyy  fqf  gqg  hqh 

$  edit/edt/command=edt3 . com/output=out . yyy  in.yyy 
$  @last3  3  out. yyy  iqi  jqj  kqk 

$  edit/edt/command=edt3 . com/output=in.yyy  out. yyy 
$  @last3  3  in.yyy  lql  mqm  nqn 

$  edit/edt/command=edt3 . com/output=out .yyy  in.yyy 
$  @last3  3  out. yyy  oqo  pqp  rqr 

$  edit/edt/command=edt3 . com/output=in . yyy  out. yyy 
$  @last3  1  in.yyy  sqs 

$  edi t/edt/command=edt3 . com/output=pmodel .fil  in . yyy 
$  ! 

$  COPY/CONCAT  BLANK. FIL, PMODEL. FIL, IVCHARP. FIL  IVCHARP. CEL 
$  copy  ivcharp. eel  ivchar.cel 
$  ! 

$  !  Modify  the  command  procedure  which  obtains  the  IV  characteristics. 
$  ! 

$  copy/concat  blank. fil , spiceiv. xxx ,plwl . dat  in.yyy 
$  edit/edt/command=edspicep . way/output=out .yyy  in.yyy 
$  @last3  2  out. yyy  gzg  fzf 

$  edit/edt/command=edt3 . com/output=in. yyy  out. yyy 
$  @last3  3  in.yyy  eze  dzd  czc 

$  edit/edt/command=edt3 . com/output=spiceiv. com  in.yyy 
$  ! 

$  !  Simulate  the  p-channel  MOSFET  IV  characteristics. 
$  ! 

$  (aspicEiv 

i 

i-    

| 

$   !     Simulate  the  optional  dynamic  circuit. 
$  ! 

$  CKT: 
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$  ! 

$  IF  (L  . EQS .  0)   .AND.    (M  . EQS .  0)  THEN  GOTO  END 
$  ! 
$  ! 

$  !     Insert  the  correct  widths  and  lengths  into  the  SPICE  file 
$  ! 

$  copy/concat  blank. fil , subckt .xxx.plw. dat ,nlw. fil  in.yyy 
$  ! 

$  @last3  2  in.yyy  hzh  izi 

$  edit/edt/command=edt3.com/output=out.yyy  in.yyy 
$  @last3  3  out.yyy  VW  VW  VW 

$  edit/edt/command=edt3 . com/output= in.yyy  out.yyy 
$  @last3  2  in.yyy  jzj  kzk 

$  edit/edt/command=edt3 . com/output=out .yyy  in.yyy 
$  @last3  3  out.yyy  VW  VW  VW 

$  edit/edt/command=edt3 . com/output=subckt . f il  out.yyy 
$  ! 
$  ! 

$  !  Simulate  the  dynamic  circuit. 
$  ! 

$  copy/concat  - 

blank. f il .nmodel . f il .pmodel . f il , subckt . f il , circuit . f il  circuit . eel 
$  ! 

$  RUN  SPICENBS 
CIRCUIT. CEL 
CIRCUIT.DAT 
CIRCUIT. PLT 
$  ! 

$  !  Plot  the  results  of  the  dynamic  circuit. 
$  ! 

$  ASSIGN  CIRCUIT. PLT  FOR011 

$  RUN  S 5 PLOT 

G 

WAAARBBBCDDD$ 
TIME  (sec)$ 
NODE  VOLTAGE  (V)$ 
$  RUN  MCVAX 
DEVICE  LPAO:  P300 
READ  IOP020.DAT 
PLOT 
EXIT 

$!   SET  TERM/BRO 
$  END: 
$  ! 

$  !  Sum  up  the  results  for  this  chip. 
$  ! 

$  submit/noprinter/notify  RMS KEY 
$  ! 
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SUXES  .XXX 


$  ! 

$   !  FILENAME:  SUXES.COM 
$  ! 

$  !  This  command  procedure  runs  SUXES  which  obtains  the  values  for  most 
$  !  of  the  SPICE  parameters. 
$  ! 

$  purge  *.* 

$  RUN  S SUXES 

6,5 

N 

Y 

1 

INPUT 
2 

77WcccRazaCgggF 

77WcccRazaCgggS 

OPT 

STRAT 

N 

F 

22 

N 

EXPAR 
Q 

$  ! 
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SPICEIV.XXX 


$  ! 

$  !  FILENAME:  SPICEIV.COM 
$  ! 

$  SET  DEFAULT  ' KEYSDIR' 
$  ! 

$  !  Run  SPICE  to  obtain  the  IV  characteristics. 
$  ! 

$  RUN  SPICENBS 
IVCHAR.CEL 
IVCHAR.DAT 
IVCHAR.PLT 
$  ! 

$  !  Massage  the  data  for  the  plotting  program  PLOTMY. 
$  ! 

$  copy  777WczcRdzdCezeF.DAT  in.dat 
$  ! 

$  copy  777WczcRdzdCezeS.DAT  in2.dat 
$  ! 

$  copy  77722fzfgzgP.DAT  in3.dat 
$  ! 

$  copy  ivchar.plt  in4.dat 
$  ! 

$  RUN  IVKEY 
$  ! 
$  ! 

$   !   Plot  the  data  where : 
$  ! 

$   !     +  =  IV  data 
$   !     *  =  SUXES  data 
$  !     -  =  SPICE  data 
$  ! 
$  ! 

$  purge  *.* 
$  ! 

$  run  plotmy 

n 

8 

out . dat 

n 

1 

out2 . dat 

n 

1 

out3 . dat 

n 

1 

out5 . dat 

n 

1 

out6 . dat 
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n 
1 

out 7 . dat 

n 

1 

out8 . dat 

n 

1 

out9 . dat 

n 

1 

1 

777WczcRdzdCeze$ 

VDS  (V)$ 

IDS  (A)$ 

2 

n 

n 

y 
+ 

+ 


$  RUN  MCVAX 

DEVICE  LPAO:  P300 

READ  IOP020.DAT 

PLOT 

EXIT 

$  ! 


RMSKEY.COM 


$  ! 

$  !   FILENAME:  RMSKEY.COM 

$  ! 

$  !     Find  the  average  error  for  the  SUXES  runs  and  gather  the  parameter 

$  !  values  for  the  file  RESULTS . DAT . 

$  ! 

$  SET  DEFAULT  'KEYSDIR' 

$  ! 

$  copy/concat  blank. fil , results . dat , switchl . dat , freq . dat , sfreq . dat  in2.dat 

$  copy  sys$login: keys . log  in. dat 

$  ! 
$  L=l 

$  IF  L  . EQS .  0  THEN  GOTO  PCHAN 
$  ! 

$  edit/edt/command=ednrms . way/output=nrms . dat  in. dat 
$  ! 

$  copy/concat  blank. fil , in2 . dat .nparspice . dat ,nrms . dat  in3.dat 
$  copy  in3 . dat  in2 . dat 
$  ! 

$  PCHAN: 
$  M=l 

$  IF  M  .EQS.  0  THEN  GOTO  CONTINUE 
$  ! 

$  edit/edt/command=edprms . way/output=prms . dat  in. dat 
$  ! 

$  copy/concat  blank. fil , in2 . dat ,pparspice . dat ,prms . dat  in3.dat 
$  copy  in3 . dat  in2 . dat 
$  ! 

$  CONTINUE: 
$  ! 

$  edit/edt/command=edres . way/output=in. dat  in2.dat 
$  ! 

$  !  Put  the  parametric  results  into  an  appropriate  format  for  each  chip, 
$  !  each  wafer,  and  each  lot. 

$  ! 

$  purge  *.* 
$  ! 

$  run  SUMKEY 
$  ! 

$  copy  out. dat  results.dat 
$  ! 

$  purge  *.* 

$  purge  sys$login: * . LOG 

$  ! 
$  N-l 

$  IF  N  .EQS.  3  THEN  GOTO  LAST 
$  IF  N  .EQS.  0  THEN  GOTO  NEXT 
$  ! 

$  !  The  last  chip  on  a  wafer  was  just  analyzed.  Place  the  switches  into  the 
$  !  command  procedure  CORKEY  which  determine  whether  or  not  the  wafer 
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$  !  maps  are  to  be  made . 

$  ! 

$  copy/concat  blank. fil, corkey.xxx.switchl.dat, switch2.dat  in.yyy 

$  ! 

$  edit/edt/command=edmap.way/output=corkey.com  in.yyy 

$  !  • 

$  print  results.dat 

$  ! 

$  submit/noprinter/notify  STKEY 

$  GOTO  END 

$  ! 

$  NEXT: 

$  ! 

$  !  Consider  the  next  chip. 

$  ! 

$  submit/noprinter/notify  KEYS 

$  GOTO  END 

$  ! 

$  LAST: 

$  print  results.dat 

$  ! 

$  END: 

$  ! 
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STKEY.COM 


$ 

$  !  FILENAME:  STKEY.COM 
$ 

$  SET  DEFAULT  'KEYSDIR' 
$  ! 

$  !  Prepare  the  parametric  data  files  for  STAT2 . 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  w.dat 
$  copy/concat  blank. fil , s tend. fil , in. dat  w.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  l.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  l.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  wp.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  wp.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in . dat  lp . dat 
$  copy/concat  blank. fil , stend. fil , in. dat  lp.dat 
$  ! 

$  edit/edt/command=edskey .way/output=in. dat  mckt.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  mckt.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  sckt.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  sckt.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  pdiff.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  pdiff.dat 
$  ! 
$  ! 

$  edit/edt/command=edskey .way/output=in. dat  ncgdo.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  ncgdo.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  ncgbo.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  ncgbo.dat 
$  ! 

$  edi t/edt/command=edskey . way/output=in . dat  nrms . dat 
$  copy/concat  blank. fil , stend. fil , in. dat  nrms. dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  nld.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  nld.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  ntox.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  ntox.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in . dat  nsub . dat 
$  copy/concat  blank. fil , stend. fil , in. dat  nsub. dat 
$  ! 

$  edit/edt/command=edskey .way/output=in. dat  nvto . dat 
$  copy/concat  blank. fil , stend. fil , in. dat  nvto. dat 
$  ! 
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$  edit/edt/command=edskey . way/output=in. dat  nkp.dat 
$  copy/concat  blank. fil , s tend. fil , in. dat  nkp.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  nganuna.dat 
$  copy/concat  blank,  fil ,  s  tend,  fil ,  in.  dat  nganuna.dat 
$  ! 

$  edi t/edt/command=edskey . way/output=in . dat  nphi . dat 
$  copy/concat  blank. fil , s tend. fil , in. dat  nphi. dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  nuO . dat 
$  copy/concat  blank . fil , s tend. fil , in. dat  nuO.dat 
$  ! 

$  edit/edt/command=edskey .way/output=in. dat  ndelta.dat 
$  copy/concat  blank. fil , s tend. fil , in. dat  ndelta.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  nvmax.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  nvmax.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  nxj.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  nxj.dat 
$  ! 

$  edit/edt/command=edskey .way/output=in. dat  nkappa.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  nkappa.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  nnf s . dat 
$  copy/concat  blank . fil , stend . fil , in. dat  nnf s. dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  neta.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  neta.dat 
$  ! 

$  edit/edt/command=edskey .way/output=in. dat  ntheta.dat 
$  copy/concat  blank . fil , stend. fil , in. dat  ntheta.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  nerr.dat 
$  copy/concat  blank . fil , stend . fil , in . dat  nerr.dat 
$  ! 

$  purge  *.* 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  pcgdo.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  pcgdo.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  pcgbo.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  pcgbo.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  prms . dat 
$  copy/concat  blank . fil , stend. fil , in. dat  prms . dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  pld.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  pld.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  ptox.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  ptox.dat 
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$  ! 

$  edit/edt/command=edskey . way/output=in. dat  psub.dat 
$  copy/concat  blank. fil , s tend. fil , in. dat  psub.dat 
$  ! 

$  edit/edt/conunand=edskey .  way/output=in.  dat  pvto  .  dat 
$  copy/concat  blank. fil , s tend. fil , in. dat  pvto. dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  pkp.dat 
$  copy/concat  blank. fil , s tend. fil , in. dat  pkp.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  pganuna.dat 
$  copy/concat  blank,  fil ,  s  tend,  fil ,  in.  dat  pganuna.dat 
$  ! 

$  edit/edt/conunand=edskey . way/output=in. dat  pphi . dat 
$  copy/concat  blank. fil , s tend. fil , in. dat  pphi. dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  puO.dat 
$  copy/concat  blank. fil , s tend. fil , in. dat  puO.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  pdelta.dat 
$  copy/concat  blank. fil , s tend. fil , in. dat  pdelta.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  pvmax.dat 
$  copy/concat  blank. fil , s tend. fil , in. dat  pvmax.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  pxj . dat 
$  copy/concat  blank. fil , stend. fil , in. dat  pxj. dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  pkappa.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  pkappa.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  pnf s . dat 
$  copy/concat  blank. fil , stend. fil , in. dat  pnf s . dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  peta.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  peta.dat 
$  ! 

$  edit/edt/command=edskey .way/output=in. dat  ptheta.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  ptheta.dat 
$  ! 

$  edit/edt/command=edskey . way/output=in. dat  perr.dat 
$  copy/concat  blank. fil , stend. fil , in. dat  perr.dat 
$  ! 

$  purge  *.* 
$  ! 

$   !     Create  the  database  for  STAT2 . 
$  ! 

$  run  crdb 
SUXES'S  VALUES 
4 

$  ! 

$   !     Run  STAT2  to  fill  the  database  and  obtain  the  correlation  coefficients. 
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$  ! 

$  run  STAT 2 
asg  w.dat 
rea  1  0 
xlt  1.0E-12 
prs 
axp 

wdb  ,  1 , 1 , 1 , 1 ,  2 
ASG  1 . DAT 
REA  1  0 
XLT  l.OE-12 
PRS 
AXP 

WDB ,1,1,1,1,3 
asg  wp . dat 
rea  1  0 
xlt  l.OE-12 
prs 
axp 

wdb ,1,1,1,1,4 
ASG  lp.DAT 
REA  1  0 
XLT  l.OE-12 
PRS 
AXP 

WDB ,1,1,1,1,5 
asg  mckt.dat 
rea  1  0 
xlt  l.OE-12 
prs 
axp 

wdb  ,  1 , 1 , 1 , 1 ,  6 
ASG  sckt . DAT 
REA  1  0 
XLT  l.OE-12 
PRS 
AXP 

WDB ,1,1,1,1,7 
asg  pdiff.dat 
rea  1  0 
xlt  l.OE-12 
prs 
axp 

wdb  ,  1 , 1 , 1 , 1 ,  8 
asg  ncgdo.dat 
rea  1  0 
xlt  l.OE-12 
prs 
axp 

wdb ,  1 , 1 , 1 , 1 ,  9 
asg  ncgbo.dat 
rea  1  0 


xlt  l.OE-12 

prs 

axp 

wdb ,  1 , 1 , 1 , 1 , 10 

asg  nrms . dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb ,1,1, 1,1, 11 

asg  nld.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb ,  1 , 1 , 1 , 1 , 12 

asg  ntox.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1, 13 

asg  nsub . dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb ,  1 , 1 , 1 , 1 , 14 

asg  nvto . dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1, 15 

asg  nkp . dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1, 16 

asg  ngamma.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1, 17 

asg  nphi.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb ,  1 , 1 , 1 , 1 , 18 


asg  nuO.dat 
rea  1  0 
xlt  l.OE-12 
prs 
axp 

wdb , 1 , 1 , 1 , 1, 19 

asg  ndelta.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb ,  1 , 1 , 1 , 1 , 20 

asg  nvmax.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb ,1,1, 1,1, 21 

asg  nxj . dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb ,  1 , 1 , 1 , 1 , 22 

asg  nkappa.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1 ,1,1, 23 

asg  nnf s . dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb ,  1 , 1 , 1 , 1 , 24 

asg  neta.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb ,  1 , 1 , 1 , 1 , 25 

asg  ntheta.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1 ,26 
asg  nerr.dat 
rea  1  0 
xlt  l.OE-12 
prs 


axp 

wdb,l,l,l,l,27 

asg  pcgdo.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb ,1,1,1,1, 28 

asg  pcgbo.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb , 1 , 1 , 1 , 1 , 29 

asg  prms . dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb , 1 , 1 , 1 , 1 , 30 

asg  pld.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1, 31 

asg  ptox.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1, 32 

asg  psub.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1, 33 

asg  pvto.dat 

rea  1  0 

xgt  -l.OE-12 

prs 

axp 

wdb , 1 , 1 , 1 , 1 , 34 

asg  pkp.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1, 35 
asg  pgamma.dat 
rea  1  0 


xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1 ,36 

asg  pphi . dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb ,1,1, 1,1, 37 

asg  puO . dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1 , 1 , 1 , 1 , 38 

asg  pdelta.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb ,1,1, 1,1, 39 

asg  pvmax.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb , 1 , 1 , 1 , 1 ,40 

asg  pxj . dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1,41 

asg  pkappa.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb  ,  1 , 1 , 1 , 1 , 42 

asg  pnf s . dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1 ,43 

asg  peta.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb  ,  1 , 1 , 1 , 1 , 44 


asg  ptheta.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb.l, 1,1, 1,45 

asg  perr.dat 

rea  1  0 

xlt  l.OE-12 

prs 

axp 

wdb, 1,1, 1,1,46 
GET  6 

SDB  ALL  COEF  0.0 
GET  7 

SDB  ALL  COEF  0.0 
GET  8 

SDB  ALL  COEF  0.0 

end 

$  ! 

$  !  Gather  and  rank- order  the  correlation  coefficient  information. 
$  ! 

$  submit/noprinter/notify  CORKEY 
$  ! 
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CORKEY.XXX 


$  ! 

$   !  FILENAME:  CORKEY.COM 
$  ! 

$  SET  DEFAULT  'KEYSDIR' 
$  ! 

$  copy  sys$login: stkey . log  in2.dat 
$  ! 

$  !  Gather  and  rank-order  the  correlation  coefficients  with  respect  to  the 
$  !  simulated  frequency. 

$  ! 

$  edit/edt/command=edsimf . way/output=in. dat  in2.dat 
$  ! 

$  run  RANKEY 
$  print  out. dat 
$  ! 

$  !  Gather  and  rank-order  the  correlation  coefficients  with  respect  to  the 
$  !  percent  difference  between  the  measured  and  simulated  frequencies. 

$  ! 

$  edit/edt/command=edpdif f . way/output=in. dat  in2.dat 
$  ! 

$  run  RANKEY 

$  copy  out2.dat  mapl.dat 
$  print  out . dat 
$  ! 

$  !  Gather  and  rank- order  the  correlation  coefficients  with  respect  to  the 
$   !  measured  frequency. 

$  ! 

$  edit/edt/command=edmeaf . way/output=in. dat  in2.dat 
$  ! 

$  run  RANKEY 
$  copy  out2 . dat  map2 . dat 
$ !   copy  out . dat  save . dat 
$  print  out. dat 

$  ! 
$  K=l 
$  N=2 

$  IF  K  . EQS .   0  THEN  GOTO  MID 
$  ! 

$  !  Modify  the  command  procedure  for  MAPKEY  to  include  the  appropriate  lot 
$  !  and  wafer. 

$  ! 

$  copy/concat  blank. fil .mapkey .xxx , lotnumber . fil ,nlw. fil  in.zzz 

$  copy/concat  blank . fil , in. zzz ,map2 . dat , mapl . dat , switchl . dat  in.yyy 

$  ! 

$  edit/edt/command=edmaps . way/output=out .yyy  in.yyy 
$  @last3  3  out. yyy  EEE  DDD  CCC 

$  edit/edt/command=edt3 . com/output=in.yyy  out. yyy 
$  @last3  3  in.yyy  BBB  QQQ  QQQ 

$  edit/edt/command=edt3 . com/output=out .yyy  in.yyy 
$  @last3  3  out. yyy  QQQ  QQQ  ZZZ 
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$  edit/edt/command=edt3 . com/output=in.yyy  out.yyy 
$  @last3  1  in.yyy  AAA 

$  edit/edt/command=edt3 . com/output=mapkey . com  in  yyy 
$  ! 

$  submit/nop r inter/notify  MAPKEY 
$  GOTO  END 
$  ! 

$  MID: 

$  IF  N  . EQS .   2  THEN  GOTO  FINISH 
$  ! 

$  !  Consider  the  first  chip  on  the  next  wafer. 
$  ! 

$  submit/noprinter/notify  KEYS 
$  GOTO  END 
$  ! 

$  FINISH: 

$  copy  results.dat  oresults.dat 
$  END: 
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MAPKEY.XXX 

$  ! 

$   !  FILENAME:  MAPKEY.COM 
$  ! 

$  !     Generate  the  FIRST  wafer  map  for  the  measured  frequency. 
$  ! 

$  SET  DEFAULT  'KEYSDIR' 
$  run  crdb 
SUXES'S  VALUES 
4 

$  ! 

$  purge  *.* 
$  ! 

$  run  STAT 2 
asg  mckt.dat 
rea  1  0 
prs 

xlt  1.0E-12 

axp 

prs 

MP4 ,8,1,1,0,0,0 

MEASURED  FREQUENCY  (MHz)   -  LOT  AAA  WAFER  ZZZ 
END 

$  PRINT/NOFEED  STAT 2 . LOG 
$  ! 

$  !  Generate  the  SECOND  wafer  map  for  the  simulated  frequency. 
$  ! 

$  run  STAT 2 
asg  sckt.dat 
rea  1  0 
prs 

xlt  1.0E-12 

axp 

prs 

MP4, 8,1, 1,0,0,0 

SIMULATED  FREQUENCY  (MHz)   -  LOT  AAA  WAFER  ZZZ 
END 

$  PRINT/NOFEED  STAT 2 . LOG 
$  ! 

$  !     Generate  the  THIRD  wafer  map  for  the  percent  difference. 
$  ! 

$  run  STAT 2 
asg  pdiff.dat 
rea  1  0 
prs 

xlt  1.0E-12 

axp 

prs 

MP4 ,8,1,1,0,0,0 

PERCENT  DIFFERENCE  (%)   -  LOT  AAA  WAFER  ZZZ 
END 
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$  PRINT/NOFEED  STAT 2 . LOG 
$  ! 

$  !     Generate  the  FOURTH  wafer  map  for  the  highest  correlated  parameter 
$   !  with  respect  to  the  measured  frequency. 

$  ! 

$  run  STAT 2 
asg  BBB.dat 
rea  1  0 
prs 

xlt  1.0E-12 

axp 

prs 

MP4 ,8,1,1,0,0,0 

BBB  -  LOT  AAA  WAFER  ZZZ 

END 

$  PRINT/NOFEED  STAT 2 . LOG 
$  ! 

$  !     Generate  the  FIFTH  wafer  map  for  the  second  highest  correlated  parameter 
$  !  with  respect  to  the  measured  frequency. 

$  ! 

$  run  STAT 2 
asg  CCC.dat 
rea  1  0 
prs 

xlt  1.0E-12 

axp 

prs 

MP4 ,8,1,1,0,0,0 

CCC  -  LOT  AAA  WAFER  ZZZ 

END 

$  PRINT/NOFEED  STAT 2 . LOG 
$  ! 

$  !     Generate  the  SIXTH  wafer  map  for  the  highest  correlated  parameter 
$  !  with  respect  to  the  percent  difference. 

$  ! 

$  run  STAT 2 
asg  DDD.dat 
rea  1  0 
prs 

xlt  1.0E-12 

axp 

prs 

MP4 ,8,1,1,0,0,0 

DDD  -  LOT  AAA  WAFER  ZZZ 

END 

$  PRINT/NOFEED  STAT2.LOG 
$  ! 

$  !  Generate  the  SEVENTH  wafer  map  for  the  second  highest  correlated  parameter 
$  !  with  respect  to  the  percent  difference. 

$  ! 

$  run  STAT 2 
asg  EEE.dat 
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rea  1  0 
prs 

xlt  l.OE-12 

axp 

prs 

MP4, 8,1,1,0,0,0 

EEE  -  LOT  AAA  WAFER  ZZZ 

END 

$  PRINT/NOFEED  STAT 2 . LOG 
$  N=2 

$  IF  N  .EQS.   2  THEN  GOTO  FINISH 
$  ! 

$  !  Consider  the  first  chip  on  the  next  wafer. 
$  ! 

$  submit/nopr inter/notify  KEYS 
$  GOTO  END 
$  ! 

$  FINISH: 

$  copy  results.dat  oresults.dat 
$  END: 
$  ! 

$   !  THE  END  ! 
$  ! 
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RCKEY . FOR 


FILENAME :     RCKEY . FOR 

This  program  chooses  the  next  chip  to  consider  for  SUXES 
and  SPICE. 


DIMENSION  IW(25),  IR(IOO) ,  IC(IOO) 

character*4  lotnuml 
character*4  lotnum2 


377 
378 
C 


478 
C 
C 
C 

11 

C 

21 


open(uni 
open(uni 
open(uni 
open(uni 
open(uni 
open(uni 
open(uni 
open(uni 
open(uni 
open(uni 
open(uni 
open(uni 
open(uni 
ifirst=l 
assign=0 


t=20, 
t=30, 
t=40, 
t=45, 
t=46, 
6=50, 
t=55, 
t=60, 
t=61, 
6-70, 
t=90, 
t=98, 
t=99, 

.0 


file='rc.fil' ,  status=' old' ) 
file='line.fil' ,  status=' old' ) 
file=' in.dat' ,  status=' old' ) 
file='nmodel . fil' ,  s tatus-' unknown ' ) 
file='pmodel. fil' ,  s tatus=' unknown ' ) 
f ile='plw. dat ' ,  s tatus=' unknown ' ) 
f ile='plwl . dat ' ,  s tatus=' unknown ' ) 
file='nlw. fil' ,  s tatus=' unknown ' ) 
f ile='nlwl . dat' ,  s tatus=' unknown ' ) 
file='freq.dat' ,  status=' unknown' ) 
f ile=' switchl . dat ' ,  s tatus=' unknown ' ) 
f ile=' lotnumber . f il ' ,  status=' old' ) 
file=' results . dat ' ,  s tatus=' unknown ' ) 


READ  THE  FILE 
M0=0 

DO  378  M0=l,7 
READ  (20,377)IDUMMY 
FORMAT  (Al) 
CONTINUE 

M0=0 

DO  478  M0=l,5 

READ  (30,377)IDUMMY 

CONTINUE 

READ  THE  PREVIOUS  ROW  AND  COLUMN 

read  (40,11)  IWLW 

format  (Al) 

write  (90,21)  iwlw 

format  (lx.Al) 

READ  (40,41)  IRLW 

READ  (40,41)  ICLW 

if  (irlw  .ne.  0)  go  to  102 
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assign=l . 0 

read  (98,734)  lotnuml,  lotnum2 
format  (2A4) 


734 

C  

c 

C  WRITE  THE  HEADING 

C 

WRITE  (99,1) 

1  format (lx, '#' ) 
WRITE  (99,2) 

2  format(lx,'*    FILENAME  =  RESULTS . DAT' ) 
WRITE  (99,4) 

WRITE  (99,3)  lotnuml,  lotnum2 
format(lx,'*    Lot  Number  =  ',lx,2A4) 
WRITE  (99,4) 
format(lx, '*') 
WRITE  (99,5) 

format(lx,'*    Wl      R      C        Wn  Ln  Wp  '  , 

+  '  Lp  ') 

WRITE  (99,667) 

format(lx,'*  MCKT  SCKT  PDIFF') 

WRITE  (99,4) 
WRITE  (99,6) 

6  format  (lx,'*  CGDO  CGSO  CGBO' , 
+'            RMS  ') 

WRITE  (99,7) 

7  format  (lx,'*  LD  TOX  SUB  ', 
+'            VTO  KP') 

WRITE  (99,8) 

8  format ( lx , ' *  GAMMA  PHI  UO  '  , 
+'            DELTA  VMAX') 

WRITE  (99,9) 

9  format(lx,'*  XJ  KAPPA  NFS  ', 
+  '            ETA            THETA  ERR') 

write(99,666) 

666        format (lx, **************************************************** , 


3 
4 
5 

667 


C--- 
C 
C 
C 

102 
800 

12 
C 

22 


FIND  THIS  ROW  AND  COLUMN  IN  RC.FIL 

1=0 
1=1+1 

READ  (20,12,END=231)  IW(I) ,   IR(I) , 

format  (Al , lx, 12 , lx, 12) 

write  (90,22)  IW(I),  IR(I),  IC(I) 

format  (lx,Al,2x,Il,2x,Il) 

if  (assign  .eq.  0.0)  go  to  222 

iwlw=iw(I) 

irlw=ir (I) 

iclw=ic (I) 

iwaf=iw(I) 


IC(I) 
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IROW=IR(I) 
ICOL=IC(I) 
go  to  880 

222        if  (IW(I)   .ne.  iwlw)  go  to  800 

IF  ((IR(I)    .EQ.   IRLW)   .AND.    (IC(I)    .EQ.   ICLW) )  GO  TO  700 
GO  TO  800 

C 

C  CHIP  NOT  FOUND  THEREFORE  START  FROM  SCRATCH 

C 

231        iwaf=iw(l)!  Therefore  first  chip 
IR0W=IR(1) 
IC0L=IC(1) 
ind=l 
ifirst=0 
GO  TO  880 

C 

C  THE  CURRENT  ROW  AND  COLUMN  HAS  BEEN  FOUND. . .FIND  THE  NEXT  CHIP 

C 

700        READ  (20,12,END=900)  IW'(I)  ,  IR(I)  ,  IC(I) 
C  write  (90,22)  IW(I),  IR(I) ,  IC(I) 

iwaf=IW(I) 

IROW=IR(I) 

IC0L=IC(I) 

C 

C  FIND  THE  LENGTH  AND  WIDTH 

C 

880        READ  (30,15,   END=900)  IIW,   IIR,   IIC,  Wn,  ALn,  Wp ,  ALp ,   FREQ,  IFLG 
15  FORMAT  (2x,Al,2x,I2,lX,I2,lX,F5.2,lX,F5.2,lX,F5.2,lx,F5.2, 

+1x,F5.2,1x,I3) 
C  write(90,25)  iiw,  iir , iic , wn, aln, wp , alp , f req , if lg 

25  FORMAT  ( 3x , Al , 2x , 12 , IX, 12 , IX, F5 . 2 , IX, F5 . 2 , IX, F5 . 2 , lx , F5 . 2 , 

+lX,F5.2,lx,I3) 

if(iwaf  .ne.  iiw)  go  to  880 

IF  ((IROW  .EQ.   IIR)    .AND.    (ICOL  .EQ.   IIC))  GO  TO  770 
GO  TO  880 

C 

C  Reasonable  W,  L,  and  FREQ  ? 

C 

770        IF((((WN. GE. 0.00) .AND. (ALN. GE. 0.00))    .AND.    ( (WP . GE . 0 . 00) . AND . 
+(ALP.GE.0.00)))    .AND.    (FREQ . GE . 0 . 00) )  GO  TO  777 

C 

C  GO  TO  THE  NEXT  CHIP 

C 

if  (ifirst  .eq.  1)  GO  TO  700 

ind  =  ind  +  1 

iwaf  =  iw(ind) 

irow  =  ir(ind) 

icol  =  ic(ind) 

go  to  880 

C 

C  WRITE  THE  VALUES  IN  NLW.FIL,  PLW.DAT,   FREQ. DAT,  NLW1.DAT,  and  PLW1.DAT 

C 
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write 

/  C  A 

(60 

1  A  \ 

Li)) 

TT  T  A  T7 

IWAF 

T  TT1  T  rpT? 

WRITE 

/  C  A 

TD  ATT 

I  ROW 

WRITE 

(60 

41) 

T  A  AT 

I  COL 

WRITE 

(60 

51) 

T  TXT 

WN 

T  Tit  T 

WRITE 

(60 

CI  \ 

51) 

A  T  XT 

ALN 

write 

/  C  A 

(50 

1  A  \ 

10) 

iwaf 

WRITE 

/  C  A 

(50 

41) 

T  n  ATT 

I  ROW 

WRITE 

/  C  A 

(50 

41) 

ICOL 

T  TO  T  T< 

WRITE 

(50 

51) 

T  TTI 

WP 

WRITE 

(50 

CI  \ 

Dl) 

A  T  T5 

ALP 

write 

(  ~IC\ 

if|\ 
J.U  ) 

iwaf 

WRITE 

(70 

41) 

IR0W 

WRITE 

(70 

41) 

ICOL 

WRITE 

(70 

51) 

WN 

WRITE 

(70 

51) 

ALN 

WRITE 

(70 

51) 

WP 

WRITE 

(70 

51) 

ALP 

WRITE 

(70 

51) 

FREQ 

c 

LN  =  IFIX(ALN*10) 

LP  =  IFIX(ALP*10) 

MWN  =  IFIX(WN*10) 

IF  (WN  .GE.   10.0)  MWN=IFIX(WN) 

MWP  =  IFIX(WP*10) 

IF  (WP  .GE.   10.0)  MWP=IFIX(WP) 

C 

write  (61,53)  IWAF 
WRITE  (61,42)  IROW 
WRITE  (61,42)  ICOL 
WRITE  (61,42)  MWN 
WRITE  (61,42)  LN 
write  (55,53)  iwaf 
WRITE  (55,42)  IROW 
WRITE  (55,42)  ICOL 
WRITE  (55,42)  MWP 
WRITE  (55,42)  LP 

41  FORMAT  (12) 

42  FORMAT  (lx,I2) 
51          FORMAT  (F5.2) 
10  FORMAT  (Al) 

53  FORMAT  (lx,Al) 

GO  TO  888 

C 

900        WRITE  (6,*)   '  NO  MORE  CHIPS  LEFT  OR  NO  MATCH' 
IROW=0 
ICOL=0 
ALN=0 . 0 
ALP=0 . 0 
WN=0.0 
WP=0 . 0 
FREQ=0 . 0 

write  (60,10)  iwaf 
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WRITE 

(60,41) 

I  ROW 

WRITE 

(60,41) 

ICOL 

WRITE 

(60,51) 

WN 

WRITE 

(60,51) 

ALN 

write 

(50,10) 

iwaf 

WRITE 

(50,41) 

IROW 

WRITE 

(50,41) 

ICOL 

WRITE 

(50, 51) 

\  ~^      9  / 

WP 

WRITE 

(50,51) 

ALP 

write 

(70,10) 

iwaf 

WRITE 

(70,41) 

IROW 

WRITE 

(70,41) 

ICOL 

WRITE 

(70,51) 

WN 

WRITE 

(70,51) 

ALN 

WRITE 

(70,51) 

WP 

WRITE 

(70,51) 

ALP 

WRITE 

(70,51) 

FREQ 

888  CONTINUE 

C 

C  IFLG=1        JUST  N- CHANNEL 

C  IFLG=2       JUST  P- CHANNEL 

C  IFLG=3  BOTH 

C 

IF  (IFLG  .EQ.   3)  GO  TO  503 
IF  (IFLG  .EQ.   2)  GO  TO  502 
WRITE  (90,511) 
WRITE  (90,512) 
write(46,39) 
39  format(lx, '*» ) 

GO  TO  1000 

502  WRITE  (90,510) 
WRITE  (90,513) 
write(45 , 39) 
GO  TO  1000 

503  WRITE  (90,511) 
WRITE  (90,513) 

1000  CONTINUE 

510  FORMAT ( ' $  L=0') 

511  FORMAT ( ' $  L=l') 

512  FORMAT ( ' $  M=0') 

513  FORMAT ( ' $  M=l ' ) 
C 

C  should  the  maps  be  generated  ? 

C 

if  (ifirst  .eq.  0)  go  to  995 
994        read  (20 , 12 , end=999)  iwaf 2 ,   irow2,  icol2 

997        read  ( 30 , 15 , end=999)  iiw,   iir,  iic,  wn,  aln,  wp,  alp,  freq,  iflg 
if  (iwaf 2  .ne.   iiw)  go  to  997 

if  ((irow2  .eq.   iir)   .and.   (icol2  .eq.   iic))  go  to  996 
go  to  997 

996        if  ((((wn  .gt.  0.0)   .and.   (aln  .gt.  0.0))   .and.   ( (wp  .gt.  0.0) 
+.and.   (alp  .gt.  0.0)))   .and.   (freq  .gt.  0.0))  go  to  993 
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go  to  994 

c 

ii  i,  lwaii    .iic,    lwai ^   go  to  jjz 

995 

write(90 , 61) !     Same  wafer... 

C 1 

bl 

rormac^  9  in=u  ) 

„  „     +_  _  QQQ 

go  to  770 

write  (,  yu ,  oi)  !     Dirrerent  water... 

81 

rormat(  9  N=l  ) 

go  to  998 

C 

999 

if  (assign  .eq.  1.0)  wnte(90,72) 

72 

format( '$  N=3 ' ) !     One  chip... 

if  (assign  .eq.  1.0)  go  to  998 

C 

write  (90,71)!     End  of  lot... 

71 

formatC  $  N=2  ' ) 

C 

998 

continue 

END 
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SUKEY . FOR 


C 

C  FILENAME:     SUKEY. FOR 

C 

C  This  program  takes  the  IV  data  from  the  ACCUTEST  and 

C  makes  it  suitable  for  the  SUXES  program. 

C 
C 

DIMENSION  VDS(IOOO),  VGS(IOOO),  VBS(IOOO),  CIDS(IOOO) 
DIMENSION  Xl(lOO),  Yl(100) ,  VBSZ(IOO) 
DIMENSION  VDSBE(IOOO),  VGSBE(IOOO) ,  CIDSBE(IOOO) 
real  x(250) ,y(250) , vds2 (500 , 50) ,beta,vt 
integer  n2(50) ,np(50) 

C 

open(unit=20 ,  f ile=' in. dat ' ,  status^' old' ) 
open(unit=31 ,  f ile=' vbin. dat ' ,  status=' old' ) 
open(unit=32 ,   f ile=' in3 . dat ' ,   status=' old' ) 
open(unit=40 ,  f ile=' out . dat ' ,  status= ' unknown' ) 
open(unit=50 ,   f ile=' out2 . dat ' ,   status=' unknown' ) 
open(unit=80 ,  file=' guest. dat' ,  s tatus=' unknown ' ) 
open(unit=90 ,  f ile=' vbout . dat ' ,  s tatus=' unknown ' ) 

C 

ES  =  103.5E-12 
EOX  =  34.5E-12 
Q  =  1.602E-19 
u0=0 . 0 
xkp=0 . 0 
vto=0 . 0 

C 

MO=0 

DO  178  MO=l,16 
READ  (20,177)IDUMMY 

177  FORMAT  (Al) 

178  CONTINUE 
C 

WRITE  (40,2) 
WRITE  (50,2) 

2  FORMAT  ('*' ,5x, 'Data  from  the  ACCUTEST') 
write  (40,3) 

write  (50,3) 

3  format  ( ' *' , 5x, ' VDS ' , lOx, ' VGS ' , 15x , ' VBS ' , llx, ' IDS ' ) 
C 

write  (6,4) 

4  format('  Width  =  ') 
read  (5,*)  W 

C 

write  (6,5) 

5  format ('  Length  =  ') 
read  (5,*)  AL 

C 

write  (40,6)  W,  AL 


in  F/m 
in  F/m 
in  C 


82 


write  (50,6)  W,  AL 
6  format  (10x, 'Width/Length=' ,F5.2, '/' ,F5.2) 

C 

C  READ  THE  DATA  AND  WRITE  IT  INTO  THE  OUTPUT  FILES 

C 

1=0 

230  1=1+1 
KEY=40 

C*600        READ  (20,*,END=231)  VDS(I),  VGS(I),  VBS(I),  CIDS(I) 
600        READ  (20,666,END=231)  IDUM,   IONE,  VDS(I) 

READ  (20,666,END=231)  IDUM,   ITWO,  VGS(I) 

READ  (20,666,END=231)  IDUM,   ITHREE,  VBS(I) 

READ  (20,666,END=231)  IDUM,   IFOUR,  CIDS(I) 
666        format  (A4,  16,  1PE13.6) 

VDS(I)  =  (AINT(100.*VDS(I)))/100 

VGS(I)  =  (AINT(100.*VGS(I)))/100 

VBS(I)  =  (AINT(100.*VBS(I)))/100 

IF  (VDS(I)    .LT.  0.00)  CIDS(I)=-CIDS(I) 

IF  (ABS(VDS(I))   .GT.  0.35)  KEY=50 

WRITE  (KEY, 401)  VDS(I),  VGS(I),  VBS(I),  CIDS(I) 
401        FORMAT  (5X,  F5 . 2 ,  8X,  F5 . 2 ,  10X,  F5 . 2 ,  5X,  1PE15.5) 

GO  TO  230 

231  CONTINUE 
NUM=I-1 

C 

C  ORDER  ALL  THE  DATA  FOR  THE  BETA  PROGRAM 

C 

450  1=0 
KC=0 

VDSLAST=0 
400        VDSY=100 . 0 
K=0 
J=0 

DO  88  K=1,NUM 

IF  ((ABS(VDS(K))    .LT.  ABS(VDSY))  .AND. 
+  (ABS(VDS(K))    .GT.  ABS (VDSLAST) ) )  VDSY=VDS(K) 
88  CONTINUE 

IF  (ABS(VDSY)  .EQ.  VDSLAST)  GO  TO  927  !  LAST  VDS  VALUE  ALREADY  USED 
1=1+1  !  CURVE  NUMBER 

C 

C  FOUND  ONE  VDS  VALUE  NOW  FIND  THE  VGS  VALUES 

C 

VGSLAST=0 . 0 
86  VGSY=100.0 
KB=0 
KA=0 

DO  89  KA=1,NUM 

IF  (ABS (VDS (KA))    .NE.  ABS(VDSY))  GO  TO  89 
IF  ((ABS(VGS(KA))    .LT.  ABS(VGSY))  .AND. 
+  (ABS (VGS (KA))    .GT.  ABS (VGSLAST) ) )  GO  TO  82 
GO  TO  89 
82  VGSY=VGS(KA) 
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89 
C 
C 
C 


C 

76 
103 


801 


C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

927 


91 
C 


KB=KA 
CONTINUE 

FOUND  A  VGS  VALUE 

IF  (KB  .EQ.   0)  GO  TO  76 

VGSLAST=VGSY 

VDSLAST=VDSY 

KC=KC+1 

J=J+1 

VDS2(J,I)=VDSY 
X(J)=VGSY 
Y(J)=CIDS(KB) 
N2(I)=J 

IF  (NUM  .EQ.  KC)  GO  TO  927 
GO  TO  86 


!  CURVE  COMPLETE  -  NO  VGS  VALUE 


!   POINTS  PER  CURVE  INCREMENT 


!  POINTS  PER  CURVE 
!  OUT  OF  POINTS 


CONTINUE 

WRITE  (90,103)  I,  N2(I),  VDS2(1,I) 
FORMAT  ('  Number  of  points  for  curve  ' , i2 
+  '  with  Vds  =  ' ,fl0.4) 

CALL  GAIN  (Y,X,N2 (I) ,5 , VDS2 (1 , I) , BETA, VT) 
XKP=B  ETA* AL /W 

IF  (ABS(VDS2(1,I))    .NE.  0.2)  GO  TO  801 

XKPO=XKP 

VTO=VT 

CONTINUE 

WRITE  (90,*)   '  KP=  ',XKP,'  VT=  ' ,VT 
GO  TO  400 

DETERMINE : 

SUB  (IN  CM- 3) 

PHI  (IN  V) 

GAMMA  (IN  V**.5) 

XJ  (IN  UM) 

CGDO  (IN  F/M) 

CGBO  (IN  F/M) 

UO  (IN  CM**2/VSEC) 
FOR  THE  INPUTS: 

TOX  (IN  A) 

LD  (IN  UM) 


READ  AND  MASSAGE  THE  DATA 
MO=0 

do  91  M0=l,14 

read  (31 , 177) idummy 

continue 


is  ' , i3 , 


!  WRONG  VDS 
!   IN  (UA/V**2) 


!  NEXT  CURVE 


phif2=.6 
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K0=0 

do  37  K0=l,4 

read  (31 , 566 , end=131)  IDUM,  IONE,  VBS(KO) 
566        format  (A4,  16,  1PE13.6) 

read  (31 , 566 , end=131)  IDUM,  IONE,  Y1(K0) 
C*  READ  (31,*,END=131)  VBS(KO),  Yl(KO) 

Xl(KO)=sqrt(abs(VBS(KO))+phif2)  -  sqrt(phif2) 

Yl(KO)  =  ABS(Yl(KO)) 

write  (90,*)  VBS(KO),  Xl(KO) ,  Yl(KO) 
37  continue 
C 

call  lnf it (3 , 2 ,X1 , Yl , slopep , xint , rmspsub) 

write  (90,*)   ' slopep ', slopep 

KO=0 

do  47  KO=5,8 

read  (31 , 566 , end=131)  IDUM,  ITWO,  VBS(KO) 
read  (31 , 566 , end=131)  IDUM,  ITWO,  Yl(KO) 
C*  READ  (31,*,END=131)  VBS(KO),  Yl(KO) 

Xl(K0)=sqrt(abs(VBS(K0))+phif2)  -  sqrt(phif2) 
Y1(K0)  =  ABS(Y1(K0)) 

write  (90,*)  VBS(KO),  Xl(KO),  Yl(KO) 
47  continue 

call  lnf it(3 , 6 ,X1 , Yl , slopen.xint , rmsnsub) 
write  (90,*)   ' slopen' , slopen 

C 

C  READ  THE  DATA 

C 

READ  (32,177)   I DUMMY 

C 

IKEY=32 
OKEY=80 

130        READ  (IKEY,*,END=131)  tox,  xld 
C 

if(vdslast  . ge .  0.0)  gamma=abs (slopen) 

if(vdslast   .It.   0.0)  gamma=abs (slopep) 

if(vdslast  . ge .  0.0)  rmssub=rmsnsub 

if(vdslast  .It.  0.0)  rmssub=rmspsub 
C  write  (90,*)   ' ikey= ' , ikey 

C  write  (90,*)    ' gamma= ' , gamma 

sub  =  (1.0e-6/(2.*es*q))*(gamma*eox/(tox*1.0e-10))**2       !   IN  (CM-3) 
C  write  (90,*)  'sub=',sub 

phi  =  2.0*.02585*alog(sub/1.45el0) 

xj  =  xld/. 8 

cgdo  =  eox*xld*1.0e-6/(tox*1.0e-10)  !   IN  (F/M) 

cgso  =  cgdo 

cgbo  =  eox*4.*1.0e-6/(2.*tox*1.0e-10) 

uO  =  (xkpo*tox*1.0e-8*1.0e-6)/(eox*1.0e-2)  !  IN  ( (CM**2)/V-SEC) 

C 

WRITE  (OKEY.375) 

375  FORMAT  (IX,  'LD') 
WRITE  (OKEY.376)  xld 

376  FORMAT  (IX,  F5.3) 
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c 


TJPTT17 

(OKEY.373) 

373 

FORMAT 

(IX,  'TOX') 

WRITE 

(OKEY.374) 

tox 

374 

FORMAT 

(IX,  F5.1) 

C 

TJP  TTT7 

(OKEY.377) 

377 

FORMAT 

(IX, 'SUB') 

WRITE 

(OKEY.378) 

sub 

378 

FORMAT 

(1X.1E9.3) 

C 

write 

(okey,903) 

903 

format 

( lx ,   ' vto ' ) 

write 

(okey,904) 

vto 

904 

format 

(lx,f6.3) 

C 

Ti7>"  1*  fo 
WL1  IC 

(okey,905) 

905 

format 

( lx ,   ' xkp ' ) 

write 

(okey,906) 

xkpo 

906 

format 

(lx,F6.3) 

C 

TJPTTT7 

(0KEY,577) 

577 

FORMAT 

(IX, 'GAMMA') 

WRITE 

(0KEY.578) 

gamma 

578 

FORMAT 

(1X.F5.3) 

C 

(0KEY.477) 

477 

FORMAT 

(IX, 'PHI' ) 

WRITE 

(0KEY.478) 

phi 

478 

FORMAT 

(1X.F5.3) 

C 

(okey,901) 

901 

format 

(lx,  'u0') 

write 

(okey,902) 

uO 

902 

format 

(lx,F6.1) 

C 

TJP  TTT7 

(OKEY.677) 

677 

FORMAT 

(IX, 'XJ') 

WRITE 

(OKEY.678) 

xj 

678 

FORMAT 

(1X.F5.3) 

C 

TJP  T  TP 
WKi  IE. 

(OKEY,777) 

777 

FORMAT 

(IX, 'CGDO' ) 

WRITE 

(OKEY.778) 

cgdo 

778 

FORMAT 

(1X.1E9.3) 

C 

WRITE 

(OKEY.877) 

877 

FORMAT 

(IX, 'CGSO') 

WRITE 

(OKEY.878) 

cgso 

878 

FORMAT 

(1X.1E9.3) 

C 

WRITE 

(OKEY.977) 

977  FORMAT (IX,  'CGBO') 
WRITE  (OKEY.978)  cgbo 

978  FORMAT  (lx,lE9.3) 
C 

WRITE  (OKEY.107) 

107  FORMAT ( IX ,  'RMSSUB') 
WRITE  (OKEY.108)  rmssub 

108  FORMAT  (lx,F5.3) 
C 

131  CONTINUE 
END 

C 
C 

SUBROUTINE  LNFIT  (N.M.X.Y, SLOPE, XINT, RMS) 

C 
C 

DIMENSION  X(l) ,Y(1) 

DIMENSION  YINT(lOl) ,YCALC(101) 

C 

1=0 

SUMX=0 
SUMY=0 
SUMXY=0 
SUMXX=0 

C 

DO  10  I=M,M+N-1 
C  write  (6,*)   'X=' ,X(I) , ' Y=' , Y(I) 

SUMX=SUMX+X(I) 

SUMY=SUMY+Y(I) 

SUMXY=SUMXY+X ( I ) *Y( I ) 

SUMXX=SUMXX+X ( I ) *X ( I ) 
10  CONTINUE 
C 

ANUMS=N*SUMXY - SUMX*SUMY 
C  write  (6,*)   ' anums= ' , anums 

ANUMI=SUMXX*SUMY- SUMX*SUMXY 
C  write  (6,*)   ' anumi= ' , anumi 

DEN=N*SUMXX- SUMX*SUMX 
C  write  (6,*)  'den=',den 

SLOPE=ANUMS/DEN 
C  write  (6,*)   ' slope=' , slope 

YINT(M)=ANUMI/DEN 

XINT= - ANUMI /ANUMS 

C 

C  DETERMINE  THE  RMS  DEVIATION  IN  THE  Y  DIRECTION 

C 

J=0 
SQU=0 

DO  20  J=M,M+N-1 
YCALC(J)=SLOPE*X(J)+YINT(M) 
SQU=SQU+(Y(J)-YCALC(J))**2 
20  CONTINUE 
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c 

RMS=SQRT(SQU/N) 

C 

RETURN 
END 

C 

C  -  

c 

subroutine  gain(Id,Vg, icount ,NDER,Vds ,beta,vt) 

C 

real  Id(250) , Vg(250) , Vds .beta, vt 
integer  icount ,nu,np,NDER 

C 

slopem=0 . 0 
slopen=0 . 0 
np=NDER/2 

C 

do  1000  I =np+l , icount -np 

x=0 . 0 

y=0.0 

xy=0 . 0 

xsq=0 . 0 

nu=0 

iflag=l  !  n- channel 

if  (vg(2) .lt.vg(l))  iflag=-l  !  Therefore  p-channel 

C 

do  391  M0=I-np,I+np 
x=x+Vg(M0) 
y=y+Id(MO)*1.0E+06 
xy=xy+Id(MO)*Vg(MO)*1.0E+06 
xsq=xsq+Vg(MO)*Vg(MO) 
nu=nu+l 
391  continue 

C 

de 1 1 a=nu*xs  q - x*x 

if  (delta  .eq.  0.0)  go  to  360 

slopen=(nu*xy-x*y) /delta 

if  (if lag  .eq.  1  .and.  slopen  .gt.  slopem)  go  to  999 
if  (if lag  .eq.   -1  .and.  slopen  .It.  slopem)  go  to  999 
go  to  1000 

999  slopem=slopen  !  slopem  =  max  slope 

yintc=(xsq*y-x*xy) /delta 
xintc=-yintc/slopem 
if  (Vds  .gt.   .90)  go  to  1000 
Vt=xintc-(Vds/2.0) 
1000  continue 

C 

beta=slopem/Vds 
return 

360  write (6,*)   '  ***  Delta  equals  0  ***' 
stop 
end 
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CAPKEY . FOR 


377 
100 
C 

C 
C 
C 

331 


FILENAME :     CAPKEY . FOR 

This  program  recalculates  the  capacitances  given  the 
optimized  parameters  from  SUXES . 


open(unit=20 ,  file=' guest . dat' ,  status=' old' ) 
open(unit=40 ,  f ile=' expar . dat ' ,  status=' old' ) 
open(unit=60 ,  f ile='parspice . dat ' ,  s tatus=' unknown ' ) 

EOX  =  34.5E-12  !  in  F/m 

READ  THE  PRESUXES  DATA 

KEY  =  20 
k=0 

do  100  k=l,25 

read  (key , 377 , end=331)  dummy 

format  (Al) 

continue 

READ  (KEY,*,END=331)  DRMSSUB 
READ  THE  AFTSUXES  DATA 
IKEY  =  40 

READ  (IKEY,*,END=231)  XLD , trashl , trash2 
READ  (IKEY,*,END=231)  TOX, trashl , trash2 
READ  (IKEY,*,END=231)  XNSUB , trashl , trash2 
READ  (IKEY,*,END=231)  VTO , trashl , trash2 
READ  (IKEY,*,END=231)  XKP , trashl , trash2 
XKP  =  XKP*1.0E+6 

READ  (IKEY,*,END=231)  GAMMA, trashl , trash2 
READ  (IKEY,*,END=231)  PHI , trashl , trash2 
READ  (IKEY,*,END=231)  UO , trashl , trash2 
READ  (IKEY,*,END=231)  DELTA, trashl , trash2 
READ  (IKEY,*,END=231)  VMAX, trashl , trash2 
READ  (IKEY,*,END=231)  XJ , trashl , trash2 
READ  (IKEY,*,END=231)  XKAPPA , trashl , trash2 
READ  (IKEY,*,END=231)  XNFS , trashl , trash2 
READ  (IKEY,*,END=231)  ETA, trashl , trash2 
READ  (IKEY,*,END=231)  THETA, trashl , trash2 

RECALCULATE  THE  CAPACITANCES 

cgdo  =  eox*xld*1.0e-6/(tox*1.0e-10) 
cgso  =  cgdo 

cgbo  =  eox*4.*1.0e-6/(2.*tox*1.0e-10) 
ATOX  =  T0X*le-4 
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c 
c 
c 

231 


60 
61 
62 
63 
64 
65 
66 
67 
68 
69 


WRITE  THE  PERTINENT  DATA  INTO  PARSPICE.DAT 


MKEY  =  60 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
CONTINUE 
END 


MKEY 

60) 

XLD 

MKEY 

69) 

ATOX 

MKEY 

62) 

XNSUB 

MKEY 

63) 

VTO 

MKEY 

64) 

XKP 

MKEY, 

60) 

GAMMA 

MKEY, 

60) 

PHI 

MKEY, 

61) 

UO 

MKEY, 

65) 

DELTA 

MKEY, 

66) 

VMAX 

MKEY, 

60) 

XJ 

MKEY, 

61) 

XKAPPA 

MKEY, 

62) 

XNFS 

MKEY, 

67) 

ETA 

MKEY, 

67) 

THETA 

MKEY, 

68) 

CGDO 

MKEY, 

68) 

CGSO 

MKEY, 

68) 

CGBO 

MKEY, 

60) 

DRMSSUB 

(lx,F5.3) 
(lx,F6.1) 
(lx,lE10.3) 
(lx,F6.3) 
(lx,F6.2) 
(lx,F8.3) 
(lx,F9.0) 
(lx,F8.4) 
(lx,lE9.3) 
(lx,F6.4) 


) 
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FILENAME :     IVKEY . FOR 


PURPOSE:     THIS  PROGRAM  MASSAGES  THE  DATA  FROM  SUXES  AND  SPICE  SO 
IT  IS  SUITABLE  FOR  THE  PLOTTING  PACKAGE 


DIMENSION 

DIMENSION 

DIMENSION 

open(unit= 

open(unit= 

open(unit= 

open(unit= 

open(unit= 

open(unit= 

open(unit= 

open(unit= 

open(unit= 

open(unit= 

open(unit= 

open(unit= 

open(unit= 


VDS(IOOO) , 
CIDS(IOOO) 
VD(IOOO) , 


■10 
=20 
=30 
40 
=11 
=21 
=31 
=41 
=42 
=43 
=44 
=45 
=46 


file= 
file= 
file= 
file= 
file= 
file- 
file= 
file- 
file= 
file- 
file= 
file= 
file= 


VGS(IOOO) 
,  CNVBS(10 
CID(IOOO) 
in.dat' ,  s 
in2 . dat ' 
in3 . dat ' 
in4 . dat ' 
out . dat ' 
out2 . dat 
out3 . dat 
out4 . dat 
out5 . dat 
out6 . dat 
out7 . dat 
out8 . dat 
out9 . dat 


VBS(IOOO) 


tatus='old' ) 
status=' old' ) 
status=' old' ) 
status=' old' ) 
s  tatus= ' unknown ' ) 
status= ' unknown ' 
s  ta tus= ' unknown ' 
s  ta tus= ' unknown ' 
s  tatus= ' unknown ' 
s  t atus= ' unknown ' 
s  ta tus= ' unknown ' 
s  ta tus= ' unknown ' 
status= ' unknown ' 


READ  THE  DATA 


JKEY=0 

JKEY=JKEY+10 
J=0 

DO  378  J-1,3 

READ  (JKEY,377)IDUMMY 

FORMAT  (Al) 

CONTINUE 

IF  (JKEY  .LT.   30)  GO  TO  600 
READ  (40,377)  I DUMMY 
READ  (40,*)  N 
N=N/6 


MASSAGE  THE  DATA  FILES  AND  THE  FILES  FROM  SUXES 


IKEY=0 
L=0 

IKEY=IKEY+10 

1=0 

1=1+1 

READ  (IKEY,*,END=231)  VDS(I),  VGS(I),  VBS(I),  CIDS(I) 

AN=0.0 

AN=AN+1 . 0 

IF  (AN  .GT.   5.5)  GO  TO  230 

IF  (ABS(VGS(I))    .LT.    ((AN  -  1.0)+. 99))  GO  TO  230 
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IF  (ABS(VGS(I))    .GT.   (AN  +  .01))  GO  TO  100 


C 

C  CORRECT  VGS  VALUE  FOUND 

C 

IF  (ABS(VDS(I))  .LT.  (.09))  GO  TO  230 
IF  (ABS(VDS(I))  .LT.  (.11))  GO  TO  777 
IF  (ABS(VDS(I))  .LT.  (.29))  GO  TO  230 
IF  (ABS(VDS(I))  .LT.  (.31))  GO  TO  777 
PF=0 . 0 
240        PF=PF+0 . 5 

IF  (PF  .GT.   5.1)  GO  TO  230 
IF  (ABS(VDS(I))    .LT.    (PF-.005))  GO  TO  230 
IF  (ABS(VDS(I))    . LT.    (PF+.005))  GO  TO  777 
GO  TO  240 

C 

C  CORRECT  VDS  VALUE  FOUND 

C 

777  L=L+1 
A=1.0 

VD(L)=VDS(I) 

IF  (VD(L)    .LT.   0.0)  A=-1.0 

CID(L)=A*CIDS(I) 

GO  TO  230 
231        WRITE(IKEY+1,4)  L 
4  FORMAT(5X,I5) 

K=0 

DO  24000  K=1,L 

WRITE  (IKEY+1,*)  VD(K),  CID(K) 
24000  CONTINUE 

IF  (IKEY  .LT.   30)  GO  TO  500 

C 

C  CONSIDER  THE  PLOT  FILES  FROM  SPICE 

C 

440  1=0 
J=0 

330  J=J+1 

IF  (J   .GE.   7)  GO  TO  331 
WRITE  (40+J,*)  N 
550  1=1+1 

READ  (40,*,END=331)  VDS(I) ,  CIDS(I) 
WRITE  (40+J,*)  VDS(I),  CIDS(I) 
IF  (I   .EQ.  J*N)  GO  TO  330 
GO  TO  550 

331  CONTINUE 
END 
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SUMKEY . FOR 


c 

C  FILENAME :     SUMKEY . FOR 

C 

c 
c 
c 
c 


PURPOSE:     THIS  PROGRAM  ADDS  TO  THE  'RESULTS .DAT' 
THE  AVERAGES. 


FILE  AND  CALCULATES 


dimension  chip(50),  wave(50),  wsum(50) 
dimension  alave(50),  alsum(50) 


common  /plomi/  MARK, 
common  map (10, 10) 


IWAF(IOO),   lotnuml,  lotnum2 


character*l 
character*l 
character*l 
character*l 
character*4 
character*4 


IWAF 
IW 

IWAFCUR 
MARK 
lotnuml 
lotnum2 


open(unit=10 
open(unit=20 
open(unit=24 
open(unit=25 
open(unit=30 
open(unit=49 
open(unit=50 
open(unit=51 
open(unit=52 
open(unit=53 
open(unit=54 
open(unit=55 
open(unit=56 
open(unit=57 
open(unit=58 
open(unit=59 
open(unit=60 
open(unit=61 
open(unit=62 
open(unit=63 
open(unit=64 
open(unit=65 
open(unit=66 
open(unit=67 
open(unit=68 
open(unit=69 
open(unit=70 
open(unit=71 


file= 
name= 
file= 
file= 
name= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 


mapl.DAT',   s tatus=' unknown ' ) 
in.dat',  status=' old' ) 
nlw.fil',   status=' old' ) 
line.fil',  status=' old' ) 
out . dat ' ,   s  tatus= ' unknown ' ) 
wp . da t ' ,   s  tatus= ' unknown ' ) 
lp . da t ' ,   s  tatus= ' unknown ' ) 
w.dat',   status=' unknown' ) 
1 . dat ' ,   s  tatus= ' unknown ' ) 
mckt.dat',  status= ' unknown' ) 
sckt.dat',   s tatus=' unknown ' ) 
pdiff.dat',  status= ' unknown ' ) 
nvto.dat',   s tatus=' unknown ' ) 
nkp . dat ' ,   s  tatus= ' unknown ' ) 
ntox . dat ' ,  s  tatus= ' unknown ' ) 
nld.dat',   s tatus=' unknown ' ) 
nsub . dat ' ,   s  t atus= ' unknown ' ) 
nrms . da t ' ,   s  t a tus= ' unknown ' ) 
nphi.dat',   s tatus=' unknown ' ) 
ngamma . dat ' ,   s  tatus= ' unknown ' ) 


nxj . dat ' , 
ncgdo . dat 
ncgso . dat 
ncgbo . dat 
nerr . dat' 
n the t a. dat 
nkappa . dat 
ndelta . dat 


s  tatus= ' unknown ' ) 
status=' unknown' ) 
status=' unknown' ) 
s  t atus= ' unknown ' ) 
s  ta tus= ' unknown ' ) 
,  s  t atus= ' unknown ' ) 
,   s  ta tus= ' unknown ' ) 
,   status=' unknown' ) 
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open(unit=72 
open(unit=73 
open(unit=74 
open(unit=75 
open(unit=76 
open(unit=77 
open(unit=78 
open(unit=79 
open(unit=80 
open(unit=81 
open(unit=82 
open(unit=83 
open(unit=84 
open(unit=85 
open(unit=86 
open(unit=87 
open(unit=88 
open(unit=89 
open(unit=90 
open(unit=91 
open(unit=92 
open(unit=93 
open(unit=94 
open(unit=95 
open(unit=96 
open(unit=97 
open(unit=98 


file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file= 
file- 
file= 
file= 
file= 
file= 
file= 
file= 
file- 
file^ 
file= 
file= 
file= 
file* 
file- 
file= 
file= 
file= 


ne  ta . da t ' ,  s  ta tus= ' unknown ' ) 
nuO . da t ' ,  s  t a tus= ' unknown ' ) 
nvmax . da t ' ,  s  t a tus= ' unknown ' ) 
nnf s . da t ' ,  s  t atus= ' unknown ' ) 
pvto . dat ' ,   s  t a tus= ' unknown ' ) 
pkp . da t ' ,  s  ta tus= ' unknown ' ) 
ptox.dat' ,  status=' unknown' ) 
pld.dat' ,  s tatus=' unknown ' ) 
psub . dat ' ,  status=' unknown' ) 
prms . dat ' ,  s tatus= ' unknown ' ) 
pphi . dat' ,  status=' unknown' ) 
pgamma . dat ' ,  status= ' unknown' ) 
px j . da t ' ,  s  t a tus= ' unknown ' ) 
pcgdo . dat ' ,  status= ' unknown' ) 
pcgso.dat',  s tatus=' unknown ' ) 
pcgbo.dat',  s tatus=' unknown ' ) 
perr.dat',  s tatus=' unknown ' ) 
p  the  t a . da t ' ,  s  t a tus= ' unknown ' ) 
pkappa . dat ' ,   status= ' unknown' ) 
pdelta.dat',  s tatus=' unknown ' ) 
peta.dat',  s tatus=' unknown ' ) 
puO.dat' ,  status=' unknown' ) 
pvmax . dat ' ,  s  tatus= ' unknown ' ) 
pnf s . da t ' ,  s  t a tus= ' unknown ' ) 
switch2.dat',  status=' unknown' ) 
pchips.fil',  status=' old' ) 
lotnumber . f il ' ,  status=' old' ) 


C 

342 


bust=0.0 

idump=0 

imark=l 

irow=0 

icol=0 

j  row=l 

jcol=4 

nodata=0 

zfil=0.0 

write  (10,342) 
format(lx, ' AAAAAA' ) 
write  (10,342) 


the  program's  doing  fine 
the  number  of  zero  locations 
the  last  MARK 
the  row  number 
the  column  number 


!  no  more  data  or  error  reading  data 
!  new  chip  AND  wafer  ? 


327 
C 

C 

C--. 
C 
C 
C 


read  (24,327)  IWAFCUR 
format  (Al) 

CALL  INTRO ( chip , wave , wsum, alave , alsum.perchip , itotchip) 


READ  THE  DATA 

lot=0  !  the  number  of  wafers 
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N=0  !   the  number  of  acceptable  chips 

1=0  !  the  total  number  of  chips 

500  1=1+1 

501  krow=irow 
kcol=icol 

read  (20 ,477 , ERR=700 , END=700)  MARK 

if  (MARK  .eq.    '#')  go  to  507  !  new  data 

if  (MARK  .eq.    '*')  go  to  588  !  average  titles 

READ  (20,66,ERR=700,END=700)IWAF(I) , IROW, ICOL, CHIP(l) ,CHIP(2) , 
+  CHIP(46) ,  CHIP(47) 

READ  (20,67,ERR=700,END=700)  CHIP(3) , CHIP(4) , CHIP(5) 
READ  (20,477)  I DUMMY 

C 

C  N-  or  P-channel  or  both  ? 

C 

if  (I  .le.  1)  go  to  503 
M=I 

if  (iwaf(I)   .ne.   iwaf(I-l))  go  to  937 
511        if  ((iwaf(I)   .ne.   iwaf(I-l))   .and.   (N  .eq.  0))  1=0 

if  (iwaf(I)   .ne.   iwaf(I-l))  CALL  WAFER ( imark, n, i, chip, 
+  wave , wsum, alave , alsum , nodata) 
if(i  .eq.  0)  i=i+l 
iwaf(I)  =  iwaf(M) 

503        read  (25 , 69 , end=502) iw, ir , ic ,xndum,yndum,xpdum,ypdum, zdum, idev 
if  (iwaf(I)    .ne.   iw)  go  to  503 

if  ((ir  .eq.   irow)   .and.   (ic  .eq.   icol))  go  to  505 
go  to  503 

502  write  (6,*)   '  The  correct  chip  was  not  located  in  LINE.FIL-1' 
go  to  909 

505  inc=0 

if (idev  .eq.   2)  inc=20  !  just  p-channel 

C 

506  read  (20 , * , end=700)  chip ( 15+inc ), chip (16+inc) , 
+  chip ( 17+inc ), chip (11+inc) 

read  ( 20 , * , end=700 )  chip ( 9+inc ) , chip ( 8+inc ) , chip ( 10+inc ) , 
+  chip(6+inc) , chip(7+inc) 

read  (20 , * , end=700)  chip ( 13+inc ), chip (12+inc) , chip (23+inc) , 
+  chip ( 21+inc ), chip (24+inc) 

READ  (20,*,END=700)  chip(14+inc) , chip(20+inc) , chip(25+inc) , 
+  chip(22+inc) , chip(19+inc) , chip(18+inc) 

if  ((idev  .eq.  1)   .or.   ( inc  . ge .  20))  go  to  920       !  just  n-channel 

read  (20,477)  I DUMMY 

inc  =  20 

go  to  506 

C 

507  if  (mark  .eq.   '#')  CALL  NEWCHIP(n, i , iw, irow, icol , chip , 
+  wave ,wsum, alave , alsum, idev, bus t , nodata, imark, zf il) 

if  (zfil  .ne.  0.0)  go  to  937 
if  (BUST  .eq.   1.0)  go  to  937 

C  -  

c 

C  WRITE  THE  DATA  INTO  OUT. DAT 
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c 

920  zfil=0.0 

if  (mark  .eq.   '-')  imark=-l 
if  (mark  .eq.   '+')  imark=l 
write  (30,477)  MARK 

WRITE  (30,41)  IWAF(I) , IROW, ICOL.chip(l) , chip(2) , chip(46) ,chip(47) 
WRITE  (30,67)  chip(3) ,chip(4) ,chip(5) 
WRITE  (30,*)   '  ' 
inc=0 

if  (idev  .eq.  2)  inc=20 
538        write  (30,42)  chip ( 15+inc ), chip (16+inc) , chip (17+inc) , chip (11+inc) 

write  (30,43)  chip ( 9+inc ), chip (8+inc) , chip (10+inc) , 
+  chip ( 6+inc ) , chip ( 7+inc ) 

WRITE  (30,44)  chip(13+inc) , chip(12+inc) , chip(23+inc) , 
+  chip(21+inc) ,chip(24+inc) 

WRITE  (30,45)  chip(14+inc) ,chip(20+inc) ,chip(25+inc) , 
+  chip(22+inc) , chip(19+inc) , chip(18+inc) 

if  ((idev  .eq.  1)   .or.   (inc  . ge .  20))  go  to  283 

write  (30,47) 

inc=20 

go  to  538 
477        format (lx,  Al) 

C  For  the  row  and  column  -  row  #1 

41  format (5x, Al , lx , 13 , 14 , lx , F5 . 2 , 6x, F5 . 2 , 6x, F5 . 2 , 6x , F5 . 2) 
C  For  row  #2 

42  format (15x, lelO . 3 , lx, lelO . 3 , lx, lelO . 3 , 2x, F5 . 3) 
C  For  row  #3 

43  format ( 16x , F5 . 3 , 3x , F6 . 1 , 7x , lelO . 3 , lx , F6 . 3 , 2x , F6 . 2 ) 
C  For  the  SUXES  parameters  -  row  #4 

44  FORMAT  ( 16X, F5 . 3 , 6x , F5 . 3 , 3x , F6 . 1 , 5x , F8 . 3 , 3x , lelO . 3) 
C  For  the  last  row 

45  format  (16x,   F5 . 3 , 3x , F6 . 1 , 7x , lelO . 3 , F8 . 4 , F8 . 3 , 4x , F6 . 2) 
C  Similar  to  41 

46  format  (5x,Al,2x,  F5 . 2 ,  2x,  F5 . 2 ,  6x , F5 . 2 , 6x, F5 . 2 , 5x, F6 . 2) 

47  format  (lx, '    ' ) 

49  format  (8x,  F5 . 2 ,  2x,  F5 . 2 ,  6x, F5 . 2 , 6x, F5 . 2 , 5x, F6 . 2) 

66  format ( 5x , Al , lx , I 3 , 14 , lx , F5 . 2 , 6x , F5 . 2 , 6x , F5 . 2 , 6x , F5 . 2 ) 

67  format(15x,F5.2,6x,F5.2,5x,F6.2) 

6  9  format ( 2x , Al , 2x , 1 2 , lx , 12 , lx , F5 . 2 , lx , F5 . 2 , lx , F5 . 2 , lx , 

+  F5.2,lx,F5.2,lx,l3) 
C    -- 

c 

C  ignore  the  averages  from  last  time 

C 

588  10=0 

do  578  10=1,24 

read  (20,477)  I DUMMY 
578  continue 

go  to  501 

C  -    

700  nodata=l 
GO  TO  707 
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283        if  (MARK  .eq.   '-')  go  to  537 
C 

708        if  (MARK  .eq.    '+')  N=N+1 
NM=N 

707        CALL  WAFER ( imark, n, i. chip, wave, wsum,alave,alsum,nodata) 

if  (nodata  .eq.  1)  go  to  937 
C  

c 

C  fill  the  wafer  map  files 

C 

C  correct  row  and  column  ? 

C 

537        if  (IWAFCUR  .NE.   iw)  GO  TO  168 

if  (map(jrow, jcol)   .eq.  map(irow, icol) )  go  to  184 
idump=map(irow, icol) -map (j row, jcol) 

C 

C      fill  with  zeros 
C 

jfill=l 

182  ifill-0 
zero=0 . 0 

do  183  ifill=l,idump 
izer=0 

do  186  izer=l,47 
write(48+izer , 189)  zero 

186  continue 

183  continue 

if(jfill   .eq.   0)  go  to  168 

C 

184  jfill=0 

if  (mark  .eq.    '+')  go  to  187 

idump=l 

go  to  182 

C 

C  

187  INC=0 

DO  600  INC=1,45 

write(50+inc,189)  CHIP(INC) 
189        format  (lx,lel0.3) 
600  CONTINUE 

write  (49,189)  chip(46) 

write  (50,189)  chip(47) 

C 

C  -  -- 

c 

C      determine  the  next  jrow  and  jcol 
168  jrow=irow 
j  col=icol 

if((jrow  .eq.  1)   .and.   (jcol  .eq.  6))  go  to  170 

if (((jrow. eq. 2) .or. (jrow. eq. 8))   .and.   (jcol  .eq.  8))  go  to  171 

if  (jcol  .eq.  9)  go  to  172 

jcol=jcol+l 
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go  to  500 

170  jrow=j row+1 
jcol=2 

go  to  500 

171  if  (jrow  .eq.  2)  jcol=l 
if  (jrow  .eq.  8)  jcol=4 
j  row=j  row+1 

go  to  500 

172  jcol=l 

if  (jrow  .eq.   7)  jcol=2 
j  row=j  row+1 
GO  TO  500 


C 

C  fill  the  rest  of  the  data  files 

C 

937        IF  (IWAFCUR  .NE.   iw)  GO  TO  427 
zero=0 . 0 

idump=map(9 , 6) -map(krow,kcol) 
kfill=0 

do  910  kfill=l, idump 
ink=0 

do  908  ink=l,47 

write(48+ink, 189)  zero 

if(kfill  .ne.   idump)  go  to  908 

write(48+ink,281) 
281        format('  *') 
908  continue 
910  continue 
427  jrow=l 

j  co 1=4 

if(zfil  .ne.  0.0)  go  to  920 

if  ((nodata  .eq.   0)    .and.    (BUST  .ne.   1.0))  go  to  511 


c  

c 

c 
c 

Will  wafer  maps  be 

generated  ? 

909 

if  ( (NM*100./itotchip)   .gt.  perchip)  , 

write(96,807) 

!     No  maps 

807 

format(' $  K=0') 
go  to  911 

C 

810 

write(96,808) 

!  Generate 

808 
C 

formate  $  K=l') 

911 

continue 

END 

c 

C  FILENAME :     INTRO . SUB 

C 
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SUBROUTINE  INTRO  (chip .wave ,wsum, alave , alsum, perchip , 
+  itotchip) 

C 

dimension  chip(l),  wave(l),  wsum(l) 
dimension  alave(l),  alsum(l) 

C 

common  /plomi/  MARK,   IWAF(IOO) ,   lotnuml ,  lotnum2 
common  map (10, 10) 

C 

character*l  IWAF 
character*l  IW 
character*!  MARK 
character*4  lotnuml 
character*4  lotnum2 

C 

C  -  

c 

C  FILL  THE  MAP  DIMENSION  TO  NUMBER  THE  CHIPS 

C 


162 
C-- 


j  map=0 

do  160  jmap=4,6 
map ( 1 , j  map ) = j  map  -  3 
map (9 ,  jmap)=jmap+59 
continue 


160  continue 
j  map=0 

do  161  jmap=2,8 
mao ( 2 . i maD ) — i ma 


map  ( 2  ,  j  map )  =j  map+2 
map ( 8 , j  map ) =j map+54 
continue 


161  continue 
j  map=0 


read  (97,477)  I DUMMY 

read  (97,*)  perchip,  itotchip 

C 

C  READ  THE  DUMMY  DATA 

C 

M0=0 

DO  378  M0=l,13 

READ  (20,477)IDUMMY 
477        format  (lx,  Al) 
378  CONTINUE 
C 

mo=0 

do  478  mo=l,5 
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478 
C 

734 
C 


281 
908 
C 

C--. 
C 
C 
C 


read  (25,477)  I DUMMY 
continue 

read  (98,734)  lotnuml,  lotnum2 
format ( 2 A4) 

ink=0 

do  908  ink=l,47 
write(48+ink,281) 
format('  *') 
continue 


1 

2 

3 
4 
5 

1006 


WRITE  THE  HEADING 

WRITE  (30,1) 
format('  #') 
WRITE  (30,2) 

format('  *    FILENAME  =  RESULTS .DAT' ) 
WRITE  (30,4) 

WRITE  (30,3)  lotnuml,  lotnum2 
format('  *    Lot  Number  =  ',lx,2A4) 
WRITE  (30,4) 
format('  *') 
WRITE  (30,5) 

formate  *  Wl  R  C  Wn  Ln 
+'  Lp') 


Wp  ', 
PDIFF' ) 


666 

C-  - 
C 
C 
C 


WRITE  (30,1006) 

format  ('  *  MCKT  SCKT 

WRITE  (30,*)  '*' 
WRITE  (30,6) 

format  ( '  *  CGDO  CGSO  CGBO ' , 

+'  RMS  ') 

WRITE  (30,7) 

format  ('  *  LD  TOX  SUB  ', 

+  '  VTO  KP') 

WRITE  (30,8) 

format ( '  *  GAMMA  PHI  U0  ' , 

+'  DELTA  VMAX') 

WRITE  (30,9) 

format('  *  XJ  KAPPA  NFS  ', 

+  '  ETA  THETA  ERR') 

write(30,666) 

format ( '  ********************************** 
+  ' ********************** ' ) 


INITIALIZE  THE  DIMENSIONS 
INC=0 

DO  100  INC=1,47 


100 


CHIP(INC)=0.0 
WAVE (INC) =0.0 
WSUM(INC)=0.0 
ALAVE(INC)=0.0 
ALSUM(INC)=0.0 
100  CONTINUE 

C  

RETURN 
END 

c 

C  FILENAME :  NEWCHIP.SUB 

C 

C 

SUBROUTINE  NEWCHIP (N , I , IW, IROW , ICOL , CHIP , WAVE , WSUM , 
+  ALAVE , ALSUM , IDEV , BUST , NODATA , I MARK , ZFIL) 

C 

dimension  chip(l) 

dimension  wave(l),  wsum(l) ,   alave(l) ,  alsum(l) 
common  /plomi/  MARK,   IWAF(IOO) ,   lotnuml,  lotnum2 

C 

character*l  IWAF 
character*l  IW 
character*l  MARK 
character*4  lotnuml 
character*4  lotnum2 

C 

C  -  

c 

C  READ  THE  NEW  DATA  IS  IT  GOOD  OR  BAD  ?   (-  =  bad) 

C 

read(20,77)  IWAF(I) 
77  format(Al) 

READ  (20,*,END=512)  IROW , ICOL, CHIP(l) , CHIP(2) , CHIP(46) , CHIP(47) 

READ  (20,*,END=512)  CHIP(3) ,CHIP(4) 

M=I 

if(I   .le.   1)  go  to  513 
if(IWAF(I)   .ne.   iwaf(I-l))  zfil=1.0 
if((IWAF(I)    .ne.   iwaf(I-l))    .and.    (N  .eq.  0))  1=0 
if(IWAF(I)   .ne.   iwaf(I-l))  CALL  WAFER ( imark.n, i , chip , 
+  wave ,wsum, alave , alsum , nodata) 

C 

C  DOES  IT  CONTAIN  N-  or  P- channel  DATA  or  both  ? 

C 

513        read  (25 , 69 , end=512) iw, ir , ic ,xndum,yndum,xpdum,ypdum, zdum, idev 
69  format(2x,Al,2x,I2,lx,I2>lx,F5.2,lx>F5.2,lx,F5.2f 
+  lx,F5.2,lx,F5.2,lx,I3) 

if(I  .eq.  0)  i=i+l 

iwaf (I)=iwaf (m) 
C  write  (30 , 69) iw, ir , ic , xndum , yndum , zdum, idev 

C  write  (30,68)  mark,  iwaf(I),i 
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C68 


512 


515 
C 

516 


C 

C-- 

C 

517 


format (lx,  Al ,   lx,  Al,lx,I2) 
if  (iwaf(I)   .ne.   iw)  go  to  513 

if  ((ir  .eq.   irow)   .and.   (ic  .eq.   icol))  go  to  515 
go  to  513 

write  (6,*)   '  The  correct  chip  was  not  located  in  LINE.FIL-2' 

BUST=1 . 0 

RETURN 

inc=0 

if(idev  .eq.  2)  inc=20  !  just  p-channel 

read  (20,*,end=512)  chip ( 9+inc ), chip (8+inc) , chip (10+inc) , 
chip(6+inc) ,chip(7+inc) 
chip(8+inc)=chip(8+inc)*l . e4 

read  (20 , * , end=512)  chip(13+inc) , chip(12+inc) , chip(23+inc) , 
chip(21+inc) , chip(24+inc) 

read  (20 , * , end=512)  chip(14+inc) , chip(20+inc) , chip(25+inc) , 
chip(22+inc) ,chip(19+inc) 

READ  (20,*,END=512)  chip ( 15+inc ), chip (16+inc) , chip (17+inc) , 
chip(ll+inc) , trash, chip (18+inc) , trash 
if  ((idev  .eq.  1)   .or.   (inc  . ge .  20))  go  to  517 


just  n- channel 


inc 

=20 

go 

to  516 

IFLG=0 

IF 

(chip(3)  . 

if 

(IFLG  .GE. 

chip (5)  =  100 

if 

<(chip(5) 

if 

(IFLG  .GE. 

if 

(idev  .eq. 

if 

((chip(6) 

if 

(IFLG  .GE. 

if 

((chip(7) 

if 

(IFLG  .GE. 

if 

<(chip<8) 

if 

(IFLG  .GE. 

if 

((chip(9) 

if 

(IFLG  .GE. 

if 

((chip(10) 

if 

(IFLG  .GE. 

if 

((chip(ll) 

if 

(IFLG  .GE. 

if 

((chip(12) 

if 

(IFLG  .GE. 

if 

(<chip(13) 

if 

(IFLG  .GE. 

if 

((chip(14) 

1) 


0.0)  IFLG=1 
go  to  607 


!  3  = 


=  100.*(ABS(chip(3)   -  chip(4) ) )/chip(3) 


!  5 


MCKT 

!  4  =  SCKT 

IFLG=2 
PDIFF 


2)  go  to  556 
.LE.  0.3)  .OR 
1)  go  to  607 

.LE.  1.0)  .OR.  (chip(7)  .GE. 
1)  go  to  607 

.LE.  250.0)  .OR.  (chip(8)  .GE 
1)  go  to  607 

.LE.  0.0)   .OR.   (chip(9)  .GE. 
1)  go  to  607 

.LE.  1.0E13)   .OR.  (chip(10) 
1)  go  to  607 


!  just  p-channel  6  =  NVTO 
(chip(6)   .GE.   1.5))  IFLG=3 


i 


.LT.  0.0)  .OR. 
1)  go  to  607 

.LE.  0.3)  .OR. 
1)  go  to  607 
.LE.  0.1)  .OR. 
1)  go  to  607 
.LE.  0.0)  .OR. 


(chip (11) 

(chip(12) 
(chip(13) 
(chip (14) 


7  =  NKP 
100.0))  IFLG=4 
!  8  =  NTOX 
1000.0))  IFLG=5 
!   9  =  NLD 
.0))  IFLG=6 

!   10  =  NSUB 
GE.  1.0E17)) 
!   11  =  NRMS 
.GT.  0.20))  IFLG=7 
!   12  =  NPHI 


IFLG=7 


GE.  0.9))  IFLG=8 
!   13  =  NGAMMA 

GE.  3.0))  IFLG=9 
!  14  =  NXJ 

GE.   1.5))  IFLG=10 
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if 
if 
if 
if 
if 
if 
if 
if 
if 

if 
if 
if 
if 
if 
if 
if 
if 
if 
if 
if 
if 
if 
if 

if 
if 
if 
if 
if 
if 
if 
if 
if 
if 
if 
if 
if 


IFLG  .GE 
(chip (15 
IFLG  .GE 
(chip (16 
IFLG  .GE 
(chip (17 
IFLG  .GE 
(chip(18 
IFLG  .GE 

(chip(19 
IFLG  .GE 
(chip(20 
IFLG  .GE 
(chip(21 
IFLG  .GE 
(chip(22 
IFLG  .GE 
(chip(23 
IFLG  .GE 
(chip (24 
IFLG  .GE 
(chip (25 
IFLG  .GE 

idev  .eq 
(chip(26 
IFLG  .GE 
(chip(27 
IFLG  .GE 
(chip(28 
IFLG  .GE 
(chip (29 
IFLG  .GE 
(chip(30 
IFLG  .GE 
(chip(31 
IFLG  .GE 


1)  go  to  607 
.LE.  0.0)  .OR. 
1)  go  to  607 
.LE.  0.0)  .OR. 
1)  go  to  607 
.LE.  0.0)  .OR. 
1)  go  to  607 
.LT.  0.0)  .OR. 
1)  go  to  607 


!   15  =  NCGDO 
(chip(15)   .GE.  1.0E-7))  IFLG=11 

!  16  =  NCGSO 
(chip(16)   .GE.  1.0E-7))  IFLG=12 

!  17  =  NCGB0 
(chip(17)   .GE.  1.0E-7))  IFLG=13 

!   18  =  NERR 
(chip (18)   .GT.  10.0))  IFLG=13 

!   19  =  NTHETA 


.LE.  -2.0)  .OR.  (chip(19)  .GE.  2.0)) 
1)  go  to  607  !  20  = 

.LE.  0.0)  .OR.  (chip(20)  .GE.  999.0) 
1)  go  to  607  !  21  = 

.LE.  -100.0)  .OR.  (chip(21)  .GE.  100 
1)  go  to  607  !  22  = 

.LE.  -2.0)  .OR.  (chip(22)  .GE.  2.0)) 
1)  go  to  607  !  23  = 

.LE.  0.1)  .OR.  (chip(23)  .GE.  1250.0 
1)  go  to  607  !  24  = 

.LE.  1.0e4)  .OR.  (chip(24)  .GE.  9.9e 
1)  go  to  607  !  25  = 

.LE.  0.0)  .OR.  (chip(25)  .GE.  4.0el4 
1)  go  to  607  !  26  = 


IFLG=14 

NKAPPA 
)  IFLG=15 

NDELTA 
.0))  IFLG=16 

NETA 

IFLG=17 

NUO 
))  IFLG=18 

NVMAX 
7))  IFLG=19 

NNFS 

))  IFLG=20 
PVTO 


1)  go  to  557 

.LE.  -3.0)  .OR.  (chip(26)  .GE. 
1)  go  to  607 

.LE.   1.0)    .OR.    (chip(27)  .GE. 
1)  go  to  607 

.LE.  250.0)   .OR.   (chip(28)  .GE 


!  just  n- channel 

-0.1))  IFLG=21 
!  27  =  PKP 
100.0))  IFLG=22 
!  28  =  PTOX 
.  1000.0))  IFLG=23 


1)  go  to  607 
.LE.   0.0)  .OR. 
1)  go  to  607 
.LE.  1.0E13)  . 
1)  go  to  607 
.LT.  0.0)  .OR. 
1)  go  to  607 


(chip(29)  .GE 
OR.  (chip(30) 
(chip(31)  .GT 


!  29  =  PLD 
5.0))  IFLG=24 
!  30  =  PSUB 
GE.  1.0E17))  IFLG=25 
!  31  =  PRMS 
0.20))  IFLG=25 
!  32  =  PPHI 


if 

((chip(32) 

.LE.  0.3)  .OR. 

(chip(32) 

.GE. 

0. 

9)) 

IFLG=26 

if 

(IFLG  .GE. 

1)  go  to  607 

i 

33 

=  PGAMMA 

if 

((chip(33) 

.LE.  0.1)  .OR. 

(chip(33) 

.GE. 

3. 

0)) 

IFLG=27 

if 

(IFLG  .GE. 

1)  go  to  607 

; 

34 

=  PXJ 

if 

((chip(34) 

.LE.  0.0)  .OR. 

(chip (34) 

.GE. 

1. 

5)) 

IFLG=28 

if 

(IFLG  .GE. 

1)  go  to  607 

1 

35 

=  PCGDO 

if 

((chip(35) 

.LE.  0.0)  .OR. 

(chip(35) 

.GE. 

1. 

0E- 

7))  IFLG=29 

if 

(IFLG  .GE. 

1)  go  to  607 

; 

36 

=  PCGSO 

if 

((chip(36) 

.LE.  0.0)  .OR. 

(chip(36) 

.GE. 

1. 

0E- 

7))  IFLG=30 

if 

(IFLG  .GE. 

1)  go  to  607 

i 

37 

=  PCGBO 

if 

((chip(37) 

.LE.  0.0)  .OR. 

(chip(37) 

.GE. 

1. 

0E- 

7))  IFLG=31 

if 

(IFLG  .GE. 

1)  go  to  607 

i 

38 

=  PERR 

if 

((chip(38) 

.LT.  0.0)  .OR. 

(chip(38) 

.GT. 

10.0))  IFLG=31 
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if  (IFLG  .GE. 


1)  go  to  607 


!   39  =  PTHETA 


if 

((chip(39) 

if 

(IFLG  .GE. 

if 

((chip(40) 

if 

(IFLG  .GE. 

if 

((chip(41) 

if 

(IFLG  .GE. 

if 

((chip(42) 

if 

(IFLG  .GE. 

if 

((chip(43) 

if 

(IFLG  .GE. 

if 

((chip(44) 

if 

(IFLG  .GE. 

if 

((chip(45) 

if 

(IFLG  .GE. 

.LE.  -10.0)  .OR.  (chip(39)  .GE.  10.0 
1)  go  to  607  !  40  = 

.LE.  0.0)  .OR.  (chip(40)  .GE.  999.0) 
1)  go  to  607  !  41  = 

.LE.  -100.0)  .OR.  (chip(41)  .GE.  100 
1)  go  to  607  !  42  = 

.LE.  -10.0)  .OR.  (chip(42)  .GE.  10.0 
1)  go  to  607  !  43  = 

.LE.  0.1)  .OR.  (chip(43)  .GE.  1250.0 
1)  go  to  607  !  44  = 

.LE.  1.0e4)  .OR.  (chip(44)  .GE.  9 . 9e 
1)  go  to  607  !  45  = 

.LE.  l.OelO)  .OR.  (chip(45)  .GE.  1.0 
1)  go  to  607 


))  IFLG=32 

PKAPPA 
)  IFLG=33 
PDELTA 
.0))  IFLG=34 
PETA 
))  IFLG=35 

PU0 
))  IFLG=36 

PVMAX 
7))  IFLG=37 

PNFS 
el4))  IFLG=38 


C 

C-.. 
C 

557 

607 
C 

G-  -  ■ 


MARK  =  '+' 
go  to  968 
MARK  =  ' - ' 
write  (30,*)   ' IFLG=' ,  IFLG 


968  continue 

C  write(30,9999)  MARK 

C9999  format(lx,Al) 

RETURN 

END 

c 

C  FILENAME:     WAFER. SUB 


SUBROUTINE  WAFER  ( IMARK , N , I , CHIP , WAVE , WSUM , ALAVE , 
+  ALSUM , NODATA) 

dimension  chip(l),  wave(l),  wsum(l) ,  alave(l),  alsum(l) 
common  /plomi/  MARK,  IWAF(100) ,  lotnuml,  lotnum2 

character*l  IWAF 
character*l  MARK 
character*4  lotnuml 
character*4  lotnum2 


DETERMINE  THE  WAFER  AVERAGES 
IF  (N  .EQ.  0)  RETURN 

IF  ((NODATA  .EQ.  1)  .and.  (I  .LE.  2))  RETURN 
IF  (NODATA  .EQ.   1)  GO  TO  791 
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IF  (I   .EQ.   1)  GO  TO  201 
IF  (IWAF(I)    . NE.   IWAF(I-l))  GO  TO  309 
201        if  (imark  .le.  0)  go  to  309 
J=0 

DO  200  J=l,47 

WSUM ( J ) =WSUM ( J ) +CHI P ( J ) 

WAVE(J)=WSUM(J)/N 
200  CONTINUE 
C 

C  

309 

791 


806 
C 


821 


'***  WAFER  AVERAGES  ***' 


Wl 


Wn 


MCKT 


Ln 


SCKT 


Wp 
PDIFF' 


if  (I  .EQ.   1)  RETURN 
if  (IWAF(I)    .eq.   IWAF(I-l))  RETURN 

WRITE  (30,4) 

WRITE  (30,*) 

WRITE  (30,4) 

WRITE  (30,4) 

WRITE  (30,*) 

+  '  Lp  ' 

WRITE  (30,*) 

WRITE  (30,*)  '*' 

WRITE  (30,6) 

WRITE  (30,7) 

WRITE  (30,8) 

WRITE  (30,9) 

WRITE  (30,806)  N,  1-1 

format  (lx,'  Num  chips/Tot  chips  =' , 12 , ' /' , 12) 

WRITE  (30,46)IWAF(I-1) ,WAVE(1) ,WAVE(2) ,WAVE(46) ,WAVE(47) 
WRITE  (30,48)  WAVE(3) ,WAVE(4) ,WAVE(5) 
WRITE  (30,*)   '  ' 
inc=0 

if  (idev  .eq.   2)   inc=20  !  just  p- channel 

write  (30,42)  wave ( 15+inc ), wave (16+inc) , wave (17+inc) , 
+  wave(ll+inc) 

write  (30,43)  wave ( 9+inc ), wave (8+inc) , wave (10+inc) , 
+  wave ( 6+inc ), wave (7+inc) 

WRITE  (30,44)  wave(13+inc) ,wave(12+inc) ,wave(23+inc) , 
+  wave(21+inc) ,wave(24+inc) 

WRITE  (30,45)  wave(14+inc) ,wave(20+inc) ,wave(25+inc) , 
+  wave(22+inc) , wave ( 19+inc) , wave ( 18+inc) 

if  ((idev  .eq.   1)   .or.   (inc  . ge .  20))  go  to  972 

write  (30,47) 

inc=20 

go  to  821 


THE  FORMAT  CARDS 

format('  *') 
format  ( '  * 
+  '  RMS  ') 

format  ( '  * 


CGDO 
LD 


CGSO 
TOX 


CGBO' , 
SUB  '  , 
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+  '  VTO  KP') 

8  format('  *  GAMMA  PHI  UO 
+'            DELTA  VMAX') 

9  format('  *  XJ  KAPPA  NFS 
+  '            ETA            THETA  ERR') 

C  For  row  #2 

42  format ( 15x , lelO . 3 , lx , lelO . 3 , lx , lelO . 3 , 2x , F5 . 3 ) 
C  For  row  #3 

43  format ( 16x , F5 . 3 , 3x , F6 . 1 , 7x , lelO . 3 , lx , F6 . 3 , 2x , F6 . 2 ) 
C  For  the  SUXES  parameters  -  row  #4 

44  FORMAT  (16X, F5 . 3 , 6x , F5 . 3 , 3x, F6 . 1 , 5x, F8 . 3 , 3x, lelO . 3) 
C  For  the  last  row 

45  format  (16x,  F5 . 3 , 3x, F6 . 1 , 7x, lelO . 3 , F8 . 4, F8 . 3 ,4x, F6 . 2) 
C  Similar  to  41 

46  format  (5x,Al,9x,  F5 . 2 ,  6x,  F5 . 2 ,  6x,F5.2,6x,F5.2) 

47  format  (lx, '   ' ) 

48  format  (15x,  F5.2,  6x,  F5.2,  8x,  F6.2) 

C   

972  CONTINUE 

write(30,666) 

666        format ('  ft*************************************************' 
+  '  '  ) 

CALL  LOTS(lot,n, i ,wave ,wsum, alave , alsum, idev.nodata) 

RETURN 

END 

c 

C  FILENAME :     LOTS . SUB 

C 

SUBROUTINE  LOTS ( lo t , n , i , wave , wsum , alave , alsum , idev , nodata) 

C 

dimension  wave(l),  wsum(l),  alave(l),  alsum(l) 
common  /plomi/  MARK,  IWAF(IOO) ,  lotnuml,  lotnum2 

C 

character* 1  IWAF 
character*l  MARK 
character*4  lotnuml 
character*4  lotnum2 

C 
C 

C  

c 

C  DETERMINE  THE  LOT  AVERAGES 

C 

lot=lot+l 
K=0 

DO  300  K=l,47 

ALSUM(K)=ALSUM(K)+WAVE(K) 

ALAVE ( K) =ALSUM(K) /LOT 
300  CONTINUE 
C 
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C  REINITIALIZE  THE  WAFER  DIMENSIONS 

C 

N=0 
1=0 
INC=0 

DO  1000  INC=1,47 
WAVE (INC) =0.0 
WSUM(INC)=0.0 
1000  CONTINUE 

IF  (LOT  .LT.   2)  RETURN 
IF  (nodata  .eq.   0)  RETURN 


C  SUM  UP  THE  RESULTS  FOR  THE  LOT 

C 

WRITE  (30,4) 

WRITE  (30,*)   '***  LOT  AVERAGES  ***' 
WRITE  (30,4) 
WRITE  (30,4) 

WRITE  (30,*)   '*    Lot  Wn  Ln  Wp  ' , 

+  '  Lp  ' 

WRITE  (30,*)    '*  MCKT  SCKT  PDIFF' 

WRITE  (30,*)  '*' 
WRITE  (30,6) 
WRITE  (30,7) 
WRITE  (30,8) 
WRITE  (30,9) 
WRITE  (30,836)  lot 
836        format  (lx, '  Num  wafers  =',I2) 

WRITE  (30,50)  lotnuml,lotnum2, alave (1) ,alave(2) ,alave(46) , 
+  alave(47) 

WRITE  (30,48)  alave(3) , alave(4) , alave(5) 

WRITE  (30,*)   '  ' 

inc=0 

if  (idev  .eq.  2)  inc=20  !  just  p- channel 

936  write  (30,42)  alave(15+inc) , alave(16+inc) , alave(17+inc) , 
+  alave(ll+inc) 

write  (30,43)  alave(9+inc) , alave(8+inc) , alave(10+inc) , 
+  alave ( 6+inc ) , alave ( 7+inc ) 

WRITE  (30,44)  alave ( 13+inc ), alave (12+inc) , alave (2 3+inc) , 
+  alave ( 2 1+inc ), alave (24+inc) 

WRITE  (30,45)  alave ( 14+inc ), alave (20+inc) , alave (2 5+inc) , 
+  alave (22+inc) , alave (19+inc) , alave (18+inc) 

if  ((idev  .eq.   1)   .or.   (inc  .ge.  20))  go  to  937 

write  (30,47) 

inc=20 

go  to  936 

C  .-  

937  CONTINUE 
write(30,666) 

666        format ('  *************************************************** , 
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4-    '  '  ) 


c 
c 
c 

4 
6 

7 

8 

9 

C 

42 
C 

43 
C 

44 
C 

45 
47 
48 
50 
C-- 


THE  FORMAT  CARDS 
format('  *') 

format  ('  *  CGDO  CGSO  CGBO' , 

+'  RMS  ') 

format  ( '  *  LD  TOX  SUB  '  , 

+ '  VTO  KP ' ) 

formatC  *  GAMMA  PHI  UO  ', 

+'  DELTA  VMAX') 

formatC  *  XJ  KAPPA  NFS 

+  '  ETA  THETA  ERR') 

For  row  #2 

format ( 15x , lelO . 3 , lx , lelO . 3 , lx , lelO . 3 , 2x , F5 . 3 ) 
For  row  #3 

formatCiex.FS.S.Sx.Fe.l^x.lelO.S.lx.Fe.S^x.FS^) 
For  the  SUXES  parameters  -  row  #4 

FORMAT  ( 16X , F5 . 3 , 6x , F5 . 3 , 3x , F6 . 1 , 5x , F8 . 3 , 3x , lelO . 3 ) 
For  the  last  row 

format  (16x,  F5 . 3 , 3x, F6 . 1 , 7x, lelO . 3 , F8 . 4 , F8 . 3 , 4x , F6 . 2) 
format  (lx, '   ' ) 

format  (15x,  F5.2,  6x,  F5.2,  8x,  F6.2) 

format  (5x,2A4,2x,  F5.2,  6x,  F5 . 2 ,  6x, F5 . 2 , 6x, F5 . 2) 


return 
END 
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RANKEY.  FOR 


C 

C  FILENAME :     RANKEY . FOR 

C 

C  PURPOSE:     THIS  PROGRAM  ORDERS  THE  CORRELATION  COEFFICIENT  DATA 

C  AND  DETERMINES  WHICH  WAFER  MAPS  TO  MAKE. 

C 

dimension  numfile(50) ,  filnl(50),  filn2(50) 
dimension  number ( 50 ) ,  cor(50),  nchips(50) 
dimension  nummax(50) ,  cormax(50) ,  nchipmax(50) 

C 

character*l  iwaf 
character*4  filnl 
character*4  filn2 
character*4  filela 
character*4  file2a 
character*4  filelb 
character*4  file2b 
character*4  lotnuml 
character*4  lotnum2 

C 

open(unit=10 ,  f ile=' in. dat ' ,  status=' old' ) 
open(unit=ll ,  f ile= ' lotnumber . fil ' ,  status=' old' ) 
open(unit=12 ,   f ile='nlw. fil' ,   status=' old' ) 
open(unit=15 ,   file='coref .fil' ,  status=' old' ) 
open(unit=20 ,   f ile=' out . dat ' ,   status= ' unknown' ) 
open(unit=30 ,   f ile=' out2 . dat ' ,   status=' unknown' ) 
open(unit=40 ,   f ile=' mapl . dat ' ,   status=' old' ) 

C  open(unit=50 ,  f ile=' out4 . dat ' ,  status= ' unknown' ) 

C 

C  read  the  lotnumber  and  wafer  letter 

C 

read  (11,6)  lotnuml,  lotnum2 

6  format (2A4) 
read  (12,7)  iwaf 

7  format (Al) 
C 

C  read  the  coref  file 

C 

j=0 

do  131  j=2,46 

read  (15 , 132 , end=733)  numfile(j),  filnl(j),  filn2(j) 
132        format  (I3.2A4) 
131  continue 
C 

C  read  the  mapl  file 

C 

733        read  (40,833)  filela,  filelb 
833        format  (2A4) 

read  (40,833)  file2a,  file2b 

C 
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C  read  the  input  data  file 

C 

j-0 

do  133  j-2,46 

read  (10,*,err=135,end=135)  number(j),  itrashl,  itrash2,  itrash3, 
+  itrash4,  itrash5,  itrash6,  itrash7,  cor(j),  nchips(j) 
133  continue 

135  jmax=j-l 
C 

C  order  the  input  data  according  to  the  correlation  coefficient 

C 

j=0 

nummax ( 1 ) =0 

do  126  j=2,jmax 

cormax(j )=0.0 

nummax (j )=0 
126  continue 
C 

j-0 

do  136  j=2,jmax 

if  (abs(cor(j))   .gt.  abs(cormax(2) ) )  nummax(2)=number ( j ) 
if  (abs(cor(j))   .gt.  abs(cormax(2) ) )  nchipmax(2)=nchips ( j ) 
if  (abs(cor(j))   .gt.  abs(cormax(2) ) )  cormax(2)=cor ( j ) 

136  continue 
C 

i=2 

137  i=i+l 

if  (i  .gt.  jmax)  go  to  800 
j-0 

do  138  j=2,jmax 

if  ( (abs(cor( j ) )   .le.  abs(cormax(i) ) )   .or.   (abs(cor(j))  .gt. 
+  abs(cormax(i-l) ) ) )  go  to  138 

if  ((number(j)   .eq.  nummax(i-l))   .or.   (number(j)  .eq. 
+  nummax(i-2) ) )  go  to  138 

if  (i  .It.  5)  go  to  602 

if  ((number(j)   .eq.  nummax(i-3))   .or.   (number(j)  .eq. 
+  nummax(i-4) ) )  go  to  138 
602        nummax(i)=number ( j ) 

nch  ipmax ( i ) =nch ip s ( j ) 
cormax ( i ) =cor ( j ) 

138  continue 
go  to  137 

C 

C  print  the  results  in  order 

C 

800        write  (20,701)  lotnuml,  lotnum2,  iwaf 

701  format(/, lx, 'Lot  ',2A4, 'Wafer  ' ,A1) 
write  (20,702) 

702  format(/,13x, 'Corr. ' ) 
write  (20,700) 

700        format  (lx,  'Parameter'  ,3x,  'Coef .  '  ,3x,  '//Chips'  ,/) 
C 
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L=0 
m=50 

filnl(50)='xxxx' 

n=0 

k=l 

802  k=k+l 

if  (k  .gt.  jmax)  go  to  900 

C 

C  find  the  proper  value  for  j 

C 

3=0 
i=0 

do  803  i=2,46 

if  (nummax(k)   .eq.  number (i))  j=i 

803  continue 
n=n+l 

C 

write  (20,801)  filnl(j),  filn2(j),  cormax(k) ,  nchipmax(k) 
801  format(2x,2A4,2x,F6.3,5x,I2) 

if  ((n  .eq.   1)    .and.    (filnl(j)    .eq.    ' MEAC ' ) )  m=j 

C 

C  DETERMINE  WHICH  MAPS  TO  GENERATE 

C 

if  (L  .GT.   1)  GO  TO  802 
if  (filnl(j)   .EQ.   'MEAC')  GO  TO  802 
if  (filnl(j)   .EQ.    'PDIF')  GO  TO  802 
if  (filnl(j)   .EQ.    'SIMC')  GO  TO  802 

if  ((filnl(m)   .eq.   'MEAC')   .and.   (filnl(j)   .eq.  filela)) 
+  go  to  802 

if  ((filnl(m)    .eq.    'MEAC')    .and.    (filnl(j)    .eq.   f ile2a) ) 
+  go  to  802 
C  write  (50,808)  filela,  filelb 

C  write  (50,808)  file2a,  file2b 

C  write  (50,808)  filnl(j),  filn2(j) 

C  write  (50,808)  filnl(m),  filn2(m) 

write  (30,808)  filnl(j),  filn2(j) 
808        format (2A4) 

L=L+1 

go  to  802 

C 

900  CONTINUE 
END 
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NLW.FIL 


Y 

0 
4 

2.46 
2.46 

RC . FIL 

* 

*  FILENAME  =  RC.FIL 

*  THIS  IS  A  FILE  OF  CHIPS  OF  WHICH  SUXES  AND  SPICE  ARE  TO  BE 

*  RUN 
* 

*  Wafer  Row  Column 
Y    2  4 

LINE. FIL 

* 

*  Filename  =  LINE. FIL 

*  Xerox  wafer 

*  Wafer  Row    Column    Wn      Ln      Wp      Lp     (Measured  Freq)  NorPorBOTH 

*  FORMAT  (2x,Al,2x,I2>lx)I2,lx,F5.2,lx,F5.2,lx,F5.2,lx,F5.2,lx)F5.2,lx,I3) 
Y      2    4    2.46     2.46    6.46    2.46    9.40  3 

LOTNUMBER . FIL 

YY027 

NPARIN . FIL 

TOX  LD 

635. 

.38 

PPARIN.FIL 

TOX  LD 

635. 

.45 
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NINPUT.XXX 


* 

* 
* 

* 
* 
* 


FILENAME :     NINPUT . DAT 

PURPOSE:     THIS  IS  THE  N-CHANNEL  INPUT  DATA  FOR  SPICE 

WHICH  GIVES  AN  INITIAL  GUESS  AND  THE  RANGE  IN 
WHICH  TO  EXPECT  THE  PARAMETER  VALUE. 
SUXES  IS  RESTRICTED  TO  THIS  RANGE. 
Date  generated  6-20-83 


parameters 
Level 
Type 
Xld 
Tox 
Nsub 
Vto 
Xkp 
Gamma 
Phi 
Uo 

Uexp 

Ucrit 

Delta 

Vmax 

Xj 

Lamda 

Kappa 

Nfs 

Neff 

Nss 

Tgate 

Eta 

Theta 


initial 
3 


1.0 

uuu 

TTT 
SSS 
RRR 
OOOd-6 
MMM 
JJJ 
KKK 
0.0 
0.0 
2.685 
7.9d5 
HHH 
0.0 
.8000 
2.28dl2 
0.0 
lell 
1.0 
.1306 
.3612 


lowbound 
1 

-1.0 

0.  1 
250.0 

1.  dl3 
0.3 
1.0d-i 

.1 
.3 


0.1 

0.  001 

1.  d4 
•l.d2 

l.d4 

.1 
l.d-6 
0.0 

0.  0 

1.  d-2 
0.0 
0.0 

•2.0 
■2.0 


upperbound 
3 

1.0 
1.0 
1000.0 
l.dl7 
1.5 
100. d-6 
2.0 
.9 

1250.0 
2.3 
l.d6 
l.d2 
9.9d7 
1.5 
l.d-1 
999.0 
4el4 
l.dl 
lel4 
1.0 
2.0 
2.0 
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PINPUT.XXX 


* 
* 
* 
* 
* 

* 
* 
* 


FILENAME :     PINPUT . DAT 

PURPOSE:     THIS  IS  THE  P-CHANNEL  INPUT  DATA  FOR  SPICE 

WHICH  GIVES  AN  INITIAL  GUESS  AND  THE  RANGE  IN 
WHICH  TO  EXPECT  THE  PARAMETER  VALUE. 
SUXES  IS  RESTRICTED  TO  THIS  RANGE. 
Date  generated  6-20-83 


parameters    initial      lowbound  upperbound 


T  1 

Level 

3 

1 

3 

Type 

-1.0 

-1. 

0 

1.0 

Xld 

UUU 

0. 

1 

2.0 

Tox 

TTT 

250. 

0 

1000.0 

Nsub 

SSS 

1. 

dl3 

1.417 

Vto 

RRR 

-3. 

0 

-0.1 

Xkp 

OOOd-6 

1. 

Od-6 

100. d-6 

Gamma 

MMM 

1 

2.0 

Phi 

JJJ 

3 

.9 

Uo 

KKK 

0. 

1 

1250.0 

Uexp 

0.0 

0. 

001 

2.3 

Ucrit 

0.0 

1. 

d4 

l.d6 

Delta 

.6749 

-1. 

d2 

l.d2 

Vmax 

3 . 24d5 

1. 

d4 

9.9d7 

Xj 

HHH 

1 

1.5 

Lamda 

0.0 

1. 

d-6 

l.d-1 

Kappa 

.20 

0. 

0 

999.0 

Nfs 

1.17el2 

lelO 

4el4 

Neff 

0.0 

1. 

d-2 

l.dl 

Nss 

lell 

0. 

0 

lel4 

Tgate 

-1.0 

-1. 

0 

1.0 

Eta 

.1210 

-10.0 

10.0 

Theta 

.2396 

-10.0 

10.0 
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NOPT.FIL 


FILENAME :     NOPT . FIL 

PURPOSE:     THIS  IS  AN  OPTIONS  FILE  FOR  AN  N-CHANNEL  TRANSISTOR 

IN  WHICH  THE  DESIRED  BIASES  ARE  SELECTED  FROM  THE  IV  DATA 
Date  generated  6-20-83 


Criteria 

Nsign=4  Maxfn=200  Delta=1.0d-9  Epsil=1.0d-6 
Linit=0.1  Lscal=2.0  Luppr=1000  Fcdsw=0.1 


RangVD 
RangVG 
RangVB 
Weight 


min=0 . 00  max=5 . 5  incr=0 . 0 
min=0 . 00  max=5 . 5  incr=0 . 0 
min=0 . 0        max=0 . 0      incr=l . 0 


POPT.FIL 


* 
* 
* 
* 
* 
* 
* 
* 


* 
* 
* 


FILENAME : 
PURPOSE: 


POPT.FIL 

THIS  IS  AN  OPTIONS  FILE  FOR  A  P-CHANNEL  TRANSISTOR 

IN  WHICH  THE  DESIRED  BIASES  ARE  SELECTED  FROM  THE  IV  DATA 

Date  generated  6-20-83 


Criteria 

Nsign=4  Maxfn=200  Delta=1.0d-9  Epsil=1.0d-6 
Linit=0.1  Lscal=2.0  Luppr=1000  Fcdsw=0.1 


RangVD 
RangVG 
RangVB 
Weight 


min=-5.5  max=0 . 00  incr=0.0 
min=-5.5  max=0.00  incr=0.0 
min=0 . 0        max=0 . 0      incr=l . 0 
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NSTRAT . FIL 


* 
* 


FILENAME:     NSTRAT. FIL 

PURPOSE:     THIS  IS  THE  STRATEGY  FILE  TO  HELP  ALLEVIATE  THE 
REDUNDANT  PARAMETER  PROBLEM. 
Date  generated  6-20-83 


parameters 
Level 
Type 
Xld 
Tox 
Nsub 
Vto 
Xkp 
Gamma 
Phi 
Uo 

Uexp 

Ucrit 

Delta 

Vmax 

Xj 

Lamda 

Kappa 

Nfs 

Neff 

Nss 

Tgate 

Eta 

Theta 


first 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
1 
0 
0 
0 
0 
0 
0 
0 
0 

1 
1 


second  third  fourth 


0 
0 
0 
1 
0 
1 
0 

1 

0 
0 
0 
0 

1 

0 
0 
0 
0 

1 

0 
0 
0 
0 
0 


0 
0 

1 

0 

1 

0 

1 

0 
0 
0 
0 
0 
0 
0 

1 

0 
0 
0 
0 
0 
0 
0 
0 


0 
0 
0 
0 
0 
0 
0 
0 

1 
1 

0 
0 
0 

1 

0 
0 

1 

0 
0 
0 
0 
0 
0 


fifth 
0 
0 
0 

1 

0 

1 

0 

1 

0 
0 
0 
0 

1 

0 
0 
0 
0 

1 

0 
0 
0 
0 
0 
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PSTRAT.FIL 


* 
* 
* 
* 
* 
* 
* 
* 
* 


FILENAME :     PSTRAT . FIL 

PURPOSE:     THIS  IS  THE  STRATEGY  FILE  TO  HELP  ALLEVIATE  THE 
REDUNDANT  PARAMETER  PROBLEM. 
Date  generated  6-20-83 


parameters 
Level 
Type 
XI  d 
Tox 
Nsub 
Vto 
Xkp 
Gamma 
Phi 
Uo 

Uexp 

Ucrit 

Delta 

Vmax 

Xj 

Lamda 

Kappa 

Nfs 

Neff 

Nss 

Tgate 

Eta 

Theta 


first 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
1 
0 
0 
0 
0 
0 
0 
0 
0 

1 

1 


second 
0 
0 
0 
1 
0 
1 
0 

1 

0 
0 
0 
0 

1 

0 
0 
0 
0 

1 

0 
0 
0 
0 
0 


third 
0 
0 
1 
0 
1 
0 
1 
0 
0 
0 
0 
0 
0 
0 

1 

0 
0 
0 
0 
0 
0 
0 
0 


fourth 
0 
0 
0 
0 
0 
0 
0 
0 
1 
1 
0 
0 
0 

1 

0 
0 

1 

0 
0 
0 
0 
0 
0 


fifth 
0 
0 
0 
1 
0 

1 

0 

1 

0 
0 
0 
0 

1 

0 
0 
0 
0 

1 

0 
0 
0 
0 
0 


117 


IVCHARN.XXX 


* 

*  FILENAME  =  IVCHARN . FIL 

*  NOTE:  Use  with  MODEL. FIL  to  make  IVCHARN. CEL 

* 

*  N-CHANNEL  MOS  OUTPUT  CHARACTERISTICS 

Ml  1  2  0  0  MODN  L=BZBUM  W=FFFUM  AD=12 . 5P  AS=12 . 5P  PD=15U 

+  PS=15U  NRD=2.0  NRS=2.0 

* 

VD  3  0 
VG  2  0 
VIDS  3  1 
* 

.OPTIONS  NODE  NOPAGE 
.DC  VD  0  5   .5  VG  0  5  1 
.PLOT  DC  I (VIDS) 
.END 

IVCHARP.XXX 

* 

*  FILENAME  =  IVCHARP.FIL 

*  NOTE:     Use  with  MODEL. FIL  to  make  IVCHARP.CEL 

* 

*  P- CHANNEL  MOS  OUTPUT  CHARACTERISTICS 

* 

M2  1  2  0  0  MODP  L=BZBUM  W=FFFUM  AD=42.25P  AS=42.25P  PD=26U 
+  PS=26U  NRD=1.0  NRS=1.0 

* 

VD  4  0 
VG  2  0 
VIDS  4  1 

.OPTIONS  NODE  NOPAGE 

.DC  VD  0  -5  -.5  VG  0  -5  -1 

.PLOT  DC  I (VIDS) 

.END 
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NMODEL.XXX 


* 
* 

*  FILENAME  =  NMODEL.FIL 

*  NOTE:  This  file  uses  parameters  from  the  process  and  SUXES 
* 

.MODEL  MODN  NMOS  LEVEL=3  LD=sqsUM  TOX=rqrUM  NSUB=pqp  VTO=oqoV 
+  KP=nqnU  GAMMA=mqm  PHI=lql  UO=kqk  DELTA-jqj  VMAX=iqi 
+  XJ-hqhUM  KAPPA=gqg  NFS=fqf  ETA=eqe  THETA=dqd 
+  CJSW=0.6E-9  CGDO=cqc  CGSO=bqb  CGBO=aqa 

+  RSH-11.63  JS=5.0E-5  TPG=1.0  MJ=.5  MJSW=.5  PB=.9  XQC=.45 
* 

PMODEL.XXX 

* 
* 

*  FILENAME  =  PMODEL. FIL 

*  NOTE:  This  file  uses  parameters  from  the  process  and  SUXES 
* 

.MODEL  MODP  PMOS  LEVEL=3  LD=sqsUM  TOX=rqrUM  NSUB=pqp  VTO=oqoV 
+  KP=nqnU  GAMMA=mqm  PHI=lql  UO=kqk  DELTA=jqj  VMAX=iqi 
+  XJ=hqhUM  KAPPA=gqg  NFS=fqf  ETA=eqe  THETA=dqd 
+  CJSW=0.6E-9  CGDO=cqc  CGSO=bqb  CGBO=aqa 

+  RSH-61.57  JS=5.0E-5  TPG--1.0  MJ=.5  MJSW=.5  PB=.9  XQC=.45 
* 

SUBCKT.XXX 

* 

*  FILENAME  =  INVERT . FIL 

.SUBCKT  INVERT  12  3 
* 

Ml  100  1  0  0  MODN  L=hzhUM  W=iziUM  AD=12.5P  AS=12.5P  PD=15U 
+  PS=15U  NRD=2.0  NRS=2.0 

M2  100  1  3  3  MODP  L=jzjUM  W=kzkUM  AD=42.25P  AS=42.25P  PD=26U 

+  PS=26U  NRD=1.0  NRS=1.0 

Rl  100  2  1332 

CI  2  0  .043P 

* 

.ENDS  INVERT 
* 
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CIRCUIT. FIL 

* 

*  FILENAME  =  CIRCUIT. FIL 

*  NOTE:  Use  with  NMODEL. FIL,  PMODEL. FIL,  and  SUBCKT.XXX  to  make  CIRCUIT. CEL 
* 

XI  1  2  3  INVERT 
X2  2  4  3  INVERT 
X3  4  5  3  INVERT 
X4  5  6  3  INVERT 
X5  6  7  3  INVERT 
X6  7  8  3  INVERT 
X7  8  9  3  INVERT 
X8  9  10  3  INVERT 
X9  10  11  3  INVERT 
REND  11  200  1110 
CEND  200  0  .044P 
X10  200  12  3  INVERT 

XII  12  13  3  INVERT 
XI 2  13  14  3  INVERT 
XI 3  14  15  3  INVERT 
X14  15  16  3  INVERT 
X15  16  17  3  INVERT 
X16  17  18  3  INVERT 
X17  18  19  3  INVERT 
X18  19  20  3  INVERT 
X19  20  400  3  INVERT 
RCIR  400  1  334 
CCIR  1  0  .073P 

RPP  1  101  1.175 
CPP  101  0  .322P 
CPROBE  101  0  .IP 

*X20  6  50  30  51  52  OUTBUF3 
* 

.IC  V(l)=5  V(2)=0  V(4)=5  V(5)=0  V(6)=5  V(7)=0  V(8)=5 

.IC  V(9)=0  V(10)=5  V(11)=0  V(12)=5  V(13)=0  V(14)=5  V(15)=0 

.IC  V(16)=5  V(17)=0  V(18)=5  V(19)=0  V(20)=5 
*.IC  V(30)=5  V(50)=0  V(51)=0  V(52)=5 

.IC  V(3)=5 
* 

VCC  3  0  5V 
*VCC0  30  0  5.0V 

*VA  24  0  PULSE(0V  5V  ONS  5NS  5NS  1000NS  2000NS) 
* 

.OPTIONS  VNTOL=10UV  LIMPTS=500  ITL4=50  ITL5=0 

.TRAN  1.0NS  200NS  UIC 

*.PLOT  TRAN  V(101) 

*.PLOT  TRAN  V(l) 

*.PLOT  TRAN  V(2) 

.PLOT  TRAN  V(4) 

*.PLOT  TRAN  V(50) 

.END 
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PCHIPS . FIL 


*  percent  acceptable  chips  to  make  wafer  maps  ...  number  of  test  sites/wafer  30.0 
65 


STEND.FIL 

4 

2 

1 

1 
J. 

6 

8 

9 

Q 

1 

1 

1 

o 
z 

9 

9 

9 

Q 

o 

4 

0 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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COREF.FIL 


2  W 

3  L 

4  WP 

5  LP 

6  MCKT 

7  SCKT 

8  PDIFF 

9  NCGDO 

10  NCGBO 

11  NRMS 

12  NLD 

13  NTOX 

14  NSUB 

15  NVTO 

16  NKP 

17  NGAMMA 

18  NPHI 

19  NUO 

20  NDELTA 

21  NVMAX 

22  NXJ 

23  NKAPPA 

24  NNFS 

25  NETA 

26  NTHETA 

27  NERR 

28  PCGDO 

29  PCGBO 

30  PRMS 

31  PLD 

32  PTOX 

33  PSUB 

34  PVTO 

35  PKP 

36  PGAMMA 

37  PPHI 

38  PUO 

39  PDELTA 

40  PVMAX 

41  PXJ 

42  PKAPPA 

43  PNFS 

44  PETA 

45  PTHETA 

46  PERR 


EDNO.WAY 


!  filename  =  edno.way 

delete  l:end-l 
exit 

EDSUXN.WAY 


!  filename  =  EDSUXN.WAY 

s/77/N/l:end 

s/22/N22/l:end 

exit 

EDSUXP.WAY 


!  filename  =  EDSUXP.WAY 

s/77/P/l:end 

s/22/P22/l:end 

exit 


EDPAR.WAY 


!  FILENAME  =  EDPAR.WAY 

!  delete  unnecessary  rows 
1 

t  'X16V 
delete  1 : . -1 
t  'Uo' 

delete  .+1 : .+2 
t  'Xj' 
delete  .+1 
t  'Xnfs' 
delete  .+1 : .+3 


delete  first  column 


sub/ 

Xld 

//I 

:  end- 

1 

sub/ 

Tox 

//I 

:  end- 

•1 

sub/ 

Nsub 

//I 

:  end- 

1 

sub/ 

Vto 

//I 

:  end- 

1 

sub/ 

Xkp 

//I 

■  end- 

1 

sub/ 

Gamma 

//I 

end- 

1 

sub/ 

Phi 

//I 

end- 

1 

sub/ 

Uo 

//I 

end- 

1 

sub/ 

Delta 

//I 

end- 

1 

sub/ 

Vmax 

//I 

end- 

1 

sub/ 

Xj 

//I. 

end- 

1 

sub/ 

Xkappa 

//I: 

end- 

1 

sub/ 

Xnfs 

//I: 

end- 

1 

sub/ 

Theta 

//I: 

end- 

1 

sub/ 
i 

Eta 

//I: 

end- 

1 

resequence 

exit 

EDSPICEN.WAY 


!  FILENAME  =  EDSPICEN.WAY 

s/  //end- 5 : end 

s/777/N/l:end 

exit 

EDSPICEP.WAY 


!  FILENAME  =  EDSPICEP.WAY 

s/  //end- 5  rend 

s/777/P/l:end 

exit 

EDLAST.WAY 


!  filename  =  EDLAST.WAY 
t  'N=' 

move  end-1  to  . 

delete 

1 

t  'M=' 

move  end-1  to  . 

delete 

1 

t  'L=' 

move  end-1  to  . 
delete 
resequence 
exit 


EDNRMS . WAY 


!  filename  =  ednrms.way 

t  'again' 
delete  . : end- 1 
t  end- 50 
t  'optimization' 
.  -3 

delete  1: . -1 
delete  .+l:end-l 
resequence 
1 

s/max  err=// 

s/%  average  err=// 

s/%  rms=// 

exit 

EDPRMS . WAY 


!  filename  =  edprms.way 

t  'again' 
t  'again' 
delete  . : end-1 
t  end- 50 
t  'optimization' 
.  -3 

delete  1 : . -1 
delete  .+1: end-1 
resequence 
1 

s/max  err=// 

s/%  average  err=// 

s/%  rms=// 

exit 

EDRES . WAY 


!   FILENAME  =  EDRES. WAY 

!  add  an  '#'  to  denote  new  data 
1 

t  '$  N=' 
delete  . -2 : . 
copy  1  to  . 
resequence 
exit 
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EDSKEY. WAY 


FILENAME  -  EDSKEY. WAY 

delete  end-1 
! 

'  *' 
'*' 

'  *' 
'  *' 
'  *' 

delete  1 : . 
i 

resequence 
exit 


EDMAP . WAY 


!  FILENAME  =  EDMAP. WAY 
1 

t  'K=' 
delete 

move  end-1  to  . 
! 

1 

t  'N=' 
delete 

move  end-1  to  . 
I 

delete  end- 2: end-1 

resequence 

exit 

EDS IMF. WAY 


!  FILENAME  =  EDS IMF. WAY 
1 

t  'GET  7' 
t  ' SDB' 
t  'ENTRY' 
delete  1:+ 
t  'NUMBER' 
delete  . : end-1 
i 

S/-JAN-//1: end-1 
S/-FEB-//1: end-1 
S/-MAR-//1: end-1 
S/-APR-//1: end-1 
S/-MAY-//1: end-1 
S/-JUN-//1: end-1 
S/-JUL-//1: end-1 
S/-AUG-//1: end-1 
S/-SEP-//1: end-1 
S/-0CT-//1: end-1 
S/-N0V-//1: end-1 
S/-DEC-//1: end-1 
i 

s/://l: end-1 
s/://l: end-1 
• 

resequence 
exit 


EDPDIFF.WAY 


!  FILENAME  =  EDPDIFF.WAY 
1 

t  'GET  8' 
t  ' SDB' 


t  'ENTRY' 
delete  1:+ 
t  'NUMBER' 


delete 
■ 

. :end-l 

s/ 

-JAN- 

//I 

end- 

1 

s/ 

-FEB- 

//I 

end- 

1 

s/ 

-MAR- 

//I 

end- 

1 

s/ 

-APR- 

//I 

end- 

1 

*/ 

-MAY- 

//I 

end- 

1 

s/ 

-JUN- 

//I 

end- 

1 

s/ 

-JUL- 

//I 

end- 

1 

s/ 

-AUG- 

//I 

end- 

1 

s/ 

-SEP- 

//I 

end- 

1 

s/ 

-OCT- 

//I 

end- 

1 

s/ 

-NOV- 

//I 

end- 

1 

s/ 

-DEC- 

//I 

end- 

1 

1 

s/ 

://l: 

end- 

■1 

s/ 

://l: 

end- 

■1 

! 


resequence 
exit 


EDMEAF . WAY 


FILENAME  =  EDMEAF. WAY 


1 


t  'GET  6' 

t  '  SDB' 

t  'ENTRY' 

delete  1:+ 

t  'NUMBER' 

delete  . : end-1 

S/-JAN-//1 

end- 

1 

S/-FEB-//1 

end- 

1 

S/-MAR-//1 

end- 

1 

S/-APR-//1 

end- 

i 
i 

S/-MAY-//L 

end- 

1 

S/-JUN-//1: 

end- 

1 

s/-JUL-//l: 

end- 

1 

s/-AUG-//l: 

end- 

1 

S/-SEP-//1: 

end- 

1 

S/-0CT-//1: 

end- 

1 

s/-NOV-//l: 

end- 

1 

S/-DEC-//1: 
i 

end- 

1 

s/://l:end- 

1 

s/://l:end- 
i 

1 

resequence 

exit 

EDMAPS . WAY 


!  FILENAME  =  EDMAPS. WAY 
1 

t  'N=' 
delete 

move  end-1  to  . 
i 

delete  end- 2: end-1 
s/      //end- 4 : end-1 
s/    //end- 4 : end-1 
s/  //end- 4: end-1 
i 

resequence 
exit 


Appendix  B  -  Modifying  SUXES's  Strategy  File 

The  strategy  file  needs  modification  if:  (1)  SUXES's  predicted  I-V  characteristics  do  not 
resemble  the  measured  I-V  characteristics,  (2)  the  user  does  not  accept  the  percent  devi- 
ation between  these  curves,  and  (3)  the  user  is  not  satisfied  with  the  percent  difference 
between  the  simulated  and  measured  ring  oscillator  frequencies.  One  such  modification 
technique  is  to  determine  the  parameter(s)  which  is(are)  causing  the  deviation  in  the  I-V 
characteristics.  To  do  this,  plot  the  measured  I-V  data  along  with  the  I-V  characteris- 
tics using  the  initial  estimates  for  the  parameter  values.  If  this  plot  is  not  qualitatively 
satisfactory  to  the  user,  the  values  for  the  initial  estimates  must  be  reexamined.  If  using 
the  initial  estimates  for  the  parameter  values  gives  I-V  characteristics  which  resemble  the 
measured  I-V  characteristics,  then  plot  the  I-V  characteristics  after  the  values  of  the  first 
iteration  of  parameters  are  optimized. 

If  the  user  accepts  these  plots,  continue  until  an  iteration  of  parameters  fails  to  produce 
an  acceptable  plot.  In  the  set  of  parameters  which  produces  the  unsatisfactory  plot, 
eliminate  one  parameter  at  a  time  in  the  strategy  file  to  determine  which  one(s)  cause(s) 
the  deviation  between  the  measured  I-V  data  and  the  simulated  I-V  data.  By  using  this 
method,  the  parameter(s)  causing  the  deviation  can  be  found.  The  parameter  KAPPA  (the 
saturation  field  factor)  is  an  example  of  the  above.  Assume  the  user  is  not  satisfied  with  the 
upward  bending  of  the  I-V  characteristics  in  the  saturation  region,  and  assume  KAPPA  is 
causing  this  bending  after  the  first  iteration,  while  the  I  V  characteristics  using  the  initial 
estimates  do  not  have  this  trait.  A  technique  to  eliminate  the  effect  of  KAPPA  is  to  obtain 
an  accurate  initial  estimate  of  the  value  of  this  parameter  and  then  avoid  optimizing  this 
parameter's  value  until  the  values  of  the  other  parameters  have  been  optimized.  To  obtain 
an  initial  estimate  for  the  value  of  KAPPA,  a  parameter  value  must  be  optimized  which 
has  very  little  influence  on  the  result  (for  example,  the  value  of  NFS).  Therefore,  using 
about  four  or  five  different  initial  estimates  for  the  value  of  KAPPA,  SUXES  can  be  run 
to  optimize  only  the  value  of  NFS  while  using  the  initial  estimates  for  the  values  of  the 
other  parameters.  An  appropriate  value  of  KAPPA  is  found  when  the  user  is  satisfied  with 
the  resemblance  of  the  simulated  I-V  characteristics  to  the  measured  I-V  data.  Using  the 
technique  described  above,  the  simulated  I-V  characteristics  will  have  a  realistic  shape  even 
though  the  rms  deviation  might  be  lower  with  other,  less  realistic  optimization  techniques. 

After  a  reasonable  initial  estimate  for  the  value  of  KAPPA  is  found,  an  appropriate  strat- 
egy file  must  be  found.  SUXES  can  operate  in  an  automated  or  a  manual  mode.  Up 
until  now  it  did  not  matter  in  which  mode  it  was  operating.  At  this  point,  it  is  best 
to  use  the  manual  mode  which  includes  a  sensitivity  analysis.  To  obtain  an  appropriate 
strategy  file,  several  parameters  (perhaps  three)  can  be  randomly  selected  and  optimized. 
These  are  potential  candidates  for  the  first  iteration  of  parameters.  These  parameters  can 
be  optimized  as  a  group  if  no  errors  occur  and  high  values  in  the  sensitivity  analysis  are 
obtained.  Low  values  in  the  sensitivity  analysis  indicate  the  existence  of  redundant  pa- 
rameters, and  the  directions  of  the  eigenvectors  indicate  which  parameters  are  redundant 
[1,16].  The  boundary  values  may  be  altered  to  decrease  the  number  of  errors  due  to  the 
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pegging  of  a  parameter  at  a  boundary.  If  the  above  criteria  have  been  met,  and  if  a  plot  of 
the  measured  I-V  data  resembles  the  simulated  I-V  data,  then  the  first  three  parameters 
have  been  found. 

After  the  values  of  the  first  three  parameters  are  found,  the  process  is  repeated  for  the 
second,  third,  etc.,  iteration  of  parameters  until  the  values  of  all  the  parameters  are  opti- 
mized and  the  user  is  satisfied  with  the  rms  deviation.  The  values  of  the  parameters  can 
be  optimized  more  than  once. 

The  optimization  technique  for  the  p-channel  MOSFETs,  which  follows,  was  found  in 
the  manner  described  above.  This  optimization  technique  is  also  used  for  the  n-channel 
MOSFETs. 


(i)  DELTA,  ETA,  THETA  (the  first  iteration  parameters) 

(ii)  TOX,  VTO,  GAMMA,  DELTA,  NFS  (the  second  iteration  parameters, 
etc.) 

(iii)  LD,  NSUB,  KP,  XJ 

(iv)  PHI,  UO,  VMAX,  KAPPA 

(v)  TOX,  VTO,  GAMMA,  DELTA,  NFS 
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Appendix  C  -  The  Procedure  to  Run  KEYS 


The  computer  procedure  KEYS  was  written  on  a  VAX-11/780  using  VAX/VMS.  Data  are 
initially  received  from  the  ACCUTEST,  then  transferred  to  the  VAX  subdirectory  which 
contains  the  computer  procedure  KEYS. 

The  ACCUTEST  data  files  have  a  format  similar  to  that  shown  in  table  C-l.  On  the  line 
beginning  with  POS,  the  chip  position  is  recorded.  The  following  two  lines  are  ignored. 
The  data  values  are  recorded  on  the  remaining  lines  starting  with  VAL.  First,  the  Vds 
value  is  listed,  followed  on  consecutive  lines  by  the  Vgs  value,  the  Vb8  value,  and  then  the 
Ids  value.  This  sequence  is  repeated  for  the  remainder  of  the  data  sets.  KEYS  converts 
these  data  into  a  format,  as  shown  in  table  C-2,  suitable  for  SUXES.  Also,  KEYS  provides 
the  plotting  package,  PLOTMY,  with  the  Vds  and  Ids  data  points  needed  to  plot  the  I-V 
characteristics.  The  format  for  the  plot  file  is  similar  to  that  given  in  table  C-3  (without 
the  date  and  time).  The  other  plotting  package,  S5PLOT,  plots  the  ring  oscillator  node 
potential  as  a  function  of  time  and  records  the  period  and  frequency.  It  is  provided  data 
similar  to  those  shown  in  table  C-3. 

The  files  listed  in  the  order  of  their  occurrence  for  the  computer  procedure  KEYS  are  given 
in  table  C-4.  The  files  that  may  be  edited  before  KEYS  is  run  are  given  in  table  C-5.  After 
editing,  KEYS  is  evoked  by  typing  " submit /noprinter /notify  KEYS." 
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Table  C-l.       Example  of  a  Data  File  from  the  ACCUTEST 


TPN  COIN 
STA  1 

DAT  19-NOV-90 

BTI  11:11:10 

ETI 

OPE 

LOT 

DEV 

WFT 

TSN 

DSN  25 
WFN  1 

DIE  25 

POS  8.000000E+00  5.000000E+00 

VAL  1  0.000000E+00 

VAL  2  1.491000E-10 

VAL  1  1.000000E-01 

VAL  2  5.000000E-01 

VAL  3  0.000000E+00 

VAL  4  1.491000E-10 

VAL  1  2.000000E-01 

VAL  2  5.000000E-01 

VAL  3  0.000000E+00 

VAL  4  1.609000E-10 

VAL  1  3.000000E-01 

VAL  2  5.000000E-01 

VAL  3  0.000000E+00 

VAL  4  1.681000E-10 

VAL  1  5.000000E-01 

VAL  2  5.000000E-01 

VAL  3  0.000000E+00 

VAL  4  1.806000E-10 


VAL  1  5.000000E+00 

VAL  2  5.000000E+00 

VAL  3  0.000000E+00 

VAL  4  1.849000E-04 

BIN  0 
EOF 
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Table  C-2.       Example  of  a  SUXES  Data  File 


Data  from  the  ACCUTEST 


*  VDS 

VGS 

VBS 

IDS 

Width/Length= 

:  2.44/  2.44 

0.10 

0.50 

0.00 

1.49100E-10 

0.20 

0.50 

0.00 

1.60900E-10 

0.30 

0.50 

0.00 

1.68100E-10 

0.10 

1.00 

0.00 

6.12000E-07 

0.20 

1.00 

0.00 

9.11000E-07 

0.30 

1.00 

0.00 

1.03100E-06 

0.10 

1.50 

0.00 

2.23000E-06 

0.20 

1.50 

0.00 

4.04000E-06 

0.30 

1.50 

0.00 

5.47000E-06 

0.10 

2.00 

0.00 

3.82000E-06 

0.20 

2.00 

0.00 

7.24000E-06 

0.30 

2.00 

0.00 

1.02300E-05 

0.10 

2.50 

0.00 

5.34000E-06 

0.20 

2.50 

0.00 

1.02700E-05 

0.30 

2.50 

0.00 

1.47600E-05 

0.10 

3.00 

0.00 

6.78000E-06 

0.20 

3.00 

0.00 

1.31700E-05 

0.30 

3.00 

0.00 

1.91300E-05 

0.10 

3.50 

0.00 

8.13000E-06 

0.20 

3.50 

0.00 

1.59000E-05 

0.30 

3.50 

0.00 

2.33000E-05 

0.10 

4.00 

0.00 

9.45000E-06 

0.20 

4.00 

0.00 

1.85900E-05 

0.30 

4.00 

0.00 

2.73000E-05 

0.10 

4.50 

0.00 

1.07500E-05 

0.20 

4.50 

0.00 

2.11000E-05 

0.30 

4.50 

0.00 

3.11000E-05 

0.10 

5.00 

0.00 

1.19700E-05 

0.20 

5.00 

0.00 

2.36000E-05 

0.30 

5.00 

0.00 

3.48000E-05 
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Table  C-3.       Example  of  a  Plot  File 


19-NOV-90  05:15:55 
201 


0.000E+00 

1.158E+00 

1.000E- 

09 

6.697E- 

-01 

2.000E- 

-09 

5.665E- 

-01 

3.000E- 

09 

4.710E- 

-01 

4.000E- 

-09 

3.861E- 

-01 

5.000E- 

-09 

3.147E- 

01 

6.000E- 

09 

2.560E- 

-01 

7.000E- 

09 

2.082E- 

01 

8.000E- 

-09 

1.683E- 

-01 

9.000E- 

-09 

1.362E- 

-01 

1.000E- 

-08 

1.101E- 

01 

1.100E- 

-08 

8.813E- 

■02 

1.200E- 

-08 

7.098E- 

02 

1.300E- 

08 

5.754E- 

02 

1.400E- 

08 

4.680E- 

02 

1.500E- 

■08 

3.767E- 

02 

1.900E-07  1.040E+00 

1.910E-07  2.005E+00 

1.920E-07  3.091E+00 

1.930E-07  4.029E+00 

1.940E-07  4.770E+00 

1.950E-07  4.915E+00 

1.960E-07  4.977E+00 

1.970E-07  4.993E+00 

1.980E-07  4.997E+00 

1.990E-07  4.999E+00 

2.000E-07  5.000E+00 
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Table  C-4.  The  Files  Listed  in  the  Order  of  Their  Occurrence  in  the  Computer 
Procedure  KEYS.  This  list  does  not  include  files  that  are  created  during  the  program. 
The  FORTRAN  files  are  provided  data  files  which  can  be  edited  by  the  editor  files.  The 
command  procedures  are  connected  together  creating  the  computer  procedure  KEYS. 


POSSIBLE 

COMMAND  FORTRAN  DATA  EDITOR 

PROCEDURES  FILES  FILES  FILES 


(1)  KEYS.COM 

(1)  ED  NO.  WAY 

+  (1)  NLW.FIL 

+  (2)  RC.FIL 

+  (3)  LINE.FIL 

+  (4)  LOTNUMBER.FIL 

(1)  RCKEY.FOR 


(2)  LAST3.COM 


(3)  INKEY.XXX 

(5)  W*R*C*VB.DAT 

(6)  NW*R*C*.DAT 

(7)  PW*R*C*.DAT 
+  (8)  NPARIN.FIL 

+  (9)  PPARIN.FIL 

(2)  SUKEY.FOR 


'10) 

i  NINPUT.XXX 

i  PINPUT.XXX 

+ 1 

I  NOPT.FIL 

+ 1 

13 

l  POPT.FIL 

+ 1 

14 

I  NSTRAT.FIL 

+ 1 

[is; 

1  PSTRAT.FIL 

2)  EDSUXN.WAY 

3)  EDSUXP.WAY 


(4)  EDPAR.WAY 
(4)  CAPKEY.FOR 

+  (16)  IVCHARN.XXX 
+  (17)  IVCHARP.XXX 
+  (18)  NMODEL.XXX 
+  (19)  PMODEL.XXX 


(3A)  SUXES.XXX 

*  (3)  SSUXES.FOR 


(3B)  SPICEIV.XXX  (5)  EDSPICEN.WAY 

(6)  EDSPICEP.WAY 

(5)  SPICENBS.FOR 

(6)  IVKEY.FOR 

(7)  PLOTMY.FOR 


+  (20)  SUBCKT.XXX 
+  (21)  CIRCUIT.FIL 
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8 


8)  S5PLOT.FOR 

9)  MCVAX.FOR 


(4)  RMSKEY.COM 


(10)  SUMKEY.FOR 


+  (22)  PCHIPS.FIL 


7)  EDLAST.WAY 

8)  EDNRMS.WAY 

9)  EDPRMS.WAY 

10)  ED  RES.  WAY 


(5)  STKEY.COM 


83 


CRDB.FOR 
STAT2.FOR 


+  (23)  STEND.FIL     (11)  EDSKEY.WAY 


(6)  CORKEY.XXX 


(13)  RANKEY.FOR 


(24)  COREF.FIL 


S3 


ED  MAP.  WAY 
EDSIMF.WAY 


14)  EDPDIFF.WAY 

15)  EDMEAF.WAY 


(7)  MAPKEY.XXX 


(16)  EDMAPS.WAY 


*  To  obtain  this  program,  consult  K.  Doganis.  Refer  to  table  1  for  details. 
+      The  user  may  edit  these  files  according  to  table  C-5. 
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Table  C-5.       The  Files  to  be  Edited  at  the  Discretion  of  the  User 


FILENAME 


POSSIBLE  MODIFICATIONS 


FORMAT 


(1)  NLW.FIL 


(2)  RC.FIL 

(3)  LINE.FIL 


The  wafer  letter,  the  row,  and  the  column  of  the  first  (Al) 
chip  to  consider  are  entered  into  this  file  on  the  first  (12 
three  lines.  If  the  row  number  is  zero,  the  comprehensive  (12 
summary  is  started  from  scratch.  This  is  done  before 
each  lot  is  analyzed.  The  fourth  and  fifth  lines  should 
not  be  modified. 


KEYS  considers  the  chips  (denoted  by  wafer  letter,  (A1,1X,I2, 
row,  and  column)  in  this  file.  IX,  12) 

The  data  lines  of  this  file  contain  the  wafer  letter,  (2X,A1,2X, 
the  row,  the  column,  the  width  in  micrometers  of  I2,1X,I2, 
the  n-channel  MOSFETs,  the  length  of  the  1X,F5.2,1X, 
n-channel  MOSFETs,  the  width  of  the  p-channel  F5.2,1X,F5.2, 
MOSFETs,  the  length  of  the  p-channel  MOSFETs,  1X,F5.2,1X, 
the  measured  ring  oscillator  frequency  in  MHz,  F5.2,1X,I3) 
and  a  designation  indicating  NMOS  technology  (1), 
PMOS  technology  (2),  or  CMOS  technology  (3). 


(4)  LOTNUMBER.FIL  The  lotnumber  is  entered  into  this  file. 

(5)  NPARIN.FIL 


(6)  PPARIN.FIL 


(7)  NOPT.FIL 


(8)  POPT.FIL 


(9)  NSTRAT.FIL 


The  n-channel  MOSFET  oxide  thickness  in 
angstroms  is  followed  by  the  lateral  diffusion 
in  micrometers. 

The  p-channel  MOSFET  oxide  thickness  in 
angstroms  is  followed  by  the  lateral  diffusion 
in  micrometers. 

The  n-channel  MOSFET  options  file  for  SUXES. 


The  p-channel  MOSFET  options  file  for  SUXES. 


The  n-channel  MOSFET  strategy  file  for  SUXES. 


(10)  PSTRAT.FIL         The  p-channel  MOSFET  strategy  file  for  SUXES. 


(11)  IVCHARN.XXX 


For  the  n-channel  MOSFETs,  modifications  to 
the  SPICE  device  parameters  (AD,  AS,  PD,  PS, 
NRD,  and  NRS)  and  output  characteristics  are 
made  here.  The  program  will  input  the  widths 
and  lengths;  therefore,  do  not  modify  this  portion 
of  the  file. 


(2A4) 
free  format 


free  format 


see  reference 
[1] 

see  reference 
[1] 

see  reference 
[1] 

see  reference 
[1] 

see  references 
[2,6] 
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(12)  IVCHARP.XXX 


(13)  NMODEL.XXX 


(14)  PMODEL.XXX 


(15)  SUBCKT.XXX 

(16)  CIRCUIT.FIL 

(17)  PCHIPS.FIL 

(18)  STEND.FIL 


For  the  p-channel  MOSFETs,  modifications  to 
the  SPICE  device  parameters  (AD,  AS,  PD,  PS, 
NRD,  and  NRS)  and  output  characteristics  are 
made  here.  The  program  will  input  the  widths 
and  lengths;  therefore,  do  not  modify  this  portion 
of  the  file. 

For  the  n-channel  MOSFETs,  modifications  to 
the  SPICE  model  parameters  LEVEL,  CJSW, 
RSH,  JS,  TPG,  MJ,  MJSW,  PB,  and  XQC  are 
made  here.  The  program  will  input  the 
other  parameters;  therefore,  do  not  modify 
this  portion  of  the  file. 

For  the  p-channel  MOSFETs,  modifications  to 
the  SPICE  model  parameters  LEVEL,  CJSW, 
RSH,  JS,  TPG,  MJ,  MJSW,  PB,  and  XQC  are 
made  here.  The  program  will  input  the 
other  parameters;  therefore,  do  not  modify 
this  portion  of  the  file. 

Modifications  to  the  SPICE  subcircuit  file  are 
made  here.  The  program  will  input  the  widths 
and  lengths;  therefore,  do  not  modify  this 
portion  of  the  file. 

Modifications  to  the  chosen  SPICE  file  are 
made  here. 

The  percentage  of  acceptable  chips  needed  to 
make  wafer  maps  is  the  first  data  value,  and 
the  number  of  test  sites  per  wafer  is  the 
second  data  value. 

Modifications  to  the  STAT2  stend  array  are 
made  here. 


see  references 
[2,6] 


see  references 
[2,6] 


see  references 
[2,6] 


see  references 
[2,6] 


see  references 
[2,6] 

free  format 


(1X,4(I2,1X)) 


f  :. 
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J3. 


4.    TITLE  AND  SUBTITLE 

Evaluating  a  Chip,  Wafer,  or  Lot  Using  SUXES,  SPICE,  and  STAT2 

5.  AUTHOR(S) 

J.  C.  Marshall  and  R.  L.  Mattis 
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